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ON THE ANALOGUES RELATING FLEXURE AND 
EXTENSION OF FLAT PLATES 


By R. V. SOUTHWELL (Ozford) 


teceived 10 January 1950] 


SUMMARY 


The displacement of a flat plate bent by transverse loading, and the extensional 
(Airy) stress-function in ‘plane strain’ or in ‘plane stress’, are governed by equations 
of identical form; and the boundary conditions have identical form when edge- 
displacements are specified in the flexural, edge-tractions in the extensional problem, 
so mathematically, in these circumstances, only a single problem is presented. This, 
the ‘first analogue’ relating flexure and extension, is well known. 

A ‘second analogue’, relating the flexural problem when edge-tractions with the 
extensional problem when edge-displacements are specified, is believed to have been 
first propounded in 1941. By introducing two quantities U and V, analogous with 
the components uw and v of extensional displacement, it permits a treatment of the 
flexural problem by any method—e.g. ‘two-diagram (relaxational) technique’ 
which yields extensional solutions of this second type. 

In this paper both analogues are combined in an inclusive statement covering the 
perforated (multiply connected) plates which were discussed in 1948. Reasons are 
stated for believing that ‘two-diagram technique’ is preferable in problems governed 


by ‘mixed’ boundary conditions. 


Nomenclature 
2h denotes the (uniform) thickness of the plate. 
p denotes the uniform density of its material. 
E denotes the Young’s Modulus of its material. 
uw denotes the Modulus of Rigidity of its material. 
o denotes the Poisson’s Ratio of its material. 
v denotes the direction of the outward normal to an edge. 
s denotes distance measured along the edge. 
(x,v) denotes the angle (measured counter-clockwise) between the 


x-direction and the direction v. 


[7 cos(x, v), sin(x, v).| 
Symbols Spe cial to the fle xural proble m: 


9 


-o* 


2 Eh’ 


D denotes the ‘flexural rigidity’ of plate | - *s 


w the transverse (flexural) displacement, 

Z the surface-intensity of transversé loading, 
¢@ a function such that DV74 Z, 

[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 3 (1950)] 
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} 
N,., N, the line-intensities of shearing action on sections perpendicular} 2, Th 
to Ox. Oy, hending 
M,, M, the line-intensities of flexural couple on sections perpendicular 
to Ox, Oy, 
H_,,, H, the line-intensities of torsional couple on sections perpendicular 
to Ox, Oy I 
[ H, H, H, say (cf. section 6)], 
M, the line-intensity of flexural couple on an edge, and the 
H, the line-intensity of torsional couple on an edge. equilibr 
N.,N,) 
[(/H-+mM, = 1H,+mM,, IM,+mH = 1M,—mH,.| 
U,V are ‘quasi-displacements’ defined in section 3. ‘ 
in whiel 


Symbols special to the extensional problem: 
A, k denote functions of o depending on the state of stress, 


u, v the component displacements of the middle surface in the directions 








Ox, Oy ia j 
; . C ov | 
A , 
Cx cy I 
© the potential of the conservative body-forees X, Y, 
: eQ P . ’ , : 
px p the component of body-force in the direction of x, 
Cox 
: 0Q , : : ; ’ : 
pk p the component of body-force in the direction of y, 
cy 


x, Airy’s stress-function, 
X,, X,, Y, the components of stress in the plane of the middle-surface, 


Cnvs Crys Cy, the components of strain in the plane of the middle-surface, 
X,, ¥, the components along Ox, Oy of tractions applied to the edge. 


[X, = 1X,+mX,, Y, = 1X,+mY,.] 


L V 


Introduction 
1. By seeking to adapt the relaxational technique to computation of | 
flexural and extensional stresses in flat plates, I came to consider (1938) 
the mathematical analogues which relate these physically independent | 
problems. One was well known though seldom utilized, e.g. to solve | 
flexural problems by photo-elastic methods; but a second (so far as I know) 
was first propounded in a paper written jointly with L. Fox (1). Here I 
attempt to present both in one comprehensive statement, and thereby to 
show that ‘slices’ and ‘slabs’—that is, the extension and flexure of flat 
plates—are problems which should not be treated separately (as they , 


commonly are) because mathematically they are identical. 
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cular} 2, The flexural displacement w is related with the line-intensities of 
ending couple (/,,.M,) and of torsional couple (H, H H) by* 


Jf y 
cular { ¢ TB Vaid Thi C~-w 
M D tg .)> H = D(l—o) ; 
Ox cy- CXLCY 
cular | ~ fal (1) 
M, = D(—— +0 3) 
cy” Cr 


nd these expressions substituted in the three equations of momental 
quilibrium yield expressions for the line-intensities of shearing action 


V..N.). together with the relation 


| DV4w = Z, (2) 
| n which Z stands for the surface-intensity of transverse loading. Z being 
specified, (2) governs w at points in the middle-surface of the plate; and 
| 1 Song . . 
vhen the boundary conditions relate to displacement, then w and its 


normal gradient éw/év have to take given values at the boundary. 


ions pa nik 
| The extensional stresses, expressed in terms of “Airy’s stress-function’ y 
Dy ; O4 ; oy ; C7y . 
| X rt. 9Q) Xx x :. = x L pQ, (3) 
Cy CXC Y CL” 
} . . . . . . . 
satisfy the two stress-equations of equilibrium which hold in circumstances 
her of ‘plane stress’ (Z Z Z. 0) or of ‘plane strain’ 
(Z Z é 0): 
| the conditions for compatibility of the accompanying strains require 
UX shall satisty the equation 
‘ace \ ty kV *(pQ2) 0, (4) 
200 here ; . 3 
-_ i l—o in circumstances of ‘plane stress | . 
se, ae (9) 
h (1 —2oe)/(1—c) in circumstances of ‘plane strain’ } 
u 1 1 . . . 
| when the body-forces are conservative and hence expressible in terms 
a potential Q (ef., e.g., (3), sections 401-6). When the tractions 
1 of | *Pplied to the edge are specified, boundary values of Gy/éx, €x/cy can be 
L . 
20 omputed from the expressions 
» } 
lent CX ‘ , :, CX ' 
xX fY, — pQ sin(a, v)! ds, x {X,,—pQ cos(x,v)} ds, (6) 
Ive co] , cy p 
Ow y and its normal gradient @y/év have to take given values on the 
re | , boundary.t Accordingly (since the forms of (2) and (4) are exactly similar), 
vy tO | that case of the flexural problem in which edge-displacements, and that case 
flat {the extensional problem in which edge-tractions are specified, reduce to one 
hey , 


Cf. « 2), section 298. This paper employs a slightly different notation. 


Cf. (1), 1941, section 7 
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and the same mathematical problem. This is the first of the analogues Clear 
mentioned in section 1. in (7), 1! 
in effec 


3. When edge-displacements are specified in an extensional problem,} shall fi 
use of the function yx is less convenient (at least in a relaxational treatment)| stability 
because its second derivatives then enter into the boundary conditions; and] case ‘B 


a better procedure is to work in terms of component displacements w and] in case 


v, related by the equations 4, In 
oA , : interna! 
dx A(V*u+ pX/p) = 0, entail n 
\}by a d 
oA ms ; . 
— + A(V7v-+ pY/p) = 0, | exempl: 
cy | 
Lang a 
. ‘ ; ou ov ie 
in which A stands for (‘ 1. and applied 
ae yy method 
—— ' ; hem i 
A (1—o)/(1+-c) in circumstances of ‘plane stress’ | | — = 
Seti: ; oe ‘} been co 
A (1—2o) in circumstances of ‘plane strain’. | 


duced ¢ 


The second analogue mentioned in section 1 relates this case of the{ graph | 








extensional problem (case ‘A’: edge displacements specified) with that} traction 


case (“B’) of the flexural problem in which the line-intensities of force and} in the ¢ 


couple (i.e. tractions) have to take given values at the boundary. It makes} but ‘Sa 


use of two quantities U and V which at all points in the middle surface! dence i 
} 


of the plate have to satisfy+ 5. N 
‘ ' : bounda 
oA l-+-oc _,,;; od 
on ' 1 ——" > 0 points 
er (9\| fronted 
oa i "v{V4 9 = which 
cy l—o cy 


relative 


s : ol oV ‘ ' P wie oe be clait 
when A’ stands for | and where ¢ is a calculable function. Theit| 


Ox ey plates 
boundary values are to be calculated from  Tegards 
. , 2 | standin 
c , ~ CD 
D = (M,— D¢)sin(x, v)—cos(x, v) (x, . fo ls { Sut 1 
cs CS 1 . 
‘ | evident 
oV : ‘ie O02 Oq | (sectio 
D (M,— D¢)cos(x, v)+-sin(x, v) | N— —_"+D p ds. a4 
os Z cs CV Moreoy 
(10)| general 
Cf. (4), equations (25)-(27)—which also carry numbers (with asterisks) referring t 
their derivation in (1). The quantities V2w,, A, G, H, N, in those equations, are her + Th 
, . . © 
replaced by ¢, A’, —M,, H,, N, respectively, and ¢ is not assumed to be zero on thi by L. Fo 


boundary. (Cf. sections 7, 8 and 11.) 
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Clearly, (7) and (9) are exactly analogous except that the constant A, 
. (7), is replaced by a different constant in (9); and this change amounts 
n effect) to substitution of a fictitious value o’ for Poisson’s ratio o. We 
shall find (section 10) that the substituted value implies mechanical 





tability of the presumed material; so nothing obstructs a solution of 


ndl case ‘B’ of the flexural problem by any method that has been successful 


n ease ‘A’ of the extensional problem. 

4. In 1948 I extended this ‘treatment by analogy’ to plates which have 
nternal boundaries (holes) either loaded or unloaded. I showed that these 
ntail no serious complication of the relaxational technique, being tractable 
yy a device (‘block relaxation’) which has long been customary; and 
‘xemplified the requisite computations by a solution, obtained by Miss E. D. 
Lang and Miss G. Vaisey, for a tension-member (‘tie’) to which load is 
ipplied through ‘pins’. The ‘isochromatic lines’ of the photo-elastic 
method were deduced and presented in a diagram ((4), Fig. 7) which gave 


them in somewhat insufficient detail. More recently further ‘lines’ have 





f the| 
_ 
tha 
> ant 
lake 


rfac 


“hei | 





been constructed by Miss Vaisey, whose more complete diagram is repro- 
duced as one part of Plate I. For the accompanying photo-elastic photo- 
graph | am indebted to Dr. M. M. Frocht: the uniform longitudinal 
traction on a half-circumference, which was the mode of loading presumed 

the computations (Plate Ia), can hardly have been exactly reproduced ; 
but ‘Saint-Venant’s principle’ justifies a comparison, and the correspon- 
lence is seen to be very close. 

5. Neither in 1941 nor 1948 was detailed consideration given to ‘mixed’ 

undary conditions (specifying displacements at some, tractions at other 
points on the edge), because no practically important case had been con- 
fronted. Two distinct relaxational techniques were available, either of 
hich could deal with ‘mixed’ conditions,+ and the question of their 
relative merits was “left to the test of experience’. But by 1948 it could 
be claimed that with their aid any problem concerned with flat elastic 
plates is soluble with sufficient accuracy for practical purposes; and as 
regards flexure this position would not have been attained without under- 
standing of the ‘second analogue’ {section 3). 

But the method employed to derive this second analogue was not made 

ident in the papers of 1941 and 1948; the significance of the function ¢ 
section 3) is not apparent, and the derivation of (10) seems artificial. 
Moreover, the two analogues are related: they can be combined in a single 


jeneralization. This, the main thesis of the present paper, is established 


The application of one of these (\ two-diagram technique’) has been exemplified 
L. Fox (5). 
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in sections 7-9. It is utilized in discussions of ‘mixed’ boundary conditions} underste 
(sections 11—12) and of ‘perforated’ plates (sections 13-14). and we § 

My thanks are due to Miss G. Vaisey for her assistance, and to the} substitu 
Clothworkers’ Company whose grant made that assistance possible. I a 
indebted to her both for the drawing in Plate I and for the computatio 
and plotting which this entailed. (The earlier diagram mentioned i 
section 4 was drawn by Mr. D. N. de G. Allen.) 


First extensional-flexural analogue (‘Case A’ of flexure and 
‘Case B’ of extension)+ 
6. The first analogue used in the paper of 1941 is based on the similarity | and (13 
of equations (2) and (4). Of these the former is derived from the conditions 
of flexural equilibrium, viz. 


agret 
Z+- Me, 2M, 0 The 
= oy conditi 
v oM, _é H, ‘ a siderati 
Ox cy analog 
Nv oH, eM, o 
Cx ey 


(in which H, = H, = H, say), by substituting for M,. M,, H from (1)., “he ¢ 


uv 





Thereby we obtain the expressions 


Secon 
; ¢ ; : é ; ‘Ca 
N, D — (Vw), N,, D — (Vw), (12 a 
Cox cy /. \ 
: . . : . (sectic 
and these substituted in the first of (11) vield (2). 
The second (extensional) equation (4) is obtained from the condition 
of compatibility , 9 and b 
= Cc” oa Sa ‘ . : 
5 Ore Cry Cy, = 9 (13) } elimin 
oy™ CLCY ox | 
by substituting for e,.., @,.,, €,,, from 
] - : 2(1+<c) . 
e.. i [X,—o(Y,+Z.)], ete., ts K } (14)| ande 
and then for X,, X,, Y,, their expressions (3). Two values (5) are taken 
by the constant & in (4) because Z, = 0 in circumstances of ‘plane stress’, | 
é.. = 0 in circumstances of ‘plane strain’. 


It can be shown (ef. (8) of section 3) that the difference between plane ) 
stress and plane strain amounts in effect to an alteration in o. On that 


The different ‘cases’ were thus denoted in (1) (section 8), but were given different letters ) 


in (4). 








itio 


and 


arit 


tio} 
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understanding only one sort of extensional system need be considered, 
nd we shall concentrate conditions of plane stress. Here (Z, being zero) 


3) for X,, X,, J 


Cll 


substitution from transforms the relations (14) to 


7 o- 4 
Ke V°x—(1-+c) x (1—o)p22, 
Ca 
Ki %(1+c)—*, (15) 
eOY 
Ke v2) o) X + (1—o)pQ, 
cy~ 
ind (13), when these expressions are substituted, becomes 


V2V2, (16) 


wreeing with (4) and (5) in circumstances of plane stress. 


The 


nditions are also similar (as they 


formal similarity of (2) and (16) is evident, so if the boundary 
are in the two cases now under con- 


sideration) the flexural and extensional problems are mathematically 


nalogous. Exact correspondence will obtain when 
Du Z - 
— const. (17) 
x | r)V*(po2) 


[he constant is not necessarily ‘non-dimensional’.) 

Second extensional-flexurai analogue (‘Case B’ of flexure and 
‘Case A’ of extension)‘ 
7. We proceed te 


» explain the derivation of the ‘second analogue’ 


section 3). Let d be a function defined by 
DYV*d Z (18) 
| by a boundary condition to be stated later. Then, when N,, N, are 
minated from (11), there results 
c*H c 
.(M.— Do) +-2 (VW, —D¢) 0, (19) 
CL CHCY cy” 
nd comparing this with (13) we observe that 
' ° cu 
(M Dd) corresponds with e.., 
: - OX 
} ° cl 
W_.—D¢) corresponds with e,,, : (20) 
oy 
P cov cu 
2H corresponds with e,, }-—. 
, CX cy 


Cf. footnote to 


section 6. 
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The correspondences suggest that we express (M,—D¢?), (M,—D¢), —2H 
eU eV V oe , j ; 
by pb. p<‘ ; vf } c respectively, thereby introducing tw 
Cx cy Cx cy 


quantities U and V which are the analogues of uw and v. In this way 
giving their expressions (1) to M,, M,, H, we obtain relations as unde 
in U and V: 


oU ctw o7u 
- +o ——¢ 
Cx ey? Cx? 
ove ew 
—. +g—.— @, (21 
ey Ox? cy? | 
eV. aU ow 
| -2(1—o) —_. 
Cx (Oy Cxoy 
8. The most general solution of the third of (21) is 
: Cw. Ob 
l (l1—o) } ’ 
CX CX 
22) 
; cw OCs 
| (l1—o) : 
cy cy 


% denoting some function of x and y as yet undefined; and when these 
expressions are substituted'in the first and second of (21) we have 
Cubs C7us 


3 3 
cy- Ox" 


V2w—d, (23) 


whence V2 


two equations and (22). 


0 and so (9) result when w is eliminated between those 
Also, from (22) again, 


wit h 


C 
So D ; 
( 


und simi 


D. 


But e 


Conseql 
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D 


and the 

As yt 
the bou 
can be 
) both d 
Comb’ 

9. V 


and CC 


and v 





U 52 
: (: : cos. v) : sin(a, v) ; | 2 ee ") 
C8 \ CX cy cy Ox CXC Y 
c . c m 
2| cos(r,») L sin(a, v) (vee ¢), by (23), 
Ox oy ‘ 
9° (V2w—¢). (24) 
Ov 
Now N, cos(x, v).V,.+-sin(x, v)N, D— Vw, by (12). 
CV 
Hence and from (24) we have 
D(“ 2 | LN, | pi) ds, 
cy Cx ‘ CV 
and from this and the expressions in section 7, 
Vv : ; y vs 
D! | (. 1 D) ds—H, p‘ | (X.+ Db’) ds+-H, 
oy J Cv Ox 2 Cv 





a term 


section 
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eu OV 
vith D M,—D¢, p* M,—D¢. 
Cx . cy 
oU U U 
So D D! cos(x,v) { sin(a, v) < 
ak | oy ox | 
estes: Gal [ (x D) as H\—sin(x,v)(M,—D¢) 
ae a o1 , jai i . 
ind similarly » (25) 
Do. see aie ~ Rd sin(v, ») | (%, pe) ds H). 
cs CV 


But equilibrium of an element at the boundary requires that 


cos(x,v)H+sin(x,v)M, cos(x, v)H,+-sin(a, v)M,, | 


(26) 
cos(x,v).M,-+sin(2, v)H cos(a, v)M, sin(a,v)H,. } 
( onsequently the relations (25) are equivalent to 
U iz | C . 
ps cos(a. v) (N, Hl, pt) ds—sin(x,v)(M,— Dd), 
cs cs CV 
. (27) 


pe! cos(x, v)(M,—D¢)+sin(a, v) | (, oH, } p<) ds, 


cs J cs CV 
nd these yield boundary values of U and V.+ 
\s yet no boundary condition has been imposed on ¢. If its values at 
the boundary are specified, its values elsewhere are determined by (18) and 
an be found (approximately) with the aid of Relaxation Methods; then 
both d and éd/év will be known 
Combination of the two analogues 


9. We have from (15) (in circumstances of ‘plane stress’), 


Cu l | Cc “y c zy 
Cry | a x + (1—o)pQ)|, 
Ca k | cy” Cx" 
aL l | o-x C *y O 28 
Cory 4 > Co = = (1 a) p&2 ; (2! ) 
cy k | ¢ _ cy~ 
Cv Ou 2(1+o) cx 
ft . * > 5 
Ca cy iad EB Cc UCY 
ind comparing (28) with (21) we see that U and V will correspond with u 
und DI whe n 
w corresponds with y/E, 
¢ corresponds with (l—o)pQ E, (29) 
o corresponds with o. 
Equations (10) and (27) are identical. We observe that N, and H,, are associated in 
term | N on a as in Kirchhoff’s explanation of the double boundary condition (Cf. (2), 


section 297), 
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Then with the left-hand side of (16) there will correspond 
E(V4w—V"d) E(V4w—Z/D), by (18), 
that is, equation (2) will correspond with (16) and so with (4); also 
aU al 


A’ will correspond with A in equations (7), 
CX cy 


Poy , l—o\. ; ” 
will correspond with A in equation (7), 
o \ lo 


(od Ch 
2| 
Ca 


QO aQ’ To 
2 will correspond with —2(1—c) p ee al A®(xX,} ), 
cy A 


cx CY r 
(3¢ 
that is, equations (9) will correspond with (7). Thus except that th 
sign of o is reversed, an exact correspondence has been established whiel 


embraces both of the preceding analogues. 


10. It remains to show that the change in the sign of o is acceptable 
having regard to the fact (cf. (3), section 123) that a material would b 
mechanically unstable whose Poisson’s ratio had a value outside of the 
range . 

: i<¢ 4; (31 
jor uf Relaxation Methods are applied to an unstable system, spe cial devices 
must be used in order to ensure that the solution will not diverge. 

We have seen that equations (9) relate to a quasi-extensional system in 
which o is replaced by a fictitious value o’ a. Clearly, in order that 

ae 1. 


o must have a value in the range 


L Se 4. (32) 





systems 


or of Us 
at parts 
involvi 
and V 
specific 
this ac 
proble: 
similat 


derive 


For 





and Ci 


in wh 


so in order that a material which itself is stable may yield a stable material 


for the quasi-extensional system, both of the criteria (31) and (32) must 


be satisfied.—that is to say, its o must have a value in the narrower range 


5<C . (33) 
Now in fact, while (31) is the range imposed by theory, (33) is satisfied 
by all materials which have practical interest. Accordingly it imposes 


no restriction on the use that can be made of our two analogues. 


Mixed boundary conditions in the flexural problem 


11. When tractions are specified at some, displacements at other parts 


of the edge (or boundary), we have a choice between treating flexural 


In t] 
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Mix 
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systems in terms of w or of U and JV, extensional systems in terms of x 


orof wand v. Working flexural problems in terms of w we have to satisfy, 


; :; . OH i = 
at parts of the boundary where /, and [N, ‘ are specified, conditions 
Ss 


‘involving second and third derivatives of w; working them in terms of U 


: ow ow 
and V we have to satisfy, at parts of the boundary where >= are 
; ox Cy 


specified, conditions involving only first derivatives of U and V. On 


this account the second procedure would seem preferable in flexural 





r), problems where the conditions are ‘mixed’, and in extensional problems, 
(2% similarly, working in terms of displacements u and v. We proceed to 
0 2 

_ | derive the requisite conditions at a boundary. 
Db tt . =< , . 
For flexure, writing (21), section 7, in the equivalent forms 
vhie 
cru c | ol 
o- | a)h oO 
ca" cy Ca 
] C-W Cc U C | 

able, | (1—o?) l—a)d hae. (34) 

ld be | cy” Cr Cy 

r the C-u ol ol 

2 rT) ‘ 
(3] CXCOY Ca cy 
ahs ud combining these with the identities 
oj{cu = ab o~t ¢ Cw) Ow o*u = 
} ] m | l -—m ' (35) 

ml C8 \« x} CXOY Cx OS \CY cy~ CXCY 

that .: 
which 7, m stand for cos(a,v), sin(a,v), we deduce that 

( Cw C V ol 

9(] ) | cos(x.v)(1 0) 

C8 \co2 Ox Cy 

(32) | oV ol 
2 sin(a v)) o (1 o)d| 
Tia oy sa 
| (36) 
LUST ( Cu ou oV 
2(l—o | 2 cos(a.1 )! oO 1 (] o)d| 

nge } cs \cy {Cx cy 

(33 ; ov au 

sin(x,v)(1+o) { ) 

fied a 

Ses In these, the wanted boundary conditions in U and V, we may give ¢ any 

he undary values which make for simplification. 
) 
Mixed boundary conditions in the extensional problem 
rts | 12. From (36), merely by using the analogue, we can deduce correspond- 


iral ing boundary conditions for extensional systems. (As before we assume 
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conditions of plane stress, for which the correspondence is as stated a 
: _ ‘ — 0 O [c : 
(29), section 9.) Then we have expressions for - ox and _X which 
cs Cx cs oy 
X 


Cc 
when we replace - 


ox , ‘ . ' 
, * by their expressions (6), becomet 
cx oy 


1—o* ,, Cu ov) 
2 x = 2 cos(x, v) Lo } 


E Cx cy 
; ov. cu 
(3 e)sin(r,»)( +7 ) 
CX cy 
— , e (37 
i , ° v v 
2 i y. 2sin(x,»)(é _ + 
E cy CX 
ov CU 
1 (1 a)eos(r»)( -_ ) 
Cx cy 


agreeing with equations (5) of a paper by L. Fox (5) which treats of 
‘mixed’ conditions in extensional problems. 





Conditions at internal boundaries (holes) in perforated plates 


. fe Cw cw ¥ ae } 
13. From (36), because and must be single-valued, it follows that 


6x cy 
the integrals 
(017 AV “(tT av 
ey 4 eB e's 24) ds, | Bcd Ts 24) ds (38) 


| os Cv A | ov Os A 


must vanish severally in respect of every interior boundary (hole): and 
then, in order that w may be single-valued, the integral 
hd c f \ { a9 a9 
cw Cow O-w ow 
) (i m ds (x Ly ds 
A cy Ox ; OSCX CSCY 
i.e. the integral 
F oU . av ove au ma—li ies 
‘ Fe — u(- _ } aa + 24) | ds (39) 
; | Os Ov Cs Cv A } 
must vanish in respect of every hole. (Here /, m stand for cos(, v), sin(, v), 
as before. A now denotes the ratio (1-+-c)/(1—c).) 
These are the conditions stated in equations (41) of Southwell (4). They 
have no analogue in the extensional problem because (ibid. section 8) 
ox OX 


Cx Cy 


x need not be single-valued. 


+ The cancellation of the terms in 2 was to be expected, since the body-forces enter 
neither into the stress-strain relations nor into the equations for transformation of stress. 
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14. In extensional solutions the rotation and the displacements must 
be single-valued, and hence (ibid. sections 12-13) the conditions 


. 


c > 
‘ {Vx 1_(] a)pQ! ds 0 
o1 


« 


. 


¢ ly ; sin(x, v) |{V2y+(1—oe)pQ} ds O +; (40) 
CV 

p | Yr COS(2,1 ) {V2y +-(] -o) pQ} ds 0 
Ov 


must hold in respect of every internal boundary. Using the correspondence 


stated in (29). we obtain analogous flexural conditions as under: 


¢ : (V7w dh) ds 0. 
CV 
ly sin(x,»)| (Wt d)ds = 0, > (41) 
OV 


. 


- : Ccos(x, v) (vee d) ds 0. 


CV 





Relations equivalent to (41) are obtained (¢ here replacing V*w,) when 
Z, L, M are eliminated from 





Z D = V0 ds, Do ly = sin(x, »)|Wee ds, | 
M D ¢ x COs(x, Vv) Vee ds. 
Lo | 
(42) 
and re 
D | V2w, ds+Z 0 
CV 
D t y ; sin(x,») [Ve da = 3 0. (43) 
CV 





1 





cos(xv.v)—a2 : ven ds+M 0. 


which are (6) and (32) of Southwell (4). They determine the constants 
1, B, C in the expression 
db A t Bx + Cy 


vhich ¢ (= Vw,) must have at every internal boundary. When they are 


satisfied, U and V will be acyclic (ibid. section 19). 
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Isochromatic’ lines for a tension-member strained by loads applied through pins, 
s determined (a) by the photo-eclastic method in the laboratory of Dr. M. M. Frocht 
Illinois Institute of Technology) and (6) by relaxational computation (G. Vaisey). 
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THE BENDING OF A UNIFORMLY LOADED 
| CLAMPED PLATE IN THE FORM OF A 
| CIRCULAR SECTOR 

By H. R. HASSE (The University, Bristol) 
Received 21 February 1950] 


SUMMARY 


deflexion of a uniformly loaded plate in the form of a semicircle clamped 
g its boundar »btained | 1 method due to Weinstein. This problem 
s the solut of the biha nic equation AAw Z, where Z is given, 
ct to the conditions that w 0 and cw/oN 0 on the boundary, N being the 
of the outward normal. The solution is expressed in the form 
i > Ae ©, 
1 
, writing &% for W or v,,, % is found by solving (in succession) two harmonic 
ms of the forms Af Z, where Z may be zero, and Ads f, where f and & 
satisfy certain boundary conditions. The constants A,, are then determined 
tisfy the boundary condition N 0. 
rical calculations show tha or six terms of the series ¥ A,, v, give a 
n=1 
pproximation to the accurat ue as judged by the closeness with which 
pproximat olution satisfies t boundary condition cw/cN 0. The proce- 
» adopted in the ea f 1 reneral circular sector and for non-uniform 
s indicated bri 
connexion between the deflexion problem and that of plane strain in which 
tl $s function itishe the juat on AAd QO, where d and cd/cN have given 
ies on the boundary, is discussed as a preliminary to the further consideration 


tter problem by a method same type. 


luis paper deals with the deflexion of a uniformly loaded plate in the form 

. circular sector whose boundary is clamped, a problem which does not 
ppear to have been investigated hitherto. The method used is that 
eveloped by Weinstein, originally for the purpose of calculating lower 
minds for the modes of oscillation of clamped plates and applied later to 
the statical case of the deflexion of a uniform clamped square plate (1). 
lhe theoretical justification of the method is given in a further paper, in 
vhich the question of convergence and that of the connexion of the method 
vith the equivalent problem in the calculus of variations are fully 
treated (2). 


1. A brief résumé of the method is as follows. The differential equation 


lor the deflexion w of a plate of area S is 


AAw = q/D, (1) 


[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 3 (1950)] 
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The ec 


where A is Laplace’s operator, q the load per unit area, and 4 
in the e 


D 2 Eh3/3(1—o?) 


in the usual notation, where 2/ is the thickness of the plate. The & 


ditions on the boundary C are is a min 
uw 0. ow 0. yt to infini 

oN ; in succe 

where NV is the outward normal to C. tions th 


The first step is to obtain the solution of (1) when the boundary ¢q).| below. 
dition éw/eéN = 0 is replaced by the condition that Aw = 0 on C. This} relaxes 
modified problem is equivalent to solving in succession two problems, ead} of J ¢0! 
dealing with the equilibrium of a loaded membrane. Denoting the soluti: Fina! 
by W and putting AW = f, the first of the above problems is that of find.| that 
ing f when Af = q’ in S, f = 0 on C, where q’ = q/D, and the second is t 
find W when AW = fin S, W = 0on C. ee 

which, 


To solve the actual problem as specified in (1) and (2) we write as th 








mth approximation to the solution 


e W+ 2 Anr,, \) for all 
n 0 
where A,, are constants, so that Usir 
“ m 
Aw = AW+ > A, >p,, where Av, = p, 
n=0 
and, since 
m so tha 
AAw = AAW+ > A, Ap,, and AAw = AAW = q/, 
n=0 for all 
we take Ap, = 0 so that p,, is a harmonic function in S. Again AAv, = in the 
and, since w = 0 and W = 0 on C, v, = 0 on C. 
It is now necessary to introduce the boundary condition éw/éN = 0 on| - , 
; : ; ae Se ; n de’ 
C, which is done by writing it in the form 
meas 
cow nin 
p, —- ds 0 (4)\  expre 
k oN , 
co W 
for all n from 0 to m. This is equivalent to ) wher 
|| p, Aw dS = 0, (5)! FB 
‘'s 
by using Green’s theorem in the form, for two dimensions, This 
by 1 
. r Cys C = 
| | (dAy—yag) as = | (5. Fede. thre 
33 } oN oN as ir 
S 94 j i 
mus 


putting ¢ = p,, 4 = w, and noting that Ap, = 0 in S and w = 0 on C. 


508 
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The condition (4) is a weaker condition than éw/eN = 0 on C, so that 
n the equivalent problem of the calculus of variations, viz. that 





J ( ( (Aw)? 2q'w| ds, 
1e eC ‘Ss 


sa minimum for the exact solution obtained from (3) by making m tend 


to infinity, the minimum values of J obtained by taking m = 9, 1, 2, 3.... 


n succession will steadily increase. Thus in these successive approxima- 


| tions the minimum value of J for the exact solution is approached from 


v cop below. If the problem were treated by the Rayleigh—Ritz method, which 
TI | relaxes the constraints as the approximation proceeds, the minimum value 
8, each} of J corresponding to the exact solution would be approached from above. 
lutic Finally the values of the constants A, are obtained from (3) and (5) so 
f fin that dt i m 
| p, (AN > A,p,) dS 0, 
d IS ft “8 s=0 
_| which, using the usual notation || PQ dS = (P,Q), is written in the form 
as tI S 
(p,, AW) > A (py. P,) 0 (6) 
Ss 0 
for all values of n from 0 to m. 
Using Green’s theorem again, putting ¢ v,, % = AW, it follows that 
|| p, AW dS = || vo, q' dS, 
‘s ‘s 
» that (6) becomes (v5  ) > A (Pes Pp) 0 (7) 
for all xn from 0 to m. This gives m equations for the A,, (m OG: &, 2.3.3 
, =" inthe successive approximations to the exact solution. 
eu 2. The case of a semicircular sector of unit radius will now be considered 
n detail. If the angle @ is measured from the axis of symmetry and r is 
measured from the centre of the straight part of the boundary, the 
(4) expression for W is, taking q l, 
HW ar ia! feos 26 192 eos 46 > F, rk 2cos ké > Gr COS ké, (8) 
where k 1. 3, 5...., and 
(5 r sin( 4k) | 0 sin( kz) k+6 
} > ’ "k 9 >\* 
z k(k-+-1)(k?—4) 7 k(k+-1)(k+2)(k?— 16) 


lhis result is a special case of that given, for a supported semicircular plate, 
by Timoshenko (Theory of Plates and Shells, p. 272) except that the first 
three terms in (8) are there expressed as a Fourier series of the same type 
isin the remaining part of (8). The functions p, are now to be chosen; they 


n must form a complete set of harmonic functions in the polar coordinates 
‘ . 
092.11 mn 
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(r,@) and the functions v, are then to be found as the solutions of the equa. 


3. In 
if (Pmt 


tions Av, = p,,v, = 0onC. The values of v,, were originally obtained by 
indirect and complicated methods and the corresponding p, were they 
deduced, but the former were found to be linear functions of comparatively | _ 
simple expressions which were therefore taken as the appropriate v,,, and It m, vi 
then the p, are the elementary harmonic functions p, = 4(n-+-1)r” cos né und mAs 
where n is zero or a positive integer.+ 

When is odd, the corresponding v; are given by 


Uk (r*§—r*+2)cos k0, (9) | When 2 


which vanish when r 1 and when @ $77, 1.e. on the boundary C of | and wh 


the plate. 


When nv is even, the v, are more complicated; for n = 0, 2, 4, 6 they are 

Vy = 1°(1+- cos 20)— ¥ B,.7r* cos k Owing 
Vo r4(cos 26-+-cos 40) pt C,.r* cos ké (10 then tl 
he r®(cos 46-+-cos 60)— ¥ D,7r* cos ké | approx 


lea * I 








Ve r°(cos 66-+- cos 80) S Erk eos ké | | 
6 he Kk 
where the infinite series are those required to make v, zero when + = 1, 
their values being clearly zero for @ }7. The constants in these series 


are 


. 4 »)2 
B sin(4k7z) a 
: 7 ”  k(k®?— 22) 


4 
C; sin(4kz) 
h os (9 (k? 22)(k2 4?) 


D, sin(tbey —(" —* _ 
, 7 " (k?—4?)(k?—6?) ‘ 
> ge le su 
E,. * aintdibe) = all - : 
; 7 ‘ (k?— 62)(k2— 8?) 
It is possible that linear combinations of the p,, could have been chosen 
as new functions (P,,, say) and that, with these P, and the corresponding and st 
V,, the approximation (3) to the exact solution might give a more accurate | the nu 
result. This could only be tested by numerical calculation and it will be ‘Un 
seen that the convergence is fairly satisfactory with the p, as given above. | #Y&? 
It will be noticed that the approximate solution (3) is composed of terms 
of the form r” cos n@ and r”** cos n@ which are elementary solutions of the 
equation AAd = 0. The effect of Weinstein’s method is to enable these 
terms, or combinations of them, to be arranged in an order of importance, 
though this cannot be done without introducing infinite series as these 
occur in the v, when » is even, and also in W. 
+ In the application to the semicircular sector k will always denote an odd 
positive integer. In all cases m, n are positive integers, including zero. 
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equa | 


3. Integrating over half the area, because of the symmetry, the value 
ed b Pm) is 16(m-+-1)7,/8, except when m 0, in which case 
the ' 
; Pos Po) l67/4. 
tivel : F 
lif m.n are both even or both odd and m + n. (pp... p 0. If m is even 
an Pm Ph 
nd nv 1s odd. 
‘OS nt . 
16(m n ) n 1 . 1 
Pps P . , COS(5M7r)sIn(5N7). 
m n Z n=—m-~ si 
») . l 
; ~sin(3?77) 
Q When nis odd (] v — 
. n+-2)(n 4) 
Vy Ce nd when 7 is even 
1. v, 0-14877 84 (i. o,) 0-03902 14 
"y at 1,4 0-01414 7] (1, v%) 0-O00687 145 
Owing to the factor 16 in (p,,. p,,) it is convenient to put A’, 16A,, and 
the values of A)...., Ag, as the solutions of (7) for the successive 
1} 
roximations given by m By Bs ks 6, are 
series <n = 
11 
; ae 
(1] 141 H 11 
847 |? I ) 
successive approximations to the solution w are then 
W, Wy W +A, %, Wy W+Ayv)+A,%, 
hosel “6 YW 1g % A, Vi + A, YG: 
adin | should make the values of V on the boundary approach zero as 
urate the number of the constants A ncreases. 
ill be On the boundary 1 the values of cw/coN for 6 O. is... 75° are 
hove, | given in the following table. 
| 4 7 
erms TABLE | ae 
»f the G @ 15 6 6 @ 75 
these 19 I 07 024 
ance 3 Xe tees 


these 
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‘imi j . | Critical 
Similarly the values of éw/én on the boundary 6 = 47, calculated fo} \* : 
r = 0, 0°125,..., 0-875, are shown in the next table. —— 
« small 
ras.e I] the sec 
1 : \ 5 ’ : | making 
ev 
aN 24 757 $25 747 2 17 039 13 095 
4658 | —-004 34 3 413 | —*001 966 | —-000 159 | +--001 75 3 342) +7003 The re! 
Ow, express 
\ 5 543 § 22 $313 " 2 5S ° I O2 *000 535 2 456 . ae | 
same V 
c 
\ I 37 I 27 } 124 765 I 
ON y 
: { 18 1g 00 385 
= 8 | —-000 094 0 054] -+° : : 079 | +*000 134 12 ; where 
c \ 
‘ case W 
; M ae ‘ al ae ers Pee cl ics ee 
ON as indi 
. c 2] ‘ > | 7 22 “Oc 7 Teter 2 : : H,isn 
- — is ("+ 
_ ’ 2 ’ , ; ; | the bo 
lhese tables show that, taking m 5 or 6, the approximation given by If 
(3) is sufficient to give results to an accuracy, in general, of about | in 1,000 
if the accuracy is measured by the ratio of cw,/eN toeW/eN. In the tables, 
. ° > . — > we f 
this ratio is somewhat smaller for @ is", = 0-715 and » 0-875.) (to tal 


indicating that the approximation becomes worse as @ approaches 90° and! functi 


r approaches 1. This is to be expected as it is known that at a sharp bend} cases) 


on the boundary w cannot be represented by an analytic function. | specia 
The maximum deflexion, which must lie on the radius 6 = 0, was} from ; 
obtained by calculating w for 7 = 0-1, 0-2,..., 0-9 and deriving from then 7 
the maximum value. For m = 5 this is 0-002 022 at r = 0-48594, and for| ae 
m — 6 it is 0-002 023 at r = 0-48587. — 
since 
4. For the general sector of angle 2x and radius a the expression for } 
corresponding to (8) is | 
9 } 
l i 1 r+ ecos 26 1 + cos 40 i. 3 F. (") "i ko Re a,() cos k6. | Hence 
64 48 cos2a 192 cos 4a i " Vel int } 
(12 
where k (2n+1)a/2a, n > eee and 
F a‘ sin ko asinka | l 4 o } This 
 Qak(k+-1)(k2—4) xk ae 2 k+i1'k a} 5, rem 
G, a‘ sin ka (k- 6) | 
2a k(k+1)(k+2)(k2?—16) 
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Critical cases arise when cos 2« = 0 or cos4a = 0. Taking for example 
.=4}n, then k 2 6. 10... so that K>ow as k->2. If k 2+-e, 
« small, then since 2ka 7 as n 0 for F,, it will be found, considering 
the second term in (11) and the first term in > F,(r/a)*+*? cosk@ and 
naking « > 0, that these terms together lead to the expression 
1374 cos 20 r* log(r/a)cos 26 r4@sin 26] 
L447 127 


lhe remaining values of F, and all the values of G, are given by the above 
expressions for these coefficients. The values of v,, are determined in the 


same Way as before from those of Pp 4(m-—-1)r™ cos mé. The result is 


m 
COS Ma ‘ r\k 
U,, r™+2) cos mé cos(m 2 1 S M,{ cos ké, 
cos(m-- 2)a oo a 

(13) 
where the H,. are to be determined so that v,, = 0 when r = a. A critical 
case Will arise when cos(m-+-2)« = 0 and has to be treated in the same way 
sindicated above. For certain values of m only one additional term in 
H,isneeded; for example, when a lor, the value of v,, for m 2, 6, 10,..., 
2__q?r™)cos m@, which is zero when r = a and when @ = }z, i.e. on 


the boundary ( 


If the loading is not uniform but can be expressed in the form 


| S r™ cos nO 


to take the case of symmetrical loading) the determination of W and the 
functions v, can be dealt with on the lines exemplified above. Critical 
ses will arise depending on the values of m and x and have to be treated 
specially, either directly in any particular case or by a limiting process 


ma general result. 


5. Returning to the case of the uniformly loaded semicircular plate and 


putting w (r4/64)+-¢4, then AAd 0 and ¢ 4/64 on C. Again, 
ce cw/oN 0 on C, it follows that ¢d/ér jr on @ t lz and 
| Cd 
Vv on a a 
r oo 
Hence ¢ is a function satisfying AAd = 0 in S, d 64 on C, 
Cah . ch 
\ jr? on 0 57, and \ YQ on 7 i. 
Ol C4 


is function can therefore be identified with the stress function in plane 


tress or plane strain, in which on the boundary 
— ~ l od l od 
iT", ré = 0 
r 00 ~=rerceé 


LA f 
l c-b l od 
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and oe 
c 2h or 


66 
or? 16 
so that the plate is acted on by the uniform normal thrust ; on the sen 
circle 7 1, and by the equilibrating normal thrust #7? on the straight 
edge of the plate. 

It is evident that the sulution of the problem of the deflexion of a clampe 
plate loaded in any given manner can be used to derive the stress functic 
for a plate acted on by a derived system of stress on the boundary. TI 
converse problem of finding the stress function ¢ for a given stress syste 
on the boundary, i.e. for given values of ¢ and ¢¢/eN on C, requires furth 
consideration; it does not appear to be in all cases directly and easil 
transformable into the problem of the deflexion of a clamped plate, load 
in a manner which is to be deduced from a knowledge of the given valu 


of ¢ and é¢/éN on the boundary. 
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ISOLATED FORCE PROBLEMS IN TWO-DIMENSIONAL 
ELASTICITY (1)t 
By M. ROTHMAN (Royal Technical College, Glasgow) 


Received 6 S« ptember 1949] 


SUMMARY 


Stevenson (1) has shown that tl tresses and displacements in two-dimensional 
sticity can be expressed in terms of two complex potentials Q(z) and w(z). In this 
er Stevenson’s method is used t ive the problems of one or more isolated forces 
the boundaries of holes of various shapes in an infinite plate, agreement being 
hed with spec¢ if Cases previously solved by different methods. 


1. Introduction 
Ir D is the complex displacement u+iv, © and ® the generalized plane- 
stress combinations x2+iyy and xx—yy+2ivy, then in the absence of 


body forces we may write 


8uD KQ)(z) —202'(Z)—@'(2Z), (¢. 5) 
20 = 0'(z)4+- Q"(2), (1.2) 
nd 20 = 20" (2)+0"(2), (1.3) 
where » is the rigidity, one of the Lamé constants (A, ~), and x = 3—4o, 
shbeing the modified Poisson’s ratio. Dashes denote derivatives. 
2. The cyclic functions of the complex potentials 


For single-valued displacements, i.e. solutions involving no dislocations, 
juation (1.1) gives 


CyfrQ(z)— 22! (Z)— a (2)! 0), (2.1) 


where Cyd denotes the change in ¢, a cyclic function, on passing once 
ind a closed circuit so as to keep the area enclosed on the left. 
Differentiating (2.1 


and using the fact that the stresses being single- 


ued, © must also be single-valued, and therefore that 


Cyf Q(z) +0'(zZ)} 0, (2.2) 
juation (2.1) gives CyQ'(z) = 0 (2.3) 
nd therefore C'y@ (2) KC'yQ(z). (2.4) 


Ifnow X+7Y and WM are the complex force and couple resultants at the 


rigin equivalent to the stresses round a closed circuit, then (1) 


XA+4Y LiCyfQ(z)+2Q'(Z)+-a'(2)}, (2.5) 
d M iM * i Yi wz) zw'(z)}, (2.6) 
Forms part of a thesis approved for the degree of Ph.D. by the University of London. 


(Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 3 (1950)] 
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so that, if F = X+ iY, then from (2.4) and (2.5) 

F hie + L)CyQ(z) or CyQ(z) 41 F'/(«-+-1), (29 
or changing the sign of 7 and using (2.4), Cyw’(z) 4ix F/(«-+-1) whence 
from (2.6) - 

Cyw(z) 4ix Fz/(e+1)+4(M+iM*), (2.8 
ARn 
/ 





3. Form of the complex potentials for the infinite plate 


We assume potentials of the form 





Q(z) = A,logz+Az+ >} A,2-" (3.1) 
1 
and w(z) = Czlogz+ B,logz+ Bz?+ > B,z-*, (3.2)] 
I 


where A is to be real since Q r+iBz, w = y+6z (f real) give rigid body'| 
displacements. 

The terms in A and B give finite stresses at infinity (but will inevitably 
give rise to infinite displacements there) and there can be no other positive 
powers of z, the negative powers having less effect as we recede farther 
from the origin. 


Assuming the force and couple resultants for the circuit at infinity to | 


be F and M, equations (2.7) and (2.8) give 271A, = 47 F/(«+-1) and there- 
fore A, 2F /a(k+1) and 27iCz+27iB, 4ix F2(x+1)-1+4(M+iM*) 
and therefore C 2«F /z(xk+1), By 21(M+iM*)/z, so that 
Se y* f 
ie Es on Fa (3.3) 
a7(«k-+-1) _ 
2« Fz log z ; B, 
and wv(2) A °8* _ (2iG/z)logz+ B2+ ¥ B,2~, (3.4) 
am(x+1) sig 


where J is the real part of G. 


= 


In particular for an unlimited plate with a force and couple nucleus 
and G at the origin, we have 


ee } 
i = (3.5 
m(«k+-1) 
2«Fzlogz 2iG a 
and w(z) 7s = log z. (3.6) } 
a(k- 7 
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4, Isolated forces on any boundary 


Taking complex potentials 


()(z) A log z, w(z) Azlogz, (4.1) 
then rr+-00 (-) ©) 4(A/z-+ A/2), 
and rr —00+-2ir6 — ® Me-20 (Z/z)® 4(A /Z+-A/z), 


which is real, so that r@ 0, and 7r—00 rr+-00, whence 06 = 0, and 
}(A/z+A/z). The potentials given by (4.1) are therefore seen to give 


a simple radial state of stress. 








It is clear from the previous results that if the origin be taken inside the 
material, the solution will not give single-valued displacements and so will 
be non-physical. If, however, the origin be taken outside the material, 
then no circuit within the material being reducible to the origin, the 
objection does not arise, and the origin can even be taken on the boundary 
of the material. 

If the tangent at the point O on the boundary of a material be given 
by 6 vand @ = w—a, then considering the equilibrium of the material 
within a radius 7 of O, the forces on the boundary APB clearly have a 
resultant passing through O (= X-+cY). The forces on the boundary BOA 
consist of a finite distribution of boundary tractions over the are AOB and 
a possible isolated force F (— Re'?) at O. Hence from the equilibrium of 
the material in the circuit APBOA we find 


7 x 
F+force resultant F’ of boundary stresses over AOB+ | rre“r dé = 0, 


x 


F+F’'+ | 3(A/z+A/2)z dé = 0, 


01 F+F'+}3 | (4+Ac%4) do = 0. 
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This is independent of 7, and therefore if 7 > 0, F’ > 0, and F+47A = ( 


whence A oF |r. (4.3 
Hence the complex potentials for such a force F on the boundary at th 
origin O must (from 4.1) contain the terms 


2F 2F 
Q)(z) log z, w/(z) zlogz. (4.3 
5. Change of origin 
The displacement D is given either by 
8uD KQ(z)—20'(Z)—a' (2), 
or by 8uD KQ,(2,)—z, (5 (2) @,(2,); 
where we have the two origins (1) z 0, and (2) 2 0 given by z ¢ 
i.e. 2, = z—c. Therefore 
KQ(z)—20'(Z)—o'(Z) = KO, (z,) —22,(Z,) — {4 (Z,) —eQ)(Z,)}, 
so that Q(z) = 2,(2,), w(z) w,(2,)—E€Q,(2,), (5.1 


and these are the formulae for the stresses for a change of origin. 


Therefore, for an isolated complex force F at a point z = z,, putting 
ar = & : . 
2X0 logl, — w,(¢) Clog ¢, 
7 7 
oF. 
so that Q)(z) log , (5.2) 
and oF oF 
w(z) —— (z—z,)log ¢ : — ioe ‘| 
2F a : 
ne zlog f+ —(2 »— Fz,)log ¢. (5.3) 
Now FZ, = (X+iY)(xv,—iy,) = Xx,+Yy,+i(Ya,—Xy,), 
so that FZ,— Fz, = 2iG, where G is the moment of the force F at - 
about the origin, and therefore 
2F P iG ‘ c 
w(z) zlog f+ "log ¢. (5.4) 


6. Circular hole in an infinite plate under an isolated force on the 

boundary 

For a complex force F = Re’? at the point A, (z = a,, a, = ae’) on 
the boundary of a circular hole in an infinite plate, equations (5.2) and 
(5.3) give the complex potentials 

2F 


7 


7 


Q(z) log(z ay) (6.1) 
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2F 2 
\ and w(z) — z log(z—a,) + —(Fa,— Fa,)log(z—a,), (6.2) 
(4.2 ri 
and these give the requisite force at A,. 
it th 
(4. 
(5.1 The stress resultants round the circle at infinity must balance the com 
» . 
plex force F at A,, hence the logarithmic part of the complex potentials 
ist contain terms obtained from (3.5) and (3.6) by changing the sign of 
tting Fand G | Li Fa, Fa,)|, i.e. 
O 2F log z (6.3) 
ScZlzZ die 
m(«-+-1) 
~ 9\ 2x kz LOg Zz l i 5 
(9--) (and w/(z) (Fa, Fa,)log z. (6.4) 
calles 1) i 
| Hence the terms Q(z) A log z and w(z) Bz log z+ C log z must be added 
the solution (6.1) and (6.2) so as to give the form (6.3) and (6.4). 
(x9 If 2 z—a,, then €/z |—a,/z, and so (6.1) and (6.2) can be written 
.») i 
oF 2F a 
()(z) log z lou(1 | 
2F es 2F ; a 
w(z) z log z (Fa,— Fa,)logz- =log( | 
( 4 » a.) 
(Fa, Fa,)log{ _ ‘). 
the | 


e terms in C/2 |—a, z vanishing at infinity. Hence, comparing the 


int solution with (6.3) and (6.4 


2F 2I 
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2F K 2F : Bes 5 
or A ; B ; ( ——(Fa,— Fa,), 
m(x-+1) a7(x+-1) T 
and so the joint solution which has the correct form at infinity and gives 
the isolated force at the boundary is 
2F, .  2Fxlogz 
9) Piogt \ Fx log 4.Q,, 
i m(k+1) — 
2F ee ; » 2lzioee 1... ; 
w zlog €4+ —(Fa,— Fa,)log ¢ ©— _ — (Fa, — Fa,)log z+-wy, 
7 T m(xk+1) @7 
where Q, and w, are expansible in negative powers of z, and are chosen to 
satisfy the boundary condition 
Q(Z)+20'(z)+w'(z) =k 
round the boundary 22 = a? = a,d,. Now ¢ = z—a,, therefore 
c Z—Ay a,(a/z—1) a, ¢/2 
so that the boundary condition becomes 
2F ( a, *) 2Fxlog(a*/z) . .( 2F 2FK | 
og - +— = st 
TT 2 a(x+1) | 7 a(k-+1)z} 
2F . Zz 2 — = 2 F (log z 1] l = 3 
(log ¢ | ] |. —_ (Fa, — Fa,) (Fa,— Fa,) t 
T C 7 a(k+1) TZ 
" Q),(Z) +20,(z)-+-w,(z) = | 
and as ¢ a,¢/z, and 
z/C a?/(Cz) (a*/a,)(1/z—1/€) a,(1/¢—1/z), 
the above condition can be written 
— 
| ; - ; K l 4F 
2F log €((1—1)+24 log 2/1 3 | 
TT | k+l «+I1) ° («+I1)loga 
_ . 2F ott 8 2 Fa? 
2 F log(—4d,) 2Fa,(- { - 
K+] a Zz («+ 1)z* 
| 2F(f-+a,)/C+-2(Fa,—Fa,)/C—(Fa,— Fa,) | | 
t (0, (2) +206 (z)-+-w,(z) = &. 
The constant terms log¢, logz, and the terms in 1/¢ cancel out and the 
constant terms together with the constant k can be replaced by a new | 
, ae { 
constant k’, giving 
Fa,+ Fa 2 Fra? ) 
= SC)’ (- ar 1 1 = p! 
Q),(Z)+-2Q,(z)-+-,(z) 4 = £'; 
TZ a(K-+ 1)z* 
this can be satisfied by taking Q,(z) 0 and 
| — = 2 Fra* 
W(2) (Fa,+ Fa,)log z+ : 
7 a(«-+1)z 
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s this gives finite stresses and displacements at infinity, and, since the 
coefficient of logz is real, Cyw(z) is a pure imaginary, and so the moment 
FIVES condition is not affected. 
Combining our solutions, we have 
2F 2 Fx log z _ 
(2)(z) log(z a,)- eg (6.5) 
7 m(«-+-1) 
a ») 
iii w(z) ~ 2 log(z—a,) — (Fa, Fa,)log(z—a,) 
2Fzlogz See 2 Fxa* 
on te = Fa, log z (6.6) 
mxk+1) 7 m(xk-+-1)z 
For a number of forces F. at z = a, round the hole, which is otherwise 
free from applied force, these complex potentials become 
2 XN p 2x log z XN p die 
()(z) F log z,+ F (6.7) 
TT kof am(k+1) a 
w(z) = = Flogz,t+= > (F,a,—F.a,)logz, 
2 log z N i 2 log z N rz 2Ka" —"» sa 
P Fa, 4 F., (6.8) 
) a(k L) Sone 7 yaw a(x«+1)z aw 
the terms containing S F. disappearing when the forces are in equilibrium. 
A 
Equal and opposite forces at the end of a chord 
If there are two forces F, = iRatz = a,,and F, Ratz = a, = a,, 
ne complex potentials reduce to 
| 
| ) ny 
O 2th log(= ay 
7 rd ay 
— es A 21 Rz z—a 
nd as N Fa iRa, l Ra, Q, WwW log| '. 
a. 7 -—Qy] 
Now 2(nn-+-8s) = Q'(z)4+'(2), and as mm = 0 round the boundary 
2 g” a,a,, 
tl the peripheral stress ss is given by 2ss Q)'(z) + Q'(a?/z), 1.e. 
1 
new SS i R/7)| (4,—a,)/ (z?—a?—a, z—d, z)]+ conjugate, 
nd if a, ae’™, z ae’” round the boundary, this becomes 
! , s 
2Rsin «sin 6 
7a(cos 8@—cos x) 
It will be noticed that the peripheral stress is infinite at the point of 
| application of the isolated forces, except in the particular case a Lar, 


when ss 2 R/ (za) constant round the hole. 
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7. Elliptic hole in an infinite plate under an isolated force on the | 497° 
boundary potent 
The space outside an ellipse in the z-plane is conformally mapped on the 
outside of the unit circle in a complex o-plane by the transformation 
z c(o+Ac—). (7.1) 
i.e 
Z since t 
1 which 
| so agT' 
— 14D 
where 
Fic. 4. The e 
We write ¢ a—o, and suppose that a complex force F Re’? acts 
at the point A, (z = 2,, o = o,, or € = 0) on the boundary of the hole. | and 


Considering first the complex potentials which give the isolated for¢ 
A, in the z-plane, we take 


Q, A log ¢, Wy Bz log €+-C log ¢, 


where the constants A, B, C will be chosen to give the required force at Aj. 


This force together with the stresses round a small shaded region 


equilibrium, so that if the stress resultants along this are are X +7Y an 
then PLX ti) 0 
and 40(F'z,— F2Z,)+M = 0 
where X+iY fH{Q(z)+-20'(Z)+@'(2Z)}, 
and M | real part of {w(z)—zw'(z) 220 (z)}"1, 
and as we are only going round a point semicircle, 
[Qy; £CyQ, giving F lor( B—A). 
Also 40(F'z,—F2Z,)+1} real part of (277C) 0 
so that if C Ci+iC,, then i, (2/7)(FZ,— F2,). 


€ al are th 


ofca 
(7.2) 

can b 
is in 
iM, g 

nas ti 

2(o)) 

wher 

the r 
q solvi 
ee! 

A 

(7.4) 
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the | As 7 > ©, logz + logo, and therefore at infinity the form of the complex 
tentials is 


| 2F log z : 2F loge ae 
1 the ()(z) —, 20. {2) ——. (7.5) 
m(«-+ 1) m(«-+-1) 
2x kz log rs | —_ ; 
e w(z) (FZ, — Fz, )log z, 
(7.1 a(xk+ 1) T : 
2x Fzloga | — ; 7 @ 
w/(z) (Fz, — Fz, )log a, (1.6) 
(kK 1) 7 
| since these will give, round the infinite circle, force and couple resultants 
| which balance those of F at A,. Therefore, taking* Q, D logo, and 
Ezlogo+G logo, the combinations at infinity will be of the form 
Q) A+ Dy)logc. ( (B+ E)z log o+-(C+ G)log o 
wreeing with (7.5) and (7.6) if 
IF ») iy ] 
; y ’ =K J i ° _ ; — 
{.D B+-E , XC.+4,) (Fz,—Fz,), (7.7) 
TK 1) m(xk+1) ss 7 7 
ere ( ( 1 i¢ 2 and G G, t iG. 
[he complex potentials 
‘ail Q) A log €(+-D log o+Q,(o), 
hole | Ww Bz log C C log C Ez logo + G log o W,(o) 
ce al re therefore taken, where Q, and w, must be developed in negative powers 
yand any rigid body displacements, i.e. terms of the type 
i (2 v+2B2 3 real), Ww yA dz, 
t A n be neglected The boundary condition 
is 10 Qo 2(o 1) dQ lz dw dz constant k, 
dM, “(a §{Q(o-) +k} +-2(0)Q'(o)+w'(o) = 0 (7.8) 
is to be satisfied Substituting for 2 and w this reduces to 
Q),(0-!)+ A log €—(A+ D)log o+k’'+ Blog (+ EF logo! 
2(o0-1){O0(o)+ A /€+- D/s'+-w,(o)+2(B/€+ E/o)+C/6+G/o = 0, 
here i’ kA log( G1) constant, and choosing 
B+ A 0 and #H—(A+D)=0 (7.9) 
the rigid body displacements in this equation can be eliminated; thence 
: ving equations (7.3), (7.7), and (7.9) we find 
(io 


A - pu. D << gu. aart (7.10) 


T 7(xk+1) 7(K-+ 


(7.4 * Here, D must not be confused with the same symbol in section 1. 
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The boundary condition then reduces to 


w(o)+ C/f€+ G/o+e(o+Ao-!)( B/E+ E/o)+{Q)(o-1)+k'tz’(a)4 


{Of (a)+ A/€+-D/o'e(a—+-Ac) 0), 


and writing o/¢ k’ so that Qi(c) = 0 and putting B A, this 
further reduces to 


wy(a)-+cfA(A + D)+(E—A)'+¢(D+AE)o-?4 


+ (C+ Az, —Az,)/€+{Go,—c(A—AA)}/(c0,) = 0. 


The terms in ¢-! and o~ in this equation can be made to disappear by 
choosing . - 

‘ Go, c(A—AA) (7.11) 
and C Az,—Az, 2(FZz,—F2z,)/7 (7.12) 


so that C is a pure imaginary (= 7C,) and this is the value given in (7.4). 








If, ho 


real, a 





Also from equation (7.7) iG, (F2,— F2z,)/7 which can be shown to be | 


consistent with (7.11). Hence, 
G 2Ac(Fo,+ F6,)/7—2F2Z,/n (7.13) 
and the boundary condition is satisfied by taking 
W,(o) c{\(A+ D)+(#H—A)'o+e(D+AE)o— 
to which may be added the rigid body displacement 
W,(o) c{A(A+D)+(H—A)(o+Ao-), 
so that . : , 
W,)(o) c§{D+AE+2(H—A)+A2(A- D)\oa 1 


= { F(«—A?)+ FX(«—1)!. (7.14) 


on(«k-+- 1 


Hence the final solution for the problem of a single isolated force on the 


boundary of an elliptic hole in an infinite plate is 


Equ 
If t 


nd z. 


80 tha 


and 


Then 


hound 


where 
and Pp 
whose 


form 


Se Al > FY.. oc - - 
Q oF og tie ee C o—0d}), (7.15 
iT a7(x+1) i 
2F poe tess ; , 2Fzlog 
w log 2 (Fz, k log Z al 
T 7 m(«x+1) 
2Ac . = 2Fz 2c 7 9 ; ) 
“ (Fo, + Fo,)—- tllog } ; { F(«—A?)+ FXA(«—1);. 
\ a a | o7m(«k-+1) 
(7.16) | 


5092.1 
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| | numbe 4 of Torce s in equilibrium 


| For a number of such complex forces F. at z z, (corresponding to 
c,), these potentials become (writing ¢. = o—o,) 
% 2 = 2x loge <- 
Q S Flog f,4 + 2. (7.17) 
this TT hod a(«+1) Ly 


22S Flog t, += 


M 


ok ; : 2zloga ~* = 
(Fz F ey log be - S F’ T 
T(K T 1) 4 


») 


A ; S Rs 
’ | ~~ (Fo Fé.) > F,2,|log a 





- r 7 ar d 
ar by 2c 2K \ = | — 
(x—A?) F. 1_A(«—1) 2 F}. (7.18) 
a om(k 1)| 7 
tacks 
| If, however, the forces are in equilibrium, then > F. = 0 and > FZ, is 
(7.12 | il, and the complex potentials become 
(7.4). | ») ‘ 
to b 7 i I log ¢,, (7.19) 
Ww A S F log e. = S (F z. F’ z,)log a T 
(7.13 ills ee 
—log ofde > (F,6,+F,0,)— > F,Z,}. (7.20) 
Equal and opposite force s at the ends of the minor axis 
If there are two forces, F, = iY and F, iY, at the points z, = ib 
ib, then 
> (Fz F Zz) 0 and > (F.o, , F’G,) 4Y, 
; 2iY o—t 
that Q) log |, 
7T Co- t] 
(7.14 
PY 2 o—t 4aY 
Ww log| ] logo. 
1T Co-+t 7 
yn the 
hen as 2(nn-+ss) Q)'(a)/2'(o) + Q'(6)/2'(&), and nn is zero round the 


indary, we get on simplification, 


(7.15 2Y(1—A) 


(ss Mb 


oundary 7c(1+-A?— 2A cos 2n)’ 
ere ¢ = cosyn+isiny. Also if a and b are the semi-axes of the ellipse 
id p is the perpendicular from the origin on to the tangent at the point 
hose eccentric angle is y the a e result can be written in the alternative 
l rm . a 
. 2Yb 2y p? 
(7.16 < (a? sin*y +-b? cos?n) mab - 


092.11 


U 
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This is a maximum at the end of the major axis where it has the valy 
2Y/7b, and putting a = b or A= 0 in the above results, ss = 2Y /7 
agreeing with the previous result for the circular hole. This result is th 
same as that derived by Green (2) and can be written in the forn 
(ss.7a/Y) = 2a/b which has the values 2, 4°33. and 8-16 for b/a 1, 0-462 
and 0-245 respectively, agreeing with results otherwise calculated }y 


Symonds (3). See Fig. 6, where h = b/a. 


Equal and opposite forces at the ends of a diameter 
In this case there are two forces F, at o,, and —F, at oy (o, o 


where F, is complex. Substituting in the above results we find 


2F o—o 
Q — log( zi. 
7 Oo- O71 


‘ “ . i 
siine = 2 R{Asin(B+y) sin(B )| . 
mc(1+-A*— 2A cos 2)sin( — a) 

where F = Re'® and o, = e**. When «= 8 = 4a this reduces to the 


previous result for forces at the ends of the minor axis. See Fig. 7, drawn 


for B 45°. 


Equal and opposite forces at the ends of a chord (i) parallel to the minoi 
axis, (ii) parallel to the major axis 


(i) The two forces are iY at o = o, and —iY at o = o, = 6G, giving 
amr 
27) o—o 
QO=- log 1, 
7 oO O71 
and writing o = e'?, o, = e'*, then 
™ 2Y|A sin 27 sin 7 sin a+sin «cos 7)| 
SS z 
7c(Cos 7 — COs a)(1+-A?— 2A cos 2n) 
(ii) The two forces are Y at o = o, and —Y ato = &G, giving 
Se 2Y|cos 7 sin a(1—A cos 27)—Asin 7 cos y sin 27] 
Ss r 


m™c(COs 7 — Cos «)(1+-A?— 2A cos 2n) 
and when A 0, both these results reduce to that for the circular hole, i.e. 
2Y sin «cos n 


7C(COS 7 — COS x) 


8. Curvilinear polygonal hole in an infinite plate under isolated 
forces on the boundary 
The conformal transformation, z = co(1-+-Ao~”) transforms the outside 
of a curvilinear polygonal boundary (Fig. 5) into the unit circle (4), where 
A < 1/(n—1). 


We assume a complex force F = Re’? to act on the boundary of the hol 
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> value} at the point A, (z = 2,), corresponding to the point o = o, on the unit 
2Y /7a,} circle in the o-plane, where o—o, = ¢ and o, = e'*. 
b is th Also o e’ round the boundary of the hole, and we call » the eccentric 
e forn}angle of a point on this boundary by analogy with the eccentric angle of 
, 0-462 | the ellipse (which is the special case n 2). 


— ot) 





— 








to th 
drawn| Fic. 5. 
Proceeding as in the case of the elliptic hole, complex potentials 
— Q, A log ¢, Wy Bzlog €-+C log f (8.1) 
re first taken and values for the constants A, B, C are chosen such that 
Ing . . , . 7 ° 
these potentials give the required complex force F at z = 2,, i.e. 
5) = 
- 45 _ 2 ome ; » 
I la7(B—A), iC, = —(F2z,—Fz,). (8.2) 
7 
Secondly potentials 
Q, D log c, We Ezlogo+Gloga (8.3) 
we taken so that the combinations Q,+Q, and w,+w, give the correct 
form at infinity and so provide the necessary balancing force and couple 
lor the equilibrium of the plate, i.e. 
. QO 2F log . 
le, i T(K- 1) 
2k Fz log ra 21 — 
W - (M t iM*)log z, 
a(x-+- 1) T 
plated; Where M is the real part of i Fz, }i( FZ,— F2,) and as in equation (7.7) 
oF 7 2x Fk ° ’ y 1 — = 
1+D , B+E ' i(C,+ 4G) (Fz,—Fz,), 
utsid a(xk+1) am(«k+1) T 
wher - Bee a ' 
nd iG, — (F2Z,—F2,). (8.4) 
1e holt 


Potentials Q Q 


-Q, and w = w»+w,-+w, are now taken, where 

















292 M. ROTHMAN ISOLAT 


‘ , : . | 
Q, and wy are expansible in negative powers of o only, so that the cop! and the 
ditions at infinity are unaltered, and Q,, w, are chosen to satisfy t} 
boundary condition round the hole, i.e. 


2'(a){Q(a-!) +k} +-2(a)Q'(a)+-w'(o) = 0 round o6¢ = 1. = 
Now, C a(1 = so that ¢ i . lor 
oO O71 C 
and the boundary condition reduces to eile 
2 (){Q29(0-1) +k’ + A log (—(A + D)log 0} +-2(0 1)\fOi(o)+A/C+D/oi4 agreein 
z'(a){ Blog (+ E log o}+4(o)+C/€4+ G/o+2(c)(B/6+ E/o) =| | : 
(a) Tec 


Choosing B+ A 0 and E—(A+D) 0, the logarithmic terms disappeai 
and these equations together with (8.4) can be solved for A, B, D, and BE) J\9): 
whence the boundary condition further reduces to 


2"(a){Qo(a-1) +-k’} +-2(a-1)Q)(0) + a9(0)+f(0) 0, 
where 
C G ’ 
ft(c) — | ¢H(1+-Ao-")—cEk(o-?-+Ao”-?)4 Now w 


Co 


CT 


can be 


; l Ao” 1 — lo A m 
ee (an Meee CA 2 . (8.5)) 
\FS S ; a w (9) 


The terms in 1/¢ and 1/o, which would only arise from the logarithmic whene 
terms already allowed for, have to be removed and f(c) expanded ir 


}(o) 
powers of o. 
Now 
l iyi 4 o 1% 
Co a,\¢ oa C ‘4 
gt-l eo Now t 
: go" 2 1 go” 3 x o” 1 o” 2 1 
4 C 
and 5 | 
8.0) D 
l tS Sy) 3 l 7) l 
Co"! o*-8 g,\C = g,0°*o,\¢ o a, o*-! (af 


: fl l i 7 bE ‘ia _ 
ot rf: } ate ot *e Panos af * «0° 
so that f(c) can be written 
C WC + Ac(o; 1+-Aot-) Ac(o,+-Ao;" 1)] |g UG -c(Ao; 1_ Ao; n+1)| 
Le(# A) +-cAAo}-*-+-co | Ex AAoj” aa 1 eg”—2\(A — Ex) 
| cAAlo, o”—3 4 oF ot 44.4 ot 3] 4 
| cAA[o; nt+3q-34....+o; 2o-"+2+-a7 10 sie! 


also C and G are chosen such that 


C = (400, = Az,—AZ, = 2(F2,—F2,)/z, 
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1€ CO} ind therefore comparing with (8.2) 
sfv = — v 
( Ws 2( Fz, Fz,)/7 (8.6) 
nd 
(0 cl 1c; I AAo; n+l) Ac(a, " 14 Ao”—)4-¢ {(c; 1 Ac} 1) 
: G = G,tiG, = 2dAc( Fo} 1+ Fet-!)/7—2F2Z,/2 
therefore 1G, (Fz, Fz,) 7 (8.7) 


voreeing With (8.4). 
\ll the constants A. B, C, D, BE, and G are therefore determined and 


t reduces to 


bppea 
and f(a) c| E—A-+-AAo} 7 Ac(A — Ex)o"-?—c( Ex—AAoj"**)o 2 
CAA ot 3g ot 4g? +... oF o” 44 g,0" 7 | 
cAA| op" *3a-4+- 0; to -4+-...0] 20 "t*+ 0, 'o ati * (8.8) 
Now w .z, Where « is a constant, gives a rigid body displacement which 
in be ignored, so that we can add the solution 
(8.5 7 n_9 ’ ‘ 
' o (E—A-+AAo}-?)z'(c) c(E—A+AAot-*){1—(n—1)Ao-"} 
“ithm vhence (3.3) becomes 
ded I y 1-2 | 3 2 ~n—4 n—4 2 n—3 
Ac(A — Ex)o"-*-+-cAA] 0, 0"? + a7 6 +... 07 “0°+ oF 3q] 
co~?( Ex—AAGt**)- crAA| oj nay 34g tig 
a] 2a-" 2+ oF 1a" +1] 4+-c(n—1)Ao-"( E— A+-AAo}t-?). 
Now taking a — 
£2, > A.o J Wy Cc > B,.o r cB 19; (8.9) 
i 1 
8.5) becomes 
Q(o-1)+kh’!+-Z(o6 1) (o) + w(c) 
| , 


A(n 
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so that from (8.8) and (8.10) the boundary condition can be satisfied } 
taking NA—Ex)+4, (0. ee 

AAot*-?+ A,—A(n—2—r)A, -»,=90 (1 r n—3), 

KE —XAo; "+?+ By+A(n 1)A, 9 0), 
AAay"+*+A(n—1)A,,_,+(r7—2)A,_.—(r—1I)B,_, = 90 (Br < n-I 
(n—1)A| E—A-+AAo} | +-A(n—1)k’+-(n—2)A,,_.+(n—1)B , (). 
k’+- B_,—(n—2)AA,,_. = 9. 








Putting k’ 0 and solving these equations: 
’, ae 2AF ’ 
° a(xk+1) 
4 2[A2F(n—2—r)oi+AF ar *] 
. n| 1—)2r(n—r 2)| - 
 —— 2 2) i’ 
B, 2F[«—(n—1)d?|] 204 on” 
a(k+1) T 
9 9)\)\2F 
B . 2(n—2)A2F 
a(k-+-1) 
, 2AF ( n—-2 ar . 
_ K } o?-2, 
m(xk+1)\ n—1 > 
rB 2AF gritty 2| FA(r—1)a} 7 -1{ 1 +-A2(n— 1)} ] 7 


. 71 A2(n—r—1)(r 1)| 


2| Fo5 1f(n 1 )r—(r—1)?} | 


= (2 r n—2 
7 1—A*(n—r—1)(r 1)| 
The case n 4 
As an example we consider the case of a curvilinear square (1.e. = 4 
when there are isolated forces at the vertices given by 
F=Z at ¢ = i, iR at o= 4, Rata l, 
and iR at o ‘. in 
Then substituting in the above equations, 
2R o—] 2.R o—t 
Q, log log a 
7 o+] 7 o+i 
Q., Q, 2, Y A,o 1+ ¥ Avo, 
2/2 = 2) 
2(A2F'o,+AFz,) 20k 
where A, ——————- ol A, 
a(1—A?) a(«+1) 
so that > A,=0, and 5 A, = 0, 
_ — - 


») 
therefore Q “log 


o ) 21R o— 


o- 
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Also z'(a) = c(1—3A0-*), A < 4, so that nn being zero, the periphera 


stress ss is given by: 
238 = Q)'(c) z'(o) 1-()’(6) 2'(6) 
4R(1+-3A) 
ac(1— 6A cos 4n+9A2)° 
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This is symmetrical about both axes and has the value 4 R(1-+-3A)/ze(1 


at the points of application of the forces. See Fig. 8; A 


I wish to express my thanks to Mr. A. C. Stevenson for suggesting thes: 
problems and for many helpful suggestions and criticisms in their solution, 
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| { NOTE ON THE HODOGRAPH TRANSFORMATION 
FOR THE TWO-DIMENSIONAL VORTEX FLOW OF 
AN INCOMPRESSIBLE FLUID 
3A} By 8. GOLDSTEIN and M. J. LIGHTHILL 
(Department of Mathematics, The University, Manchester) 
r TNese ? . 
ution Received 11 October 1949] 
SUMMARY 
| In the hodograph transformation of the two-dimensional vortex flow of an incom- 
Roy. | ssible, inviscid fluid, branch lines arise, similar to those previously encountered 
the theory of the irrotational supersonic flow of a gas. The result is of importance 
the calculation of vortex flows with ‘free’ streamlines, and in particular for the 
— ulation of jets with vorticity present. As a simple mathematical example the 
» Pi lograph plane is described corresponding to the flow of a stream with uniform shear 
sta circular cylinder; the hodograph plane is a Riemann surface of six sheets. 


|, Introduction 

ly Kirchhoff’s classical solution (1) for the two-dimensional motion of an 
ompressible fluid in a jet issuing through an aperture in a plane wall, 
motion was taken to be irrotational, and the same is true of recent 
extensions (2). In a real fluid, however, a boundary layer will be formed 
mg the wall, and probably the most important effects will be due to 
e vorticity in the boundary layer and not to the viscosity of the fluid. 
\. T. Lord therefore began a study of two-dimensional vortex motion in 
ets. Solutions of such problems may have practical applications in the 
lesign of slots for sucking away a boundary layer. As a beginning, it was 
ught that a solution might be sought for a fluid possessing regions of 
nstant vorticity, or, as an easier alternative, for the case of a single 
neentrated vortex sheet near the wall upstream of the aperture and 
sing into the jet. As for problems of irrotational motion in jets, the 
dograph transformation would be used, the components, u, v, of the 
id velocity (or logg and #, where q is the magnitude of the resultant 
id velocity and # the angle the velocity makes with a fixed direction) 
eng taken as independent variables. It was found, however, that the 
transformation from the physical to the hodograph plane is not a one-to- 
ie transformation, and the hodograph plane is a Riemann surface. In 
e case of a concentrated vortex sheet the surface folds back over itself 

| ita winding point (as in analytic-function theory), but with distributed 
rticity it folds along a curve, on the whole of which the Jacobian of the 

| tansformation vanishes. These curves are ‘branch lines’ of the same type 


(Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 3 (1950)] 
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as those previously encountered in the theory of the irrotational supersonic 
flow of a gas (3, 4, 5). Because of the occurrence of branch lines, it will 
not be a simple matter to use the hodograph transformation for the study 
of flows with vorticity, and problems with ‘free’ streamlines will be diffi- 
cult. Of the cases previously mentioned, the problem of the jet with a 
small constant vorticity is manageable, and has been studied further by 
W. T. Lord; the case of a concentrated vortex sheet may also be manage 
able, and will be studied further. 

If x and y are rectangular cartesian coordinates, and uw and v the corre 
sponding components of the fluid velocity, then for a two-dimensional 


motion, with vorticity «, of an incompressible fluid, w and v satisfy the 
equations 

cu ' ov 0. ov cu ie 

COX cy Cx oy 
When ¢« = 0 any singularities of the transformation from the (x,y) plane 
to the (u,v) plane are isolated points; when ¢ is not zero, branch lines (on 
which ¢(u, v)/e(x, y) 0) may occur for any value of ¢«, no matter how 
small. 

In order to demonstrate the occurrence of branch lines in the hodograph 
transformation of the two-dimensional vortex flow of an incompressible 
inviscid fluid we work out in the next paragraph the hodograph plane 
corresponding to the flow of a stream with uniform shear past a circulai 
cylinder. This is a simple mathematical example, of no physical interest 
in itself, intended purely to exhibit the kind of mathematical difficulty 
that may arise in using the hodograph transformation for other problems 
involving vortex motion. 

In his investigation of the geometry of the hodograph transformation, 
Craggs (5) showed that branch lines, defined as lines along which the 
Jacobian of the transformation vanishes, cannot occur in irrotational 
subsonic flow, and it has not previously been pointed out that they do 
occur in vortex flows. Craggs also showed that a folding back of the 
hodograph plane is usual along such lines, and that in general branch lines 
are cusped in the hodograph plane; these statements are verified in the 
example below. 


2. The flow of a stream with uniform shear past a circular cylinder 
With rectangular Cartesian coordinates (x, y), let the undisturbed values 
of the velocity components w, v of a stream of inviscid incompressible fluid 


be u U—ey, v = 0, where ¢ is a constant. Let the stream flow past 
the circular cylinder, x?+-y? 1, of unit radius (Fig. 1). The motion 
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rsonic ; . . P , ‘ 
anger Helmholtz’s theorems; the fluid being incompressible there will be a stream 
it wil ; x i o 
finetion Ww. such that u Cus /OCy, v us /Ca, which satisfies V7 —€, 
study 
. dite Y 
e cliff —s 
with a 
} ' U 1 
her b U-ey 
=» 3 
nage a 
—————-> 
—_—}> 
Corre 
isiona Fic. 1. 
fy the nd which is constant on the cylinder. The solution giving the correct 
disturbed motion is, with r and 6 as polar coordinates, 
Uy _ 2 
Uy cy Wee et ae i. 
‘ x“-+-y¥* (2*--Y*)* 
plane 
es (on | Ul) |sin 6 ter? ie(r2 3) 008 28, (1) 
r how | 
U ‘ oe 
| Cu cos 26-+-er sin 8 - sin 34, (2) 
. ye D3 
grapl 
rea] l . € 
ssibli » sin 24 . = COS 36. (3) 
plan . rr 
reular The branch lines are the lines along which @(u, v)/é(x, y) vanishes. Now 
terest Clu. 2) l clu. v) 402 UU. G6eU . , Qe?  3e? 
! : sin 39 + —— sin @— —— 4 cos 40. (4) 
iculty clr. wy) r olr.@) ry? aa y 4r° 2rt 
blems| The curve é(u, v)/ (a, y) 0) is a real curve. To proceed further, assume 
that e« is small; then 7 must be large on the curve whose equation is 
ation : 
spproximately er ; J 
wre PI ~" r? sin 36 2U /e, (5) 
tiona ind on which 7 (2U’/e)). The curve has three branches, lying between 
ev d the asymptotes @ = 0 and 6 = 7/3, 6 = 27/3 and @ = a, and 6 = 47/3 
f the and @ 57/3, respectively (Fig. 2). 
1 lines To the same approximation, the last terms on the right in (2) and 
in the 3) may be omitted. 
In the hodograph plane take coordinates X, Y, defined by 
Y l u, Y v. (6) 
inder Then from (2), (3), and (5) the parametric equations of the branch line in 
alues| the hodograph plane are, to the same approximation, 
» fluid 2] 
aa X = e(— }* (sin +3 cos 28 sin 38), (7) 
ed esin 30 ’ 
otio! n 36\2 
’ €sin 30\4%,, . 
c, by . }° U sin 20. (8) 
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In 0 < 6 < 7/3, X goes from 0 to +c and Y rises from zero and falls to 
zero again. In 27/3 < 6 < a, X goes from +c to 0 and Y falls from zero 


and rises to zero again. In 47/3 < 0 < 52/3, X goes from © back to 


©, remaining always negative, and Y rises from zero, falls through zero 














rie. 2. 


to a negative value, and rises to zero again. dY /d@ vanishes when sin 56 = 0, 
the relevant solutions for @ inside the above ranges being 
a 40 Tn 87 
= -,—,—,—. (9) 
& 56 &° $ 
dX /d@ vanishes when cos 36 = 0 or sin5@ = 0. The zeros of sin 50 given 
in (9) are therefore zeros both of dX/d@ and of dY/d@, and are cusps. 


X Y = 0 is also a cusp. X as a function of Y has stationary values at 
those zeros of cos 36 for which 6 = 7/6, 57/6, and 37/2, and it is easily 


shown that the first two of these are maxima and the third a minimum. 
From these results the branch line in the hodograph plane may be sketched 
(Fig. 3). There are five cusps. 

The mapping of the physical plane on to the hodograph plane may be 
studied by considering the behaviour of lines 6 = constant, which (from 
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6), and from (2) and (3) without the last terms on the right) are the cubics 


sithout points of inflexion in the hodograph plane whose equation is 
. o(U sin 20\3 
X Y cot 26+. sin | y J’ (10) 


lhe details of this investigation are omitted; the results are as follows. 
The hodograph plane is a Riemann surface of six sheets. Three sheets. 
‘ich we shall call sheets 1, 2, and 3, do not contain the cusped figure 

RKFHG: the other three sheets—sheets 4, 5, and 6—do not contain the 
sped figure OA BDC. 


Y 
A 





py 





Y! 


Fic. 3. 


Sheet 1 is the whole plane, cut along 4B and CD. Sheets 2 and 3 each 
nsist only of the interior of OA BDC, cut along OX. The right-hand 
lges of the cuts 46 and CD on sheet 1 are joined to sheet 3. The left- 
ind edges of these cuts are joined to sheet 2; thus at these joins the surface 
Ids back over itself. Sheets 2 and 3 are joined along OA and along OC; 


thus at these joins the surface again folds back over itself. 


Sheet 4 is the whole plane, cut along X’OX, EF, and GH. The upper 
lge of the half-cut OX on sheet 4 is joined to the lower edge of the cut OX 
nsheet 3, and the lower edge of OX on sheet 4 is joined to the upper edge 
n sheet 3. Sheets 5 and 6 each consist only of the interior of EF'GH, not 
tatall. The left-hand edges of the cuts E F and GH on sheet 4 are joined 


to sheet 6. The right-hand edges of these cuts are joined to sheet 5; at 


these joins the surface folds back over itself. Sheets 5 and 6 are joined 


long FH, which is another fold. 
The cuts OX on sheet 2 and X’O on sheet 4 are impassable barriers on 


the Riemann surface. Neither corresponds to any point at all in the 


| hy sical plane. 
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The part of the plane inside each cusped figure is covered four times by 
the Riemann surface. The rest of the plane is covered twice. 

The complexity of the hodograph plane in this simple example is a } REF 
warning of the difficulty that may be encountered in using the hodograph | , \ 
transformation to study flows with vorticity present. 
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isa) REFLECTION AT A LAMINAR BOUNDARY LAYER OF 
zraph | s WEAK STEADY DISTURBANCE TO A SUPERSONIC 
STREAM, NEGLECTING VISCOSITY AND HEAT 
CONDUCTION 
By M. J. LIGHTHILL 


(Department of Mathematics, The University, Manchester) 





Received 25 February 1950] 


SUMMARY 


\ perfect gas, with constant specific heats and zero viscosity and thermal con- 
tivity, is considered to flow in a steady two-dimensional pattern which is only 
ghtly disturbed from a basic flow parallel to a solid plane wall. The Mach number 


the basic flow varies continuously and monotonically from zero at the wall to a 


ersonic value at a distance 6 from the wall, beyond which it takes this supersonic 
ie uniformly. The disturbance is considered to be due to a weak stationary 
ple plane wave incident on the ‘boundary layer’ from outside. The streamline 


ttern and pressure coefficient distribution, resulting from a given incident wave, 
found to depend only on the Mach number distribution in the basic flow (not 


the individual velocity and temperature distributions). Any singularity in the 
: | : ' ‘ 


ident wave is propagated along characteristics in the supersonic part of the layer 
reasing in strength) and is reflected from the sonic line modified in character. 
is a discontinuous increase (in pressure or pressure gradient) is reflected as a 
sitive logarithmic infinity of the quantity (which result of a linearized theory may 
rrespond to a ridge in the real values). The upstream influence of the disturbance 


the layer is shown (for all Mach number distributions) to be completely negligible 

1 distance of one layer thickness ahead of the point of incidence of the wave 
nt on the outside of the layer. (It is explained in section 8 why this conclusion, 
the limiting case of a discontinuous Mach number distribution considered by 
n and Finston (1), is at variance with their results.) For this reason the theory 

ist be completely inadequate as a theory of shock laminar boundary layer interac- 
n,and the significance of this conclusion is discussed in section 1. Four radically 
ferent possible mechanisms for the interaction (of a more complicated nature) are 
suggested, none of which ean be discarded in the light of the data known to the author. 
Some readers may be assisted by the concise account in section 4 of the more useful 


sults of Langer’s theory (2) of differential equations with a large complex parameter. 


|. Introduction. Reflection of shocks from laminar boundary 
layers 

Iv is fully established (3, 4) that, in steady supersonic flow, a wall with 
minar boundary layer attached behaves, for purposes of shock reflection, 
ery differently from the solid boundaries of the ordinary inviscid theory. 

Some diffusion of a disturbance, in the process of reflection by the laminar 


[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 3(1950)] 
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layer was, of course, to be expected. But the experimental results, that} yiscou 
upstream influence of the shock extends a distance ‘tens of boundary layer} Mach 


thicknesses’ ahead of the point of incidence, and that a shock is reflected 


it cau 
locally as an expansion wave, are somewhat surprising. Given those results except 
the others accord with previous concepts. The boundary layer thickness} |arge 
increases ahead of the shock (over and above any normal thickening by} magn 


diffusion of shear) to a peak at the point of incidence, whence it declines} the la 
(or at least ceases to rise). As would be expected, compression waves from] the s 
the region of thickening are found to combine with compression waves comp 
from behind the expansion wave to annul the latter, and form two shocks] the s 
(as envelopes of Mach lines) which join into one (the reflected shock of the} other 
ordinary theory) at a distance of order a hundred boundary layer thick Th 
nesses from the wall (4). But this effect can naturally be erased by convex} jn th 
wall curvature. The fundamental phenomenon is that the reflected wave | impr 
consists of a rather sharp expansion in the midst of a band of compression,}  subs« 


which band, especially upstream of the point of incidence, is remarkably] sheet 





wide. the v 

Two radically opposite views as to the correct basis for an explanation | (dist 
are held; neither basis has yet been developed into a full theory. Lager-| ofth 
strom, Cole, and Trilling (5), following Liepmann, Ashkenas, and Cole (4), the e 
see the problem as a purely supersonic one, bound up intimately with the distu 
action of viscosity near the base of the shock, whose non-zero thickness| Nx 
needs to be considered. They magnify the possibilities of ‘interaction| influ 
between longitudinal and transverse waves’ in viscous compressible flow,} mak 
citing as an instance the known flow pattern for supersonic flow past) Fins 
a blunt-nosed plate parallel to the stream, which, however, is usually} This 
regarded as fully explicable on adiabatic theory (except in so far as} Thier 
boundary layer separation may be involved for certain nose shapes). the’ 
They construct a linearized theory of viscous gas flow without heat con- 
duction, which replaces the wall boundary condition of zero velocity by 
one that the velocity take a constant value slightly less than that of the 


stream. This they propose to apply to the interaction problem. (1 
On the other hand, Tsien and Finston (1), following Howarth (6), hold 

that the subsonic portion of the boundary layer is the chief factor governing | 

the anomalous shock reflection, which is usually absent for a turbulent 

layer (when the reflected wave is only displaced by a distance of the order (i 

of the layer thickness, and is not altered in character) because for this 


the region of subsonic mean speed is usually microscopic. They expect 
viscosity and heat conduction to affect only slightly the local modification 
of the laminar layer by a shock, and this view has more reason behind it . 
than might first appear. The air in the boundary layer is not ‘more | soni 
0 
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viscous’ than that outside; it merely has a greater vorticity and lower 
Mach number. Viscosity takes time to influence any fluid element—thus 
itcauses the layer to grow only slowly at the higher Reynolds numbers— 
except inside a shock, but here the velocity gradient is quite extraordinarily 
large and, on the other hand, the overall change is uninfluenced by the 
magnitude of the viscosity. The shock thickness is small compared with 
the layer thickness except where its strength tends to zero on approaching 
the sonic line (more precisely, where the reciprocal of its strength is 
comparable with the Reynolds number based on layer thickness, i.e. with 
the square root of the Reynolds number based on surface length); on the 
ther hand, here it is a less significant phenomenon. 
The theory to be given here is a modification of that of Tsien and Finston, 
n the direction that they themselves suggest as likely to give greatest 
mprovement. They represent the boundary layer as a region of uniform 
subsonic flow, thus replacing its vorticity by two concentrated vortex 
sheets, one at the interface with the supersonic main stream and one at 
the wall. Here the same work will be undertaken with a layer of finite 
listributed) vorticity, giving continuous decrease of velocity from that 
fthe main stream to zero at the wall. As in the work of Tsien and Finston 
the equations of plane, steady, adiabatic flow are used, and the squares of 
listurbances to the parallel flow neglected. 
Now the results found below, relating to the distance of upstream 
influence of the shock, are in complete disagreement with the experiments, 
aking this distance far too small, and much smaller than in Tsien and 
Finston’s theory. Yet the model of the flow has undeniably been improved. 
This indicates a possible element of chance about the good agreement of 
lsien and Finston’s work in this respect, but does not necessarily annul 
the validity of their general concepts. Four possibilities exist: 
(i) Lagerstrom, Cole, and Trilling are right in ignoring the subsonic 
layer and ascribing the results principally to viscous diffusion in 


the supersonic layer; 07 


the theory of this paper could be corrected, by taking viscosity 
into account after the manner of the theories of boundary layer 
stability, especially near the wall where the full viscous boundary 
condition could then be applied; or 

iii) viscosity is not important, and the present theory is correct for a 
really weak disturbance, but the pattern changes rapidly as the 
shock amplitude is increased, so that to predict the effect theoreti- 


cally the full non-linear adiabatic equations would be needed;* or 


* This view appears less improbable when the focusing of disturbances as they near the 


sonic line (section 6) is taken into account. 
5092,11 
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(iv) there is a bubble of dead air at the wall near the base of the shock 
(perhaps slightly in reversed flow) of considerable length but small 
breadth; this has been formed by repeated separation until the 
pressure distribution was such as to permit attachment of the main 
laminar layer (which is outside the bubble). 

Hypothesis (iv) as well as (iii) would imply considerable variation in the 
distance of upstream influence with shock strength; the evidence is 
negative on this point but not conclusively so. If the same qualitative 
behaviour is found for very weak shocks as for stronger ones, then the 
fundamental explanation must be (i) or (ii), and the effects described in 
(iii) and (iv) can only modify the phenomenon gradually as the shock 
strength is increased. But this is not yet certain. The author proposes to 
investigate hypothesis (ii) mathematically, and also (unless stronger 
evidence against them becomes available) hypotheses (iii) and (iv), but 
decided first to publish the present investigation without delay as an 
essential datum to the understanding of this much discussed topic. 

The evidence of the theory of this paper regarding the local reflection 








of the shock by the laminar layer may possibly be of value even though 
that relating to the distance of upstream influence is not. This evidence 
is that it is from the sonic line that the shock is reflected as a rather sharp | 
disturbance, and that the latter is in the nature of a pressure ridge, such | 


as might be produced if the boundary layer turned rapidly but con- 
tinuously through an angle to produce a band of compression waves, and 
then turned instantaneously back in the opposite direction (giving a 
Prandtl—Meyer expansion). The latter at least is observed (see above). 
What compression waves there are ahead of it interact with the 
expansion further from the wall (this could not be reproduced by the 
linear theory) and help to annul it and form a shock as an envelope of 
characteristics. 

Unless hypothesis (ii) is the correct one, the theory of this paper is at 
least of interest as elucidating the reflection of a smooth weak disturbance. 


| Note added in proof. Experiments conducted by Mair and Bardsley at 
the Manchester University Fluid Motion Laboratory have now established 
that the conclusions of the present theory are correct for the reflection of a 
weak shock by a turbulent boundary layer. At a Mach number just less than 
2, any shock, weak enough to produce a deflexion < 5°, is found to be reflected 
locally as a pressure ridge, shown on the Schlieren records as a black line 
immediately adjacent to (and ahead of) a white line; the latter (representing 
the expansion) appears equally strong where both emerge from the boundary 
layer, but becomes weaker further out by interaction with the compression 
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waves behind the ridge. There is no visible upstream influence or boundary 
layer separation; and thus all the predictions of the theory agree with 
observation in this range of strengths. (With stronger incident shocks the 
reflected shock and expansion wave lie further apart, and the reflected 
shock springs from upstream of the point of incidence, and is so strong 
that the layer separates.) 

It is not perhaps surprising that the theoretical conclusions are un- 
affected by the turbulence in the boundary layer, as one might expect 
the reflection of a stationary shock to depend principally on the mean flow. 
The significant deduction appears, therefore, to be that the present theory 
iscorrect for velocity profiles typical of turbulent layers, but incorrect for 
those typical of laminar layers. Preliminary calculations, based on hypo- 
thesis (ii) above, support this conjecture, for they indicate that the low 
Mach number layer very near the wall can probably send disturbances far 
ipstream by a mechanism involving viscous stresses; now this inner layer 
yould presumably be of negligible thickness in the turbulent case—which 

ight be expected to nullify the effect. ] 


2. Mathematical formulation of the idealized problem 


Let a perfect gas, 


gas, with constant adiabatic index y and zero viscosity 
nd thermal conductivity, flow in a steady two-dimensional pattern, 
the 


Let the pressure p have 


disturbance from 
QO. 


vhich differs only by a small a basic flow in 


direction parallel to a solid plane wall y 
,uniform value P in the basic flow; and let the velocity vector (wu, v) and 
the density p therein be functions {U(y), 0} and R(y) of y alone, which are 
Let A(y) fyP/R(y)* be 
the value in the basic flow of the local velocity of sound a = (yp/p)', and 


ntinuous for all y and constant for y > 6. 


suppose that the function M (y) U(y)/A(y) is non-decreasing in 0 < y <6 
fom a value M(0) 0 (where M’(0) > 0) to a value M(8) M, > 1. 





\eccents are used to denote derivatives. 
Then the equations of the disturbed flow, expressing constancy of 


itropy along a streamline and conservation of mass and momentum, are 


Cc ; Cc ; »f p ( ) C ; C 
ut pit alu yl). (pu) +-— (pv) Q, 
Ca oY Ca cy CXL C y 
(1) 
Cu Cu , Op ov cv . Op 
pu pi - 0, pu + pu —-+t 0. 
Cx cy Cx Cx cy cy 


-U(y), v, p—Ry), p—P, 


Vhen linearized, on the assumption that u 
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a—A(y) are small, equations (1) become theorv 
: . able tu 
: Cp = . Cp . ‘i 
U(y) u AMy)(L (y) cf + R'(y)v), The u 
Ox Cx 

theless 
U(y) P+ Ry) = + — {Rye} = 0, | discor 
Cx Ox «oy be hoy 

‘ t,o 
R(y)U(y) - : + R(y)U'(y)vs = = @, - 
Ox Cx functl 


Ov Cp 


R(y)U(y)< 0. 
ox oy 

the be 

Now elimination of @p/éx and éu/éx from the first three of equations (2 

gives 
a op ; CO ta, id L(y) ép ; 7 

R(y)U'(y)v4 U(y)— {Ryo} + —S = - U(y)R'(y)v, (3) provi 
Ox cy A*(y) 6x 
tinuit 
or — {1—M*(y)} = Riy)U%y) = ( Ms } (4); impo 
Ox cy\U(y) conti 


But actually R(y)U2(y) = yPM%(y). This suggests that the transforma. ‘Ht 


tions phras 
o—P a v 7 (5 incon 
UD, @ 7) 
> y a 
vi U(y) lo 
the v 


be made (so that @w is a pressure coefficient and 7 is the approximat 


5 : ; .  Hene 
deflexion—i.e. angle made by a streamline with the x-axis); then equatior ; 
. - upst! 
(4) and the last of equations (2) become 

is Cw e on Cw a con ; It w 

f1— M*(y)} M?(y) —, —— = — Sy) —, (6 ; 

Ox oy cy Cx wy 

It 


Thus the distribution of pressure coefficient @w, and the streamline the « 
pattern (which is deducible from the distribution of deflexion »), depend jnye: 
not on the separate velocity and temperature profiles in the layer, but for y 
only on the distribution of Mach number M(y). The equation for w alone (9) t 


is obtained, by eliminating 7, as the 
02 J. 2) ’ ( A o2 t - 
: * r ty) oe mT Mxy)r° nd 0. (7 
ey? M(y) cy ~ 


It is considerably simpler than that for » alone. 

The cause of the disturbance will be taken to be a simple plane wav 
incident on the layer 0 < y < 8 from outside. In it the pressure coefficient 
a@ must take the approximate form f(x«+8y), where B = (M?—1)!, by the 
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theory of small disturbances to a uniform supersonic stream. It is reason- 
ble to suppose that if f is uniformly small @w must be uniformly small. 
fhe use of (7) is therefore justified unless the derivatives of w are never- 
ltheless not uniformly small; they probably are not if f possesses a 
liscontinuity (as would be needed to represent a shock); however, it will 
be hoped that some information can be obtained by solving equation (7). 
The reflected wave in y > 6 must’ take the form w = g(x—fy), and the 


tion g is principally what we want to find. Since 


os fla Sy q(a By) for y oO, 


e boundary condition at y 0 1s 

tions Des Dont p 

3 | 26f (x Bd), (8) 
CA CY | 

rovided that the derivatives are continuous across y 5. But discon- 
muity of @ or of 7» across an approximate streamline like y 5 is 
(4 possible by then physical significance. Hence CT CL and con cx are 
tinuous across Y¥ 6, and since .V(y) is continuous, equations (6) show 
lleaies Ca cy and én cy are also continuous. This argument has been 


wed to avoid assuming .W’(y) continuous at y = 6, which might be 
mvenient when model velocity profiles are chosen. 


To complete the formulation of the problem the boundary condition at 


the wall y 0 appropriate to an inviscid fluid is needed. This is v 0. 
Oximate ,. oe ; ’ ; - eels 
Hence .W(y)en cy 0, so by (6) Cw/cx 0. Since w 0 at infinity 
(uatlol 3 ’ . 
nstream this gives 
—_ ( 
oo), 9 = 0. (9) 
It will be seen later that this with equation (7) implies that wa = O(y*) 
y—> VU.) 
It is the boundary condition (9) which may, in a sense, be the cause of 
eamline the disagreement with experiment. Its modification would be the aim of 
depen estigation (ii) of section 1. But it is incontestable for an inviscid fluid, 
ver, but for which therefore equations (7) must be solved under conditions (8) and 
aw alone (9) to get the pressure field, whence the reflected wave will be derived by 
tne equation : 
CT Cw | K 
|p 2Bq'(~— sd). (10) 
_ | C2 7 
3, Solution as a Fourier integral 
If f(a) can be expressed as a Fourier integral 
1e wave Q 


efficient 
' by the 


f(x) ( etkr F(k) dk, (11) 
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it is reasonable to seek solutions for w and g in the form 


the ze 

(x,y) [ ek] T(k, y) dk, g(x) = [ eke G(k) dk. (12 om 

a - upstre 

Then equation (7) becomes an ordinary differential equation for II, 4, Ré 

val 2M"(y) dil k2{M%(y)— 111 = 0, (13) The 

dy? M(y) dy form | 

and the boundary conditions (8) and (9) become first f 

IT, (k, 8)+- ikl (k, 8) 2Bike**P> F(k), H(k,0) = 0. (14) A bri 

as the 

Choose Q(k,y) to be any solution of equation (13) with Q(k,0) = 0 te 
Then by (14) 
I(k, y) = 2Biket*83 F(t) Oh y) . (15 

Q,(k, 5)+- BikQ(k, 5) 
Further, by equation (10), et 
: e P\ 
G(k) = [BikIl(k,8)—I1,(k, 8)]/2Bike-*85 Let 4 


e2xikBd PF I) BikQ(k, 5) 0,(k, 5) Then 


} : ~— (16) ¢ a 
BikQ(k, 8)+Q,(k, 8) for a 


Actually, on the assumption made in section 2 that .W’(0) > 0, th 
point y = 0 is a regular singularity of equation (13), at which the expo 
nents are 0 and 3. Therefore the simplest form of Q(k, y), and that which 
will be adopted, satisfies 

Uk,y)~y as yd (17 
for all k. In terms of this solution of (13), which (17) defines uniquely,, wher 
equations (15) and (16), with (12), provide a complete formal solution t 
the problem of section 1. H 

Two uses will be made of this formal solution. In section 5, by studying! pow 
the asymptotic behaviour of Q(k, y) as |k| > 00, those of II(k, y) and Gk) used 
will be found. These will be used in section 6 to find what singularities in w js al 
and g are produced by a specific singularity, such as a discontinuity, inf. and 
For it is well known how the singularities of a Fourier integral at finite 
points can be deduced from the asymptotic behaviour of the integrand. 
(This investigation also enables the convergence of the Fourier integrals| !° 
(12) to be demonstrated. ) 

Secondly, in section 7, the zeros of the common denominator of (15) and, > ‘ 
(16) will be considered. Since the integrals (12) can each be expressed in 
terms of the residues of the integrand at poles with im k > 0 for x > 0, and 
in terms of those at poles with imk < 0 for x < 0, knowledge of the zeros }, 
of the said denominator with least |imk| will enable the asymptotic dis 
turbance to the boundary layer for large |2/8| to be found. In particular, 
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the zero —ix,, with numerically smallest negative imaginary part, will 
determine the asymptotic form (proportional to e“'* as x» —oo) of the 


(12 . 
| upstream disturbance. 
IL 4, Résumé of Langer’s theory 
(13 The asymptotic behaviour of Q(k,y) as |k| > 0 (which is of different 
| form in the subsonic and supersonic regions) will be studied by the method 
first fully developed by Langer (and foreshadowed by Jeffreys and others). 
(14| A brief statement of the relevant conclusions of this theory is first given, 
1s these are somewhat hard to extract from Langer’s original paper (2). 
() 
In the equation 
dzu du (7.2 ) 
(15 pat PY) | egy) +r(y)ju = 9, (18) 
( y ( ¥ 
t ply). g(y), and r(y) be twice differentiable* in the interval yy < y < yp. 
tay) << Oiny, < y < y, andq(y) > Oiny, < y < Yo, while q'(y,) > 0. 
Then two asymptotic solutions (valid as |k| > 00 uniformly in yy < y < Ye 
H 
for all arg /) are 
uv 
0, the / , 4 \ : . 
ex] us(y) = exp(—3 | p dy)|s'\q\H(kes)\_(ies),/(22), 
7 ? 0 . 
y ( 
whicl . (19) 
Uy(Yy) exp| Lip dy} sig —!sgn(y—y,)(ks) tJ, (ks), (27), 
(17 
iquely here S | vgqdy~|q'(y)b3(y—my)? as yoy, (20) 
tion ft 
Here s is imaginary for y < y,, but (ks)-’J,(ks), involving only even 
idying} powers of s, is real, and also independent of which square root of q is 
d Gh sed in the integral for s; hence both uw, and u, have these properties. It 
‘Sin isalso differentiable with respect to y, even at y = yj. But also |s|*\q|~! 
y, Mj. and s '\q -tsgn(y—y,) are differentiable even at y = y,, and as y> y, 
finite 
— 8 la’ (y,)]%(2)8| y—y, |4 2? 1 s\8 2 3 
orand ae M1 peAaV IT ? | \*. -sen(y—y,) ~ (— ®a(y), 
egrals Lg GP i¥— 41" Od (Y1) : 3q (Y) 
(21) 
5) and , that 
2 1 23, /(277) , 
sed it us(Y) | )° ad Us(Y3) 0. u(y,) = 9, 
3q'(y;)} (—4) 9 
0, and JI 3 ee (22) 
, 2\8/,’ “4H (277) 
’ ZETOS MalYa) = (ST (Ys —G yy 
ic dis ill 
niat * The theory may be extended to the case when at a finite number of exceptional points 
ic i 


t at y,) the functions are merely continuous, not twice differentiable. 
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When |ks| is large (e.g. for fixed y ~ y, as |k| >) the functions 
2 J, (2), (277) in (19) may be replaced by their universal asymptotic formula 
as |z| > 00, namely 
iz iz 


Cc é 


(izyrrd | ( iz) +h? 


which is valid for all argz provided that the principal value of both 
(v-+-4)th powers be taken. Hence, as ks| > ©, 








/ uy \! git etks e- tks 
u,(y) ~ ex ( 1 ”) dy) Isc ToxTe 
WY) ~ XP\—3 } POY Tl (kaye t (—ikey 
; (24) 
° s 'son(y Y1) eiks e~ tks 
u(y) ~ ex | 1 »dy : : : - af : 
2 | 2 | I e q , (iks)* ( iks)* 





Thus for large |ks! the general solution is some combination of solutions 
of the form 


exp| 4 | Pp dy) q iettks | (25) 
W 
For real k these are oscillatory when y y, and exponential when y < 4. 


[f a special solution is determined by some boundary condition in either 
of these sub-intervals, then the asymptotic form of this boundary con 


dition as |k| + oo will determine what combination of the forms (25) is 


required. To permit continuation into the other sub-interval (or into the 
1 


intermediate range where & but not ‘4s is large) expression of the 
resulting asymptotic solution as a combination of the forms (24) is needed 
(being careful to take principal values of the (4)th and (?)th powers inside 
the square brackets); when this is done the special solution has been 
expressed in terms of uw, and uw, and so knowledge about it for large |/| is 
complete. 

For example in section 5 a solution will be found to have the asymptotic 
form as |k| + o, for y < y,, 

y 
“um exp| ) | p dy) q — cosh| k(a—|s!)], (26) 
"1 

where co > 0. If it is asymptotically f,(/)u,(y)+f.(k)us(y), uniformly as 
k| > oo, then by (24) 


s kis eks ‘ 


ek é ks . 
i(k o(k ~ cosh[k(o—}s|)], (24) 
Sil | ar 5 Je | , ma cosh| k(o )| 
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rallargk. The solution is easily verified* to be 
| ie l . ‘ , 
~ —. {ekokt + e-ko(_k) folk) ~ —j {ekekt+-e-*o(—k)§}. (28) 
2V3 2V3 
fhe behaviour of uw as |k| > oo for y y, is then deduced, from (24) with 
sitive 8, as 
1~ exp , | p dy)4q e'** cosh(ka— fat) ‘ks cosh(ka-+-}7i)], (29) 
rall arg kh. 
5, Asymptotic behaviour of (, G, and II 
Equation (13) does not satisfy the conditions of the theory of section 4 
the whole interval 0 y 5, since p(y) 2M'(y)/M(y) is not 
fferentiable (not, indeed, finite) at y 0. But if € is a small positive 
mber the conditions are satisfied in « y <— 6. The y, of section 4 must 
the solution of M(y,) |. Thus y, is the thickness of the subsonic layer. 
| when « y y,, for example, the asymptotic form of Q(k,y) as 
© must be some combination of the forms (25), which for equation 
become 
uf 
M (y)| 1—M*(y)|-*e#**, s= | [M*(y)—1]! dy. (30) 
But to deduce what combination is required, from the boundary con 
tion (17) for Q(k, y), information about the asymptotic behaviour right 
wn to y = 0 is needed. This can be obtained from the fact that the 
ond-order equation for d?Q dy? does satisfy the conditions of Langer’s 
heory near y = 0 (as will not be surprising to the reader when he considers 
t in the general solution to (13) the values of d?II dy? and d3I1 /dy? are 
bitrary, indicating an ordinary point of the equation for the second 


erivative). Therefore the combination of the asymptotic forms (30) 
juired to represent Q(k, y) must represent it correctly down to y 0 as 
is coefficients of powers of y higher than the first are concerned. 

Thus, by (17), the coefficient of y° in the asymptotic form must tend 


l, and that of y? must tend to 0, as |k| > oo. This fixes it as 


M (y)[ 1—M*(y)|-! cosh h ( {1—M*(y)|* dy /| sk?M'(0)], (31) 
we /i2 
that, as |k| + 00 for « y y, and all arg k, 
Lk? M (O)Y k, y) ~ M y)| 1 M*(y)| ‘ cosh| k(o 8 )I, (32) 


\3, whether the upper or lower signs be taken. 
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as in equation (26) above, where 


v1 


o | | 1— M?(y) |! dy s(0)!|. (33 


0 
The transition to the asymptotic behaviour of Q in the supersonic 
region has been carried out in section 4. By (29), as |k| > 00 fory, < y <§ 
and all arg k, 





4k2M'(0)Q(k, y) | 
~ M(y)| M*(y)—1]-*[e** cosh(ko— 4) +e-"* cosh(ko- t7t)| (34 | 
M(y)| M?(y)—1]-'[e*? cos(ks—4)+-e-*? cos(ks+ 47) |. (35) | 


Forms (34) and (35) are both useful. For real & the latter is split up into 
a part important only for positive & and another important only for | 
negative k. 

The common denominator of (15) and (16) can now be asymptotically 
estimated. Since 8 (M?—1)! (ds/dy) 
equation (34) gives 


and putting (s) S, 


y 5 y=o0 


$k?.M'(0)| Q,,(k, 8) + BikQ(k, 8)] ~ M, B-*2Bike™S cosh(ka—i). (36 


Here similar terms, involving e+**S but no factor k, have been discarded; 
this is valid (since the terms retained are of higher order) except when 
the imaginary part imk of & is large and positive, when e-‘*S cannot be 
neglected in comparison with ke’. Similarly for the numerator of (16), | 


when imk is not large and negative, 
$k?M’'(0)[ —Q,(k, 8) + pikQ(k, 8)] ~ MM, B-2pike-*S cosh(ko+}ai). (37) 
Thus for bounded im, and in particular for real k, as |k| + 00, (16) becomes 


.cosh(ko- 


: 1 : 
G(k) is e2tkBd Fk) ¢ 2ikS ! ) m prik(ps SF(k)i sen k. (38) 
cosh(ka— 477) 


As for II(k,y), when |k| +o for bounded imk and ¢« < y < y, by 

(15), (32), and (36), 
I(k, y) ~ e'*82-S) F(k)M; 181M (y)| 1—M?(y)]-! 
O{F(k)e—"**}, (39) 


cosh| k(o—|s|) | 


so that, if the integrand in (11) is bounded at k = -.oo, then the integral | 
(12) for w(x#,y) is convergent by a large exponential factor, indicating 
regularity as a function of x. 


cosh(ka— 477) } 
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| But as |k| +o for real k and y, < y < 34, by (15), (34), and (36), 
| : , oer - 

(ky) ~ F(k) My 1B M (y)[ M?(y)—1] if etk(B5+s—S) + eik(Bd s-Siisonk]. (40) 


(33 . , ‘ 
| For ¥ = Y; or, more generally speaking, whenever |k| is large but |ks| is 
t, neither asymptotic form of II is valid. Then it is necessary, in 


Psonic | estimating $k2.M’(0)Q(k, y), to use its expression in terms of the functions 
Y 19) with the coefficients (28). In particular, by (22), 


bk? M'(O)Q(k, yy.) ~ fylk)uy (yy) 


| 9 194 (9 
ie I fpkopt1 ko, mr rie \*5 ~_ 
3.2M'(y,)} (—4)! 


7 ka Jed ’ ; 4] 
(Oe {2M (y,)}*34 3)! , | 


p into| s that, as |k| + 00 for real k, by (15) and (36), 
ly fol 943 5 eo . 
. . ik(Bd—S) F/ ]-\ 14 = »\tpt7isgnk ) 
(hk, y,) ~ « Fk M; p M's ya “1)! k et ‘ (42) 
HE (Wy )s°(—3)- 
tically 
N 
] 
6. Propagation of a discontinuity 
(36 If the incident wave f(a+By) possesses a singularity, such as a discon- 
tinuity (as would be needed to represent a shock), this will be associated 
urded 


vith a specific asymptotic form of F(k) as |k| + 00 for real k. This in turn 
wher vill produce specific asymptotic forms for G(k) and II(k,y), by (38), (39), 


ot be 10), and (42), which must correspond to singularities of g(a—fBy) and of 


(16 . . 
r (1 a(z,y) which can be determined. 
Thus if f(a) includes a discontinuous increase A at x = 0 (corresponding 
»-, | tashock of ‘strength’ yA on the line x+By = 0), then F(k) includes a 
(37 ; ' ' ; ; 
tem —iA/2zk for real k of large numerical magnitude. It is only asymp- 
comes totic to this if the discontinuity is the ‘worst’ singularity for finite x. But 
suming this, . 
: iA ‘ 
- F(k)~ (43) 
(vs 2ark 
s k|-> oo for real k. 
a 
yy, by Hence by (38) G(k) ~ e2#Pe-S) : (44) 
‘ 2m\k 
(he exponential factor merely indicates that the singularity of g(~a—By) 
' ison the line «—By+2(B5—S) = 0, which will be explained physically 
(39 


elow. But the factor 1/(27\k|) indicates that the character of the singu- 


— wity has changed, and that 


cating 


A 
g(x— By) log| x By + 2( 


7 


83— S)|+-continuous function. (45) 

















316 M. J. 





LIGHTHILL 


Thus, on the linearized theory, the shock is reflected as a positive 
logarithmic infinity of pressure. Of course this cannot be exactly true, 
nor does the above constitute a proof that it be true since the conclusion 
contradicts the hypothesis of small disturbances. But in recent years 
aerodynamics has provided a large number of examples of a solution of a 
linearized problem taking infinite values locally. In many of these the 
quantitative conclusions were found to be sufficiently accurate away 


from the singularity, while qualitatively the values near the singularity 


were larger than elsewhere, though not infinite. Thus it may be permissible 


to conclude that a shock is reflected locally from the boundary layer as a 
pressure ridge, that is as a rapid compression followed by a rapid expan 
sion. The compression waves might be expected to form a shock-like 
envelope which would gradually encroach on the fan-like expansion and 
weaken it. 

The reflection of the discontinuity is better understood by observing 
the singularities of w(x,y) in the supersonic part y, < y <6 of the 
boundary layer. By (40) and (43) 


A M? iy FA — 
[l(k. y) oe I(y)| l (y) I | pik(Bs s a | iA | eik(pd s—S) A (46) 
; M, B 2ark 2ar| k 
Hence w(x,y) has a discontinuity of magnitude 
M(y)| M2(y)—1)-3 . 
y)| 7 | A ono :z pd+-S—s, (47) 
4 1? a 

: A M(y)| M?(y)—1 | ‘ 

while a(x, y)-4 — log(a+ps— S—s) (48) 
7 M, B 

is continuous at x pd+S-+s. 


The curve along which the singularity is propagated and reflected is 
shown plain in Fig. 1, where the sonic line and the edge of the boundary 
layer are shown as broken lines. The shock, after entering the layer at 
(— 88, 8), follows a characteristic, on which da/dy = {M?(y)—1 and so, 
by (30), x d+ S—s. It is reflected from its point of incidence on the 
sonic line (—fd+-S,y,) along the other characteristic through this point, 
x B5+-S-+-s. Finally it leaves the layer at (—fd+2S,8), and travels 
along the characteristic x»—By-+-2(85—S) = 0, as was already clear 
from (45). 

But on reflection from the sonic line it changes its character from a 
pressure ‘step’ to a pressure ‘ridge’. The magnitude of both singularities 


also increases, as the sonic line is approached, by a factor 


M(y)[ M*(y)—1}4/.M, B+. 
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sitive | This ‘focusing’ is not very intense. The author attempted a more exact 

| theory of the refraction of a shock entering a region of parallel and other- 
wise undisturbed steady flow at monotonically decreasing supersonic Mach 
number, which indicates that such an increase of strength occurs until a 
se the | the refraction theory does not apply, but the shock presumably weakens; 


his phenomenon does not appear on the present theory since no region 


n ofall state is reached in which the flow behind the shock is subsonic; after this 
away 


larity fsubsonic flow behind the shock is admitted.* 
issible 

\ 

Tas a ~\ 

Xpan- | * 

lk like 


n and %\ 
Pr'ving \ 7 ] 
. 1 \ 


S ecaner an mk 














a It a re ae Hr OS > 2 
Kia. 1 
18 The point of incidence on the sonic line (—fd+S,y,) is an algebraic 
nfinity of pressure in the theory, since both the step and the ridge become 
nfinitely tall, like |.47*(y)—1]|-!, on approaching it. The nature of the 
aie he nfinity, as approached along the sonic line can be deduced from equations 
iia 12) and (43). In fact 
ver at A 243 isen(x+B-S) sl Tae Nal 
- w(x, y,) ~ 5M; 16 M'(y, (A)! (—2)! |a+ps—S|-* (49) 
an ANY) 3): 
on the 
point sx approaches the singularity p5+-S. But in the subsonic layer, as 
ravels as been seen, the solutions are regular. 
cleat To conclude this section, consider an incident wave whose worst 
} singularity is a discontinuity in pressure gradient. Then no infinities of 
— pressure arise. The propagation of the singularity is predicted most easily 
arities by observing that @a/éa satisfies exactly the same differential equation 
* Though this deficiency of the theory is serious, it is difficult to believe it capable (as 
othesis (iii) of the introduction suggests) of producing the utter disagreement with 


experiment regarding distance of upstream influence. 
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and wall boundary condition as w. Hence if, for y > 8, @w/éx jumps by A 
at the line x+fy = 0, its resulting singularities in the layer and in the 
reflected wave will be exactly those given for w itself by equations (45), 
(47), (48), and (49) above. Those for wa may be obtained by integration, 

An incident wave consisting of a rapid but not discontinuous rise in 
pressure could be represented by a jump A in @a/éx followed closely by 
a jump —A in éw/éx. This would be reflected from the laminar boundary 
layer (according to the theory) as a local region of high positive éa/éx 
followed by a local region of high negative @w/éx—in other words by a 
ridge. This may be so even if the theory is inadequate when a itself has 
discontinuities. Again, a Prandtl-Meyer expansion, in which the pressure 
drops rapidly but not discontinuously, would be reflected as a pressure 
trough, i.e. as an expansion wave followed closely by a compression wave 
(with the latter probably forming a shock and bordering the expansion 
region). This prediction is to be tested in the Manchester University Fluid 
Motion Laboratory. 


7. Upstream influence of the disturbance 

Assume now that the front of the incident wave is the line x+fy = 0 
shown in Fig. 1. (Then f(~+ fy) = 0 for x+By < 0.) It will be investi- 
gated how far to the left of this line the disturbance in the boundary layer 
extends. More precisely the asymptotic behaviour of w(x, y) and g(a—By) 
as «-> —oo will be investigated, and the effective distance of upstream 
influence estimated principally from their logarithmic decrement. 

Since f(x) = 0 for « < 0, F(k) = (27)— | e-***f(ax) dx is regular in the 

0 

lower half-plane. Hence the only singularities of Il(k,y) or of G(k) in the 
lower half-plane, by (15) and (16), are poles at the zeros of their common 
denominator. Hence whenever x < 0 the functions w(x, y) and g(x), by 
(12), may be expressed as —2z7i times the sum of the residues of II(k, ye” 
and G(k)e** respectively at the zeros of their common denominator in the 
lower half-plane. 

Since in this section only the lower half of the k-plane will be considered, 
it is convenient to put k ix, and 


D(x) = Q,(—1t, 8)+BxQ(—ik, 8) (50) 


for the denominator of (15) and (16). D(x) is real for real « since Q involves 
k only through its square. If the residues of II(k,y) and G(k) at the nth 
zero x, of D(x) in the half-plane are R,,(y) and r,, respectively, then for 
z<0 : 

w(x, y) 271 > R,,(y)ek*, g(x) = —2m > PA a (51) 


n 
n=1 n=1 
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The zero x, with least real part gives the logarithmic decrement of w and 


as «> —oo; they decrease by a factor approximately e-“ per unit 
listance. 

Now as |x| 00, by (36) which is valid in the whole of the lower half 
plane, D(x) ~ 48! {M'(0)}-'«-1e*S sin(xa—}7). (52) 


Hence the zeros are asymptotically at the points (m+-})zo~! for integer m. 
These zeros for large |x| must be exactly real, since there can only be one 
ithe neighbourhood of each zero of (52), while complex zeros must occur 
pairs. 
It seems unlikely physically that any zero should be complex, as this 
would mean that the upstream disturbance was wavy. Also it is more 
kely that D(«) should be zero for real«x, where im {D(x«)} is automatically 
10, than for complex « where re{D(«)} and im{D(«)} must vanish simul- 
taneously. But this point will be returned to later. In the meantime only 
eal zeros of D(x) are considered, and a lower bound for them determined. 
It is probable that zeros of D(x) do exist near to the zeros 


5aqg-1, %7o-!, ana eae 


52). Since a good estimate of o is }7y, (its value, by (33), when M(y) 
sa linear function of y in the subsonic layer) these correspond to dis- 
turbances with logarithmic decrements of about 5yy! and more. These 
vould be effectively damped out in a single thickness y, of the subsonic 
yer. Clearly then the only hope of a result in agreement with the 
xperiments is from a zero of D(«) corresponding to the zero }7a~ of its 
symptotiec form (51), and preferably a good deal smaller still (say 
bout ;o7'). 

It might reasonably be expected that such a zero exists, if it were not 

it the zero at $zo~! of the asymptotic form of D(x) is one at which 
hat of D’(«) is 

D'(«) ~ —4B8210,{.M'(0)} ne tea < 0. (53) 
[his indicates that probably D’(«) < 0 at the zero of D(x) near x = }zo-!. 
nthe other hand, D(0) = Q,(0,5). But by the equation and boundary 
nditions for Q in the simple case k = 0, Q(y) = 3M?(y)/.M’?(0), so that 


D(0) = 3M?2/M’2(0) > 0. (54) 


his shows that there can only be an even number of zeros between 0 and 


the one near 37o0~!; and there may well be none. 


The possibility of small real zeros of D(x) will now be disproved. In fact 
twill be shown that if 

J 3 1 

« <|[[ ry — ay [ ae@) a:| 


0 0 
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then Q(—ix,y) is an increasing function for 0 < y < 8. Hence (50) cer 
tainly cannot vanish. 

To prove the result, suppose that there were values of y at which dQ/d) 
was not positive, and let y, be the least such. Then Q’(y) > 0 for 0 <y <y, 
but Q’(y,) 0. 

First assume ¥) < y,. Then by (55) 


Yo y 
1—x*? | {M-*(y)—1} dy | M*(z) dz > 0, (56 
0 0 
since the integrand is positive when y, < y < y,. Now by equation (13) for 
@, and its boundary condition (17), 
1Q/d j 
dQ/dy 3 rae 
2 | {M-2(2 L'O(y) dy. (57 
M%y) M’(0) - ‘ sl alli : 
0 


Hence for y < yy (when Q(y) > 0, M(y) < 1) 


TT] 
dQ 3 3 7 
a M*(y), Q : M?(z) dz. 5 
dy <0) YS Fpa@ | HO) sa 
0 
Hence by (56) 
Yo 
1—4$.M"*(O)q? | {M-2(y)— 1} Q(y) dy > 0, (59 
0 


so that by (57) Q’(y) 0, a contradiction. 
If alternatively y, y,, then by (57) the inequalities (58) hold wher 
¥Y < ¥, since in this region Q(y) > 0 and M(y) < 1. Hence by (55) 


Ix? | (M-*(y)—1} $20) Qly) dy > 0. (60 


But, since M(y) 1 and Q(y) > 0 when y, < y < yp, (60) implies (59 








and hence again (by (57)) the contradiction Q’(y)) > 0. This completes 


the proof that no zero of D(x) satisfies (55). 


The lower bound (55) for zeros of D(«)‘always exceeds yy1v2, since 
Ya y Yi i] 
| {M-*(y)—1} dy | M?*(z) dz - ar 2(y) | M?(z) dz) dy 
0 0 0 0 
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5) cer.| This smallest value y;'v2 is attained in the limit of a Mach number 
distribution which is small until just before y, and then suddenly shoots 
1 dQ/dy) up to |. For any ordinary distribution it would be close to yy!v12, its 
y value for a linear distribution of Mach number in the subsonic layer 
M(y) = y/y, for 0 << y < y,). But even a logarithmic decrement yy!v2 
for the upstream disturbance, corresponding to a fourfold decrease per 
nit thickness y, of the subsonic layer, is in complete disagreement with 

_.| experiment. 

{pf 


| Actually it has been shown that under condition (55) Q,, is essentially 
positive, so that a zero of D(«), by (50), is still remote. This seems to the 


(13)1 wthor overwhelming evidence that the first zero is that near 370, 
| though, as the last paragraph shows, the qualitative result is independent 


of this conjecture. 


= However, an argument which supports the conjecture that the only 


; | eros. even complex ones, are at points near 
a-!, tra-!, Baro-'.,... (62) 
s based on following the quantity 

(58 A(k) arg] Q,(k, 8) ‘ BikQ(k, 8)| (63) 
round a large semicircle in the lower half-plane with centre the origin. 
[his quantity is known for large / in the lower half-plane by (36), and for 
(59)! k= 0 it is zero or some multiple of 27 by (54). Also for small |k| the 
solution obtained by solving (57) by successive approximations would 
ppear to be fairly smooth. It might be assumed, therefore, that the 
| when) snown value of A(k) at k = 0 is joined on to those for large positive and 

iia irge negative / somewhat smoothly. 
) Now it may be shown, by equating the change in (63) on following the 
semicircular path, from its lowest point in the positive sense back to 
(00) the lowest point, to 27 times the number of encircled zeros of D(«), that the 
njecture that there is one zero of D(x) corresponding to each of (62), and 
es (59 that there are no others, is equivalent to a behaviour of A(k) near k 0 
npletes asin the broken line in Fig. 2. (Here the plain lines indicate the asymp- 
totes for large |.) But if an additional zero in the lower half-plane exists, 
nee the behaviour would have to be somewhat like that of the dotted line. 


, And if D(«) has complex zeros (which must occur in pairs), a behaviour 
ike that of the dash-dotted line would be necessary. But actually, since 
Q,(0,5) and Q(0,8) are positive, A(/) must increase with k at k = 0, by 
| (63). Hence the dash-dotted line must be replaced by a curve with the 
(g])| ‘me asymptotes cutting the axis at the same place but with positive 
slope. The mathematical unlikelihood of this, combined with the physical 

} 5092.11 Y 
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unlikelihood of the upstream disturbance being wavy, convinces the 
present author that the only zeros are in one-to-one correspondence with 
those in (62) and approach them more and 


P™ more closely as K —> OO. 
, Finally the actual term in the series of 
Ak) 7 exponentials for the reflected wave g(2— By) 
_ fox corresponding to the first zero of D(x), is 
c. cS; estimated from the asymptotic formulae, in 





order to quash the possible objection that 
though the exponential decreases very rapid- 
ly its coefficient may be uncommonly vast. 
Only the reflected wave need be considered, 
since it is this, not changes in the boundary 
layer, which is actually observed (by the 
Schlieren process). It should also be noticed 
that further terms need not be considered 
for no large rate of rise of the coefficients in 
the series (51) is permissible since each must 
converge whenever x < 0. 

When the denominator of G(k) is first 


zero, very roughly x = }zo~! and its deriva- 





tive is given by (53). The numerator of 
G(k) is not given by (37), which in the lower 





half-plane neglects a larger and predominant 
term O(k-2e'*S+ke), At the zero this is 
O(«-*e**). Hence by (16) the residue of G(k) at this zero is only a 
moderate multiple of 


Fic. 2. 


f 
. 


A . din — . ; ee J, : 
p57 B8/20 Ff ( é p57 P8/20 | o—5rz 10 f (x) dx eo Bd 20()( ). (64) 


Hence, by (51), we have roughly 


OTT 


g(x) ~ exp (x-+ 288) 
Oo 


O(f). (65) 








Thus reflection along lines upstream of the reflected wave front 
x—By-+-2(pd5—S) = 0 

(which is illustrated in Fig. 1) by more than about 2S-+o would be 

negligible. Now 2S-+o is about two boundary layer thicknesses. In any 

case the 2S portion already corresponds to signals sent back by the 

incident wave on hitting the outside of the layer. 

The conclusion that upstream influence of a weak disturbance on a 
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‘boundary layer of inviscid fluid’ is practically negligible, as stated and 
discussed in the introduction, is therefore inescapable. 


§, Examples. Connexion with Tsien and Finston’s theory 

The author first carried out the work of sections 2-7 for the case of a 
boundary layer with the Mach number varying linearly across it. For this 
case M(y) = y/y, (0 < y < 8), and 6 = y, M, and Q(k, y) is found to be 


the confluent hypergeometric function 
O(k, y) ye NBG —jtky,; 3; tky?/y,). (66) 


However, he later found that just as much could be discovered in the 
general case. 

More interesting is an attempt to achieve the Tsien and Finston theory 
1) as a limiting case of the present theory. In their theory M(y) is a 
onstant .V/, | in the whole layer 0 < y <6. This value does not 
satisfy the conditions of continuity down to M(0) 0 and up to M(d) = M, 
mposed in this paper (section 2). But it is a limiting case, as e > 0, of the 


istribution 


M(y u,” (O<y <e), M, (e<y< 85—e), | 
: (67) 
Ve —(M,—M,)°—4  (8—e <y <8) 
which for « 0 does satisfy the conditions. 
Now in 0 y e, for this distribution, by (57), 
dQ 3 M2) ) . se * , 
. 14+ O(e2)) = By2f1+ O(e2)}. (68) 
dy M’?(0) 
Hence at y € 
(e) e? + O(e°), Q'(e) 3e2+ O(e4). (69) 
But in « y d—e, by (13) 

QO" (y)—k*B5 Oly) = 9, (70) 
vhere B, (1—M3)'. Hence Q(y) is a combination of sinh(kf, y) and 
sh(k8, y) in this interval, which by (69) must be 

Vy) 3e7(K8,)-! sinh(kB, y)+ O(e*). (71) 
Hence. at y oO e. 
9 2/ ].4 hint QQ S\1 3 
Q(d—e) = 3e*(kB,)~! sinh(kB, 6)+ O(e*), (72) 


Q'(d—e) 3e* cosh(kB, 6) + O(e*). 
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But in d—e < y <4, by (13), on the 
_ < re | ance 
Q'(y Q'(5—e) , 15 , Q'(5—e) . 
|. 2 {M- 1\¢ 1 O(e8 73\) has al 
NEC ‘ Ne i (y)—1;O(y) dy NE (e*).  (73)) has al 
o—eE 

Hence 
Q'(8) = (Mi/M})Q' (6—e) + O(e?) = 3(M?/M3)e? cosh(kB, 5)+ O(e), | | : “e 
; 2 (74 2, Tee 
(0) Q(d—e)+ O(e*) 3e* *(kBo) 1 sinh(kB, 5) | O(e?), | 3.9.2 
Thus, by its definition (50), | 4 A. 
D(k) 3(.M?/M3)e? cos(KB, 5) + Bxde2(«B,)-! sin(«B, 8)+- O(e), vad G 
. B 5. P. 
and D() 0 in the limit as « — 0 when y 
M? B, cos(«By5)+ M3 B sin(«B,5) = 0. (76)| 6 1. 


Tsien and Finston introduce a parameter @ which in the notation of this 
paper satisfies 


M?B, 
tan 36 raat (77 
ai MzB 
Then (76) becomes sin(4@-+-«f, 8) 0, giving for the first zero of D(x 
Ky (ar +6) dB,. (78 


Since, by (33) in the limit as « > 0, 58, is o, this root, in lying between 
to! and zo~!, agrees with the general theory of section 7. (The lower 
bound (55) has become 2!o-1.) 

But the logarithmic decrement of the disturbance in Tsien and Finston’s 
theory (1) is (4%—46)/58,, which can be quite small, by (77), if M? B,/ M35 
is large. 

The cause of this is essentially that they take Q(y) proportional t 
cosh(kf, y), not to sinh(kf, y), in their layer 0 < y < 6 of uniform subsoni 
Mach number J/,. This is derived from the condition that v = 0 at the 
wall y 0, which (since U(y) ~ 0 there) gives that the deflexion 7 = 0 
and hence by (6) éw/éy = 0. 

3ut as a result of the present limiting process, by (12), (15), (69), and | 
(75), the boundary conditions at the side of the subsonic layer nearest to 
the wall are of the form w = 0, éw/éy ~ 0. Thus n + 0, but the pressure p 
takes a constant value P. Interpreted physically this means that the very 
thin vortex layer separating the main subsonic layer from the wall is not 
of uniform thickness as a result of the disturbance, but is at uniform | 
pressure. 

[t seems pointless to argue about the relative reasonableness of this 
condition or of Tsien and Finston’s, since the model is in any case artificial 
The author’s purpose has merely been to show why, though the Tsien and 
Finston theory is a limiting case of the present theory, the results it gives 
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n the distance of upstream influence are so very different. The signifi- 
ance of the divergence to the theory of shock boundary layer interaction 


has already been fully discussed in section 1. 
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ON SOME GENERAL THEOREMS IN THE LINEARIZED 
THEORY OF COMPRESSIBLE FLOW 


By F. URSELL and G. N. WARD 
(Department of Mathematics, The University, Manchester) 


[Received 6 December 1949] 


SUMMARY 

The linearized equations of motion for steady irrotational isentropic flow past ; 
thin body with a sharp trailing edge are used to evaluate the leading terms in th 
components of the force acting on the body when there is a vortex wake. It is show: 
that with certain additional assumptions required by the appearance of infinities at 
the edges of the wake the linearized theory is adequate for the calculation of the lift 
and drag on a thin body, but not for the calculation of the lateral force, except in 
special cases, e.g. when the wing is nearly plane. 

A quadratic identity involving the drag is used to show that there is at most on« 
solution with given boundary conditions. The corresponding bilinear identity is used 
to obtain relations between the flow past a thin body and the reversed flow past th 
same body, when the roles of leading and trailing edges are reversed. These relations 
hold for both subsonic and supersonic flow and include all flow-reversal theorems 
obtained by earlier writers. 


1. Introduction 

GREAT progress has been made in recent years in the theory of the steady 
two-dimensional flow past thin wings with a sharp trailing edge. For 
subsonic motions it is found experimentally that the action of viscosity 


establishes a flow which is smooth near the trailing edge where a thin | 


wake leaves the body. When the flow has become steady, the velocity 
gradients are not large anywhere, so that ultimately the action of viscosity 
is small. This situation is idealized by the assumptions that the viscosity 
vanishes, that the wake is infinitely thin, and that the velocity at the 
sharp trailing edge is not infinite. These conditions, the last of which is 
due to Kutta and Joukowski, determine uniquely the flow past the wing. 
They define’a cyclic irrotational motion. The velocity is continuous every- 
where, i.e. there is no wake; but the velocity potential is many-valued 
since there is a circulation about the wing. For supersonic flow there is not, 
in general, a velocity potential since particles can acquire vorticity by 
passing through a shock, but when the wing is thin the shock is weak and 
the vorticity can and will be neglected, so that there is no wake even in 
this case. 

In many cases the incidence, thickness ratio, and camber of a wing are 
small. A linearized two-dimensional theory has been constructed in which 
the exact equations and boundary conditions are replaced by linear 
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pproximate equations and boundary conditions which, it is hoped, will 
sive the leading terms in the dynamical variables as incidence, thickness 
ratio, and camber tend to zero. Where comparison with the exact theory 
is possible, it has been found that this is the case except near the leading 
edge where linearized theory gives an infinite velocity. In spite of this 
liserepaney the linearized theory gives the correct leading term for the 
forces acting on the wing, for these are given by the momentum integral 
taken over any curve surrounding the wing, and at a distance from the 
wing the linearized theory is applicable. 

The theory of the flow past a three-dimensional body with a sharp 
trailing edge is in a less satisfactory state. It is again assumed that 
viscosity enters only through the condition of Kutta and Joukowski, that 
the effect of shocks can be neglected, and that the wake leaving the sharp 
trailing edge is infinitely thin. However, it can no longer be deduced that 
there is no wake; instead it is assumed that in the ideal flow a vortex sheet 
aves the sharp trailing edge. At a distance from the body this sheet is 
supposed to roll up more and more near the edges until it has ultimately 
rolled up into two vortex cores of finite diameter. No exact solution of the 
non-linear equations has yet been found. 

When the incidence, thickness ratio, and camber are small it is reason- 
ble to attempt a linearization of the equations. The linearized equations 
f motion obtained in this way have been solved for linearized boundary 
conditions in the case of the nearly plane circular wing in incompressible 
subsonic flow (1, 2) and for a wide class of wings in supersonic flow. It is 
found that infinite velocities occur at the leading edge and at the edges 
fthe vortex sheet. 

In the present paper an attempt is made to evaluate the leading terms 
nthe components of the force acting on the body when there is a vortex 
vake. It will be shown that with certain additional assumptions (required 
by the appearance of infinities at the edges of the wake) the linearized 
theory is adequate for the calculation of the lift and drag on a thin body, 
but only so for the calculation of the lateral force in special cases, e.g. for 
inearly plane wing. For these calculations it has been found convenient 
to derive and discuss afresh the linearized equations of motion and the 
boundary conditions. The latter cannot be linearized unless the local 
thickness ratio of the body is small in one direction only. The term ‘body’ 
shere used to include a system of isolated bodies. 

It is easily seen that the leading term for the drag force is non-negative; 
this property is used to show that the linearized boundary conditions, 
including the Kutta—Joukowski condition, make the solution (if any) 
unique. From the drag formula an integral identity is derived which may 
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be used to connect a flow past a given body with the reversed flow in which 
the velocity at infinity is equal in magnitude and opposite in sign so that 
the role of leading and trailing edges is reversed. Several special cases are 
known for supersonic flow (3, 4, 5, 6, 7, 8, 9). The following simple case 
is typical: for any thin, nearly plane wing in a supersonic flow, at the 
trailing edge of which the Kutta—Joukowski condition is satisfied. a 
reversal of the stream direction leaves the initial slope of the curve of lift 
against incidence unaltered provided that the Kutta—Joukowski condition 
is also satisfied at the trailing edge for the new flow. This theorem is here 
proved under wide conditions and for subsonic flow. Examples can be 
multiplied indefinitely, a typical more complicated one being that, for 
any thin body, the rolling moment due to a given rate of yaw is equal to 
the yawing moment in the reversed flow due to the same rate of roll, the 
Kutta—Joukowski condition being satisfied at the trailing edges in both 
flows. 

The importance of these theorems lies in the fact that (in supersonic 
flow at any rate) it is frequently easier to compute the reversed flow 
rather than the original flow. For instance, to find the lift on a flat delta 
wing with supersonic leading and trailing edges in supersonic flow, it would 
be convenient to consider the lift in the reversed flow in which the wing 
advances with its base forward. Some of the theorems mentioned in this 
paper have also been given in a recent paper by C. E. Brown (9) who 
generalizes a method due to M. M. Munk (6) which involves physical 
considerations substantially equivalent to the method of the present paper. 

The infinite vortex strength at the edges of the wake requires special 
mention. Its physical meaning is not yet clear. For any finite incidence 
the strength at the edges in the exact flow must clearly vanish since 
otherwise there would be finite suction forces acting on the edges. It may 
be that the maximum vortex strength tends to zero more slowly than the 
incidence, and that it is located nearer and nearer to the edge, as incidence, 
thickness ratio, and camber are made to tend to zero at the same rate. 
The removal of this difficulty requires further study; it would then be 
possible to justify or refute the additional assumption which is now 
required for the calculation of the force. The calculation of the lateral 
force also requires further elucidation. This force is of the second order, 
and the present paper suggests that the first approximation is inadequate 
for its calculation except in special cases, e.g. when the body is nearly 
plane. It is not known how far this apparent inadequacy represents the 
true state of affairs. 

In all this work, much of the detail required for a rigorous mathematical 
treatment is omitted from considerations of space, but all essential steps 
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FLOW 329 
in the arguments are given. For the same reason the final stages in the 
proofs of the uniqueness theorem are only indicated and appropriate 
references are given. Vector methods are used throughout. 


2. The equation of motion 
Let p, p, and u be respectively the pressure, density, and particle 
velocity at a point in a steady flow of an inviscid gas, then the equations 
fmotion, expressing the conservation of mass and momentum, are 
V .(pu) 0 (1) 
u.Vu-+(1/p)Vp (0). (2) 


‘hese equations are too general for the present purpose, and it is assumed 
that the motion is irrotational and isentropic. The condition of irrota- 
tionality is expressed by 

VAu 0, (3) 
s that u is the gradient of a scalar point function, the velocity potential. 
Since the flow is isentropic there is a functional relationship between the 
pressure and the density, and with the aid of (3), (2) can be integrated to 


give Bernoulli’s equation 


fu? = constant, (4) 


| dp/p+3 
which is an expression of the conservation of energy. Also the local 

“oustic velocity, a, at any point is given by 

a* dp dp. (5) 

The velocity u at any point is now expressed in the form 

u U-+yv, (6) 
where U is the velocity a long way ahead of the body, and it is assumed 
that squares and higher powers and products of v/U and its derivatives 
in be neglected in the equations of motion. Thus the departure of the 
stream from uniform constant conditions is expressed as a small perturba- 
tion, and the linearized equations for steady irrotational flow of a com- 
pressible fluid are obtained. 

Let Py, py, and U be respectively the pressure, density, and velocity of 
the undisturbed stream, then, for sufficiently small departures from 
inform flow, 

p p 


(p y) .. ae 
| P | \ Po ' ' dp 2} 0 , 


ay 


— py 1f[p—po\"(dp\ . a 
p—} f Pe) (2) i (7) 
Po - Po dp 0 
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From (4), (5), and (7), they 
I—p 1 [p—p,y\" Hi ees 
cS... (! Ps) +... = fU2—(U-+vy 
Po =49\ Po 
—U.v—}v? (8 
and therefore 
P—Po_ _y.y—jy24 (U.v)?+... wher 
Po s 2ap suffis 
U.v—iv.w- (9) | (23): 
where w = v—UU.v/a =@.y, say, (10 “ 
, 
® being the self-conjugate dyadic (cf. 10) I—UU/a?. The vector w is pa 
proportional to the perturbation mass flux vector, as is shown later. The 
density p is expanded: 
. 
P = pot(p—Ppo)(dp/dp) y+... (11) | 
which by virtue of (9) becomes 3 T 
p Po( 1 U.v ar t cae (12 In 
‘ hence . - . ail 
und hen pu = p,(l1—U.v/a?+...)(U+v) ‘ : 
TTT 9 to tn 
p(U+v—UU.v/a2-...) " 
an 
J-Lw)_ 3) 
p,y(U+w) (13) ee 
Thus, on linearization, equation (1) becomes alle 
V.w (). (14) For | 
Equation (13) gives the physical interpretation of w, since pu is the mass __ plan 
flux vector and, to the first order, bodi 
U+w (p/P) )U. (15) “a 
. ; ; th ; nw 
The vector U-+w is the mass flux vector divided by the undisturbed 
. . ; a : an 
stream density, so that w is the perturbation flux vector. The introduction " 
’ 
of the vector w is due to A. Robinson (11), who called it the reduced : 
. = ; ‘ ‘ ; : some 
current velocity. Thus the linearized equations for steady irrotational | ' 
, ‘ ; ‘ Se ekrs : actu 
isentropic flow of a compressible inviscid fluid are Tl 
ie 
u U+v, (16) hody 
VAv = 0, (17) | unit 
7.9 = @, (18) | othe 


w = @.v = v—UU.v/az. 


Equations (16) to (19) are independent of any particular coordinate 


(19) | tion: 
(j 


frame of reference, but, by introducing a right-handed system of ortho- 


gonal Cartesian coordinates x, y, 
undisturbed stream, and a potential 4 defined by 


v = Vd, 


z with x-axis in the direction of the | ii 


(20) 
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they can be written in the form 
u = (U+¢,, ¢,. 9%); (21) 
w = {(1—M")4,, ¢,, ¢.}, (22) 
" (1 M*)¢$,., I Pyy T $.. 0, (23) 
where M U'/a, is the Mach number of the undisturbed stream, and 
suffixes denote partial derivatives. If the motion be subsonic (M < 1), 
(9)| (23) or (17) to (19) are elliptic, and if the motion be supersonic (M > 1), 
(jo) | the equations are hyperbolic and real characteristic surfaces exist. 
| The self-conjugate dyadic ® can be expressed in matrix form for this 
ee ordinate system as 
r. The 1— m2 0. (0 
Q, 2 0 
(11 0, 0, l 
3, The boundary conditions 
(12 In formulating boundary conditions for linearized theory it is necessary 
to distinguish between bodies whose local thickness in any direction normal 
the undisturbed stream is small (compared with the length of the body), 
7 nd bodies whose local thickness normal to the undisturbed stream is small 
(13) | in one direction only at any given point. Bodies of the first kind are here 
lled slender bodies and bodies of the second kind are called thin bodies. 
(14) | Forexample a pointed body of revolution is a slender body, while nearly 
‘mass plane wings and quasi-cylindrical bodies are thin bodies. There are some 
dies, such as a pointed body of revolution with attached thin wings, 
(15) Which are slender in some part (the body) and thin in another (the wings). 
ileal Inwhat follows only thin bodies of finite overall dimensions are considered, 
earre nd the following boundary conditions apply only to such bodies. 


neal By definition the boundaries of thin bodies lie at a small distance from 


some surface whose generators are parallel to U, and the projection of the 


ee , actual body surface on to this surface is called the mean body surface. 
lhe boundary conditions on the body surface are applied on this mean 
(16 ge Or ; 
m ody surface. If n* be the unit normal to the body surface and n be the 
(17 init normal to the mean body surface, both being inward normals; or, for 
(18 ther surfaces, if n be the appropriate normal, then the boundary condi- 
(19) | tions to be satisfied by the perturbation velocity are 
dinate i) on the mean body surface 
ortho n.v n*.U (n.U = 0), 
of the ll) on a mean surface corresponding to a free boundary outside which 
the gas is at rest with given pressure, 
(20) | 


.¥ known constant (nn. U = 0), 
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(iii 


at sharp trailing edges of the body Vv is finite, 


(iv) in the wake of the body, any vortex sheets extending downstream 


— 


from the sharp trailing edges are surfaces whose generators are 
parallel to U, with edges also parallel to U, at which n.v and U.y 
are continuous, 

(v) in supersonic flow, on characteristic surfaces at which v is dis- 
continuous, nA Vv and n.v are continuous, 


(vi) at all other points in the fluid, v is continuous, 

(vii) in subsonic flow, if r be distance from the nearest point of the body 
or the vortex sheet (if any), r?v is bounded as r — 0, and, in super- 
sonic flow, V = 0 on any plane normal to U and upstream from the 
body, 


(vili) in subsonic flow, if R be distance from the nearest point of the body, 
R°?U.v is bounded as R > ~, 

(ix) the over-all aerodynamic force on the body is finite (see section 4). 

These boundary conditions are not all independent, but it is convenient 
to have a complete list for later reference. 

The justification of the conditions (i) to (ix) is as follows. 

(i) Since the fluid under consideration is inviscid, it can slip against the 
boundaries, and the boundary condition at the body surface is 

n*.(U-+-v) 0, (24) 

which is the mathematical expression of the assumption that the fluid 
does not separate from the body. Since the body is thin and the body 


turbed stream, n*.U must be small compared with U, and the difference 
between n*.v on the body and n.v on the mean surface is a small quantity 
of the second order which may be neglected to give the linearized boundary 
condition (i) from (24). The solution of the linearized equations of motion 
with this boundary condition can lead to infinite values of v at sharp edges, 
and the condition v < U is clearly violated. This difficulty is discussed in 
section 4. This complete linearization of the boundary condition at the 
body surface is not possible for slender bodies. 

(ii) On free boundaries where the gas is at rest on one side, the condition 
of constant pressure may be expressed by the linearized approximation to 
Bernoulli’s equation in the form 

U.v = known constant, (25) 
and, if the free boundary extend upstream to infinity, the constant is zero. 
This boundary condition can be applied on the mean surface of the free 
boundary whose generators are parallel to U, to give a completely linearized 
boundary condition as before without causing loss of accuracy. 
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iii) At sharp trailing edges it is necessary to apply the condition that 

must be bounded there in accordance with the Kutta—Joukowski 
hypothesis in order that the mathematical solution shall agree with 
physical experience. It is found that with this condition the velocity at 
harp leading edges becomes infinite in general, as has been mentioned 
inder (i) above 

iv) A result of applying this Kutta—Joukowski condition is that sur- 
faces on Which a tangential component of velocity is discontinuous appear 
n the mathematical solution. Such surfaces of discontinuity, or vortex 
sheets. are not incompatible with the assumption that the fluid is inviscid, 


it certain physical conditions must be satisfied on them. These con- 





| ditions are that the velocity must be purely tangential, and that the 
ressure must be continuous through the vortex sheet. If v,, v, be the 
| perturbation velocities on either side of the sheet, and n be the normal to 


the sheet. then these conditions are 


n.(U-+v,) n.(U-+-v.,) 0, (26) 


P 
nd, if the perturbation velocities be small compared with U 


U.v,+4v,.w, U.v,+4Vv,.W. (27) 


rom (10). Equation (26) shows that U.n must be small compared with U, 
that n is approximately perpendicular to U. Equation (27) may be 


ritten 


(v,—v,).{U+3(w,+w,)} = 0, (28) 


ifw be the surface vorticity 
n/ w.{U+3(w,+w,)} = 0. (29) 


thus, either w — 0 and there is no vortex sheet at all, or n,w, and U are 
pproximately coplanar. Since n and w are orthogonal, and n and U are 
pproximately so, must be approximately parallel to U. It is usual to 
lopt Prandtl’s assumption that the vorticity w is exactly parallel to U at 
| points on the vortex sheet, so that the vortex sheet is a surface with 
generators and edges parallel to U and which extends from the trailing 
edge causing it to go to infinity downstream. On this assumption, only the 

hear approximation to (27), namely 
U.v, = U.yv,, 

8 satisfied on this surface in general. 
The boundary conditions on the new surface, the linearized vortex 
sheet. are os 
¥,.@ = ¥..m, (26°) 


equation (27’), and w constant on lines parallel to U. 
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[t is found that v becomes infinite in general at any edges of the vortex 
sheet on these assumptions, so that the condition (27’) is not valid in this 
form near to and at the edges, and the above assumptions become difficult 
to justify. In fact, the wake cannot extend to infinity as a vortex sheet 
as has been assumed, since (26) and (26’) show that the streamlines do not 
lie in the assumed vortex sheet even on linearized theory. In more exact 
theory it is assumed that a steady motion exists in which the vortex sheet 
rolls up behind the body to form vortex cores. It seems reasonable to 
assume further that the rate of rolling up depends upon the vortex strength: 
at a fixed distance behind the body the vortex wake will ultimately have 
the linearized form as the strength tends to zero, while on the other hand, 
for a given strength, a distance can always be found at which the vortex 
sheet has rolled up to a given extent. On the present assumptions the 
linearized theory represents an approximation valid in a region which 
increases as the strength decreases, but not valid uniformly in the whole 
space. 

(v) In supersonic flow real characteristic surfaces exist on which v may 
be discontinuous. Consideration of the conservation of mass and momen- 
tum through such a surface leads to the conclusions that the normal 
component of w and the tangential component of Vv must be continuous 
there. Conservation of energy shows that the constant in Bernoulli’s 
equation (4) is continuous so that the expansion (10) is valid on both sides 
of such a surface. 

(vi) Discontinuities of v can occur only on vortex sheets and real 
characteristic surfaces, both of which have been considered. 

(vii) From Kelvin’s circulation theorem, the total vorticity in all the 
vortex sheets must be zero so that, at large distances from the body and 
the vortex wake, the flow must approximate to that of a semi-infinite 
vortex pair. The given conditions for subsonic flow follow from this if the 
vortex pair be of finite strength. 

In supersonic flow, disturbances cannot be propagated upstream from 
the body so that the flow there must be that of the undisturbed stream 
and v = 0. In fact, since v and its derivatives can be discontinuous only 
on characteristic surfaces, v 0 everywhere upstream from the envelope 
of downstream characteristic semi-cones with vertices on the leading edges 
of the body. 

(viii) This condition follows from the condition (iv) since the vortex lines 
are of constant strength and parallel to U. 

(ix) This condition is clearly necessary if the solution is to be physically 
significant. Its consequences are discussed in section 4. For subsonic flow 


it is implicitly contained in conditions (vii) and (viii). 














4, Th 
Ad 


Jinear 
infinit 
Berne 
priate 

Fo 


tum t 


Here 
S con 
sides 
is fini 


since 
the | 
vorte 


Henc 


n be 
chan 
veloc 
ont 
by (; 
that 

If 


com) 


the 
Ther 


appl 


If y 


con 





Vortex 
in this 
ifticult 
sheet 
do not 
exact 
- sheet 
ble to 
noth 
have 
hand, 
vortex 
is the 
which 


W hole 


\ may 
omen 
ormal 
nuous 
oulli’s 


1 sides 
| real 


ll the 
y and 
ifinite 
if the 


from 
tream 
s only 
“elope 


edges 
x lines 


‘ically 


c¢ fle WwW 








LINEARIZED THEORY OF COMPRESSIBLE FLOW 





4, The aerodynamic forces 
A difficulty arises in the calculation of the aerodynamic forces on the 
jnearized theory at points where the perturbation velocity becomes 
infinite, since the usual approximation (either linear or quadratic) to 
Bernoulli’s equation is certainly not valid at such points, and the appro- 
priate form for the pressure has to be determined from other considerations, 
For steady inviscid flow with no impressed force field, Euler’s momen- 


tum theorem for any closed surface inside which wu is not singular is 
| (pn-+puu.n)dS = 0. (30) 


Here the closed surface is taken to consist of the body surface &*, a surface 
‘completely enclosing the body, and a surface o* which lies close to both 
sides of any vortex sheet in the flow between &* and S. If the velocity 


sfinite on the body, the contribution from =* is 





F | pndsS (31) 
since u.n* — 0 on the body surface; thus F is the aerodynamic force on 
the body. On the surface o*, since n.u and p are continuous through a 

rtex sheet, the integral vanishes when taken over both sides of the sheet. 
Hence : 
F | (pn-+puu.n) dS, (32) 
n being the outward unit normal to S. This integral is invariant for 
hanges of S. When there are sharp edges on the body, at which the 
elocity becomes infinite, there are finite suction forces on such edges. 
(ontinuity considerations indicate that the total force F is to be defined 
32) rather than by (31). In this paper it will therefore be supposed 
hat the aerodynamic force is given by (32). 
Ifa surface S exists on which the perturbation velocity is uniformly s’aall 
mpared with U, as is assumed in references (3) and (12), the integrand 
32) may be approximated by using the expansion (9) for the pressure, 
the expansion (12) for the density, and neglecting third-order terms. 
then, on taking account of mass conservation inside S, (32) becomes 


pproximately 


F = p, | {(U.v+4v.w)n—|(1—U.v/a?)U.n+v.n|v} dS 
Po ( ((U.v+iv.w)n—(U+w).nv} dS. (33) 


iv and w satisfy the linearized equations of flow and the boundary 
nditions of section 3, this is the momentum integral for linearized theory 
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as given in (4) and (11), and it is shown there that this integral is invarian} 
for changes of S such that S always encloses the same singularities, sine 
the divergence of the integrand vanishes in linearized theory. However 
in a flow with a vortex wake satisfying the boundary condition (iv) of 
section 3, the integral in (33) is not invariant for all changes of S such that 
S encloses the body, as is assumed in (3) and (12), since integrals of the 
form (33) taken over two surfaces S and S’ differ by the value of a similar 
integral taken over that part of the linearized vortex sheet which lies 
between S and S’. This last integral does not vanish in general, because 
the velocity at the edges of the vortex sheet is infinite on the linearized 
theory, and there is a contribution which depends upon the quadrati 


> 


terms in (33). Thus the integral (33) is not a valid approximation to the 
integral (32), and does not give the forces correctly in linearized theory 


These difficulties which arise in the application of the momentum method 





to the calculation of the aerodynamic forces are mostly due to the assump- 
















tion of the idealized form taken for the vortex sheet in linearized theory 
and it is necessary to consider the approximation to the integral (32) ir 
more detail. 
Let v* and w* be the exact perturbation velocity and the exact per- 
turbation flux vector respectively for flow past a thin body, so that 
u U-+v*, pu p(U 1 wy*), (34 
Then (32) can be written 
F | {[p9 U. v*—(p—p yt pyU.v*)|n—p,(U+w*).nv*! dS. (35 
8 
[f ¢* is the exact perturbation potential for subsonic flow or is the per 
turbation potential which gives v* to the second order in supersonic flow 
then, either exactly in the former case, or to the second order in the latter 


case, 





(U.v*n—U.nv*) dS |} UA (nA v*) dS 
S RS 
UA | nAVd*dS UA | d*ds, (36 
where the line contour y* is both sides of the curve in which the vortex 
wake intersects S for the exact flow. Equation (35) then becomes 
F py U | d* ds | {(p—po+p, U.v*)n- po V*w* ni dS. (37 


S 


[t is not difficult to show that (33) can be written in a similar form, namely, 


F = p,U/ [ dds—p, { (vw.n—3v.wn) dS, (38 
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THEORY OF 


where ¢ is the linearized perturbation potential and y is both sides of the 


urve in which the linearized vortex sheet intersects NS. 


The linearized theory tries to represent the limiting flow when the 


thickness, camber, twist, incidence, etc., of the body all tend to zero at 


proportional rates, but the approach to the limit is not uniform all over 


the fluid, as has been pointed out already in the discussion of the vortex 


sake in section 3. The perturbation velocities become small compared 
th U except at leading subsonic edges which become sharp in the limit, 
nd at the edges of the vortex sheet. The perturbation potential becomes 
small uniformly over the whole field in the limit. At points on S which 
re not close to the edges of the vortex wake, the integrand of the surface 
ntegral in (37) tends to the value of the integrand of the corresponding 
surface integral in (38), as is seen by using the expansion (9) for the 

pressure: that is 
py V*w* .n+(p—po+p, U.v*)n > p,(vw.n—}v.wn), (39) 


0 
t this is not so in neighbourhoods of the edges of the vortex wake. If it 
n be assumed that, in the limit, the contributions of the two expressions 


39) are equal (and it is not known whether this assumption is justified), 


py V*¥w* .n-+-(p— pot py U.v*)n} dS > py | (Vw.n—}v.wn)dS, (40) 


the difference between the two integrals being of the third order. 
The difference between the two line integrals in (37) and (35) is of the 


econd order (since the difference between the line elements on y* 


and y 
sof the first order), and this difference is thus of the same order as the 
irface integrals 
In general the integral over y* is difficult to evaluate except when 
[*, where ['* is both sides of the trailing edge of the body. The 
surface S must then be a surface S* which passes through the trailing edge 


{the body. Finally, then, the force on the body to the second order is 


F = p, U | dds o, | (Vw.n bv.wn) dS+p,U / ( (d*—¢) ds, (41) 


rO 


I Se nL 


F = p, | {(U.v+4v.w)n—(U+w).nv}dS+p,UA | (¢*—¢)ds. (42) 
ed a 

[he surface integral in (42) is that in (33) for the special case when S S*., 

‘ince the divergences of the integrands of the surface integrals in (41) and 

12) vanish, the surface S* can be transformed into any surface which 


ncloses the body and passes through its trailing edge: in particular S* 


in be the actual body surface =* in (41) and (42), or the mean body 
092.11 


Z 
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surface & plus a surface enclosing the sharp leading edge (where thy 
assumptions of linearized theory are certainly not valid) in (41). Thus, o1 
neglecting third-order quantities, 


F = —p, | (U.v+4v.w)n* dS+p,UA [| (6*—¢) ds 
=* T° 
‘(p—p,n* dS+p,UA { (¢*—¢) ds, (43 





so rs 


where (p—q) is the pressure given by the quadratic approximation (10) t 
Bernoulli’s equation, or 


F po US 


$a | (vw.n—jv.wn) dS-+ p, U 


1% : r 


“(b*—d) ds. (44 








In general the calculation of the force correct to the second order 
requires a knowledge of the discontinuities of ¢* correct to the second 
order. In equations (41) to (44) for the force, the term involving ¢* has 
been separated so that it occurs as a second-order term at most, the first- 
order forces being given by the other terms in these equations. For some 
special bodies, the component of the term involving ¢*—d¢ which is 
perpendicular to the first-order force can be shown to be of the third 
order, so that the lateral force is obtained correctly to the second order 
from the linearized solution alone. Such a case is that of a nearly plane 
wing, for which the component of ds on I'* normal to the mean plane is 
of the first order, so that the lateral force parallel to the mean plane is 
given correctly to the second order by linearized theory. 

In all cases the drag force parallel to U does not involve the line integrals 
and is always a second-order force which is given from the linearized solu- 
tion alone. 

When evaluating the second-order forces, care must be taken to include 
the force at any sharp edges where the velocity is infinite. If 7 be distance 
measured normally from a sharp leading edge where the velocity is infinite, 
then the boundary condition (ix), which requires the aerodynamic force 
to be finite, limits the velocity at almost all points of the edge to be 
~ Kr-*+o(r-*) as r-0 where K is a finite quantity independent of r. 
If ys be the angle made by the tangent to the edge with the undisturbed 
steam, then it can be shown (cf. Appendix to ref. 13) that the contribution 
to F per unit length of edge is 


7 p,y K?/(1— M? sin*s)! (45) 


directed normally outwards from the edge. This force occurs only on 
leading subsonic edges, where M?sin*s < 1. 
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A similar analysis can be carried out for the moment T of the aero- 
dynamic pressures about an origin of radii vectores R, and leads to the 
results: ; : 

T | (p po)R n* dS Po (d*—d)RA(U ds), (46) 


v * 





T = Po | (U.y iv.w)R \n—(U+w).nR, v} dS- 


 (6*—)R/ (U,Ads), (47) 





Po 
T* 


T = p, | PRA(UAds)—p, | (RA Vw.n—4v.wR (An) dS- 


Po | (6*—d)RA(UAds), (48) 
I 
po) in (46) being the quadratic approximation for the pressure. 
It is stressed again that the identification of these expressions (41) to 
18) with the aerodynamic force and moment requires the assumption 
bout the contributions from the flow near the edges of the vortex sheet 
to the integral over S* which was made earlier in this section. 
For later references the following results are required. The first-order 
terms in the expressions 14) and (48) for the force and moment are 


espectiv ely 


7. = VU ( dds Po ( U.vndsS ( (p,—p,)n ds, (49) 
I >> >» 
nd 
T, = pp | ¢RA(U Ads) p, | U.vRAndS = | (p,—p,RAndsS, 
. : (50) 


where [’ is the trailing edge of the mean body surface Y, and (p,—-pp) is 
the pressure given by the linear approximation to Bernoulli’s equation. 
[hese are easily derived from (44) and (48), the difference between the 
ntegrals over ['* and T° being of the second order. That part of the drag 
force (in the direction of U) which neglects the leading edge suction force, 


7 


), say, is given by 


UD, = py, | U.vw.ndS pp | U.vU.n* dS 
U. | (py—p,)n* dS (51) 
since W.n = V.n U.n* on &, and where (p,—p,) is the pressure 


given by the linear approximation to Bernoulli’s equation as before. 
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5. Some integral identities for linearized compressible flow Von 


Consider the quadratic terms in the: momentum integral: —* 
| (vw.n—tv.wn)d\J. (52 
; ‘fs, “, a (TT ; The 
If v and w satisfy the linearized equations (17), (18), and (19) and are 
maa ks ’ supers 
bounded and continuous on, and inside, a closed surface NS, then it can be ip 
5 os ° ’ ‘ (1) 4 
shown that this integral (52) vanishes over S, as has been mentioned 
before, in section 4. It will now be shown that this result is still true if : 
“ 
in supersonic flow, v be discontinuous at a characteristic surface of the - 
linearized equations. The integrand in (52) may be written 
0 
Ivw.n+twA(vAn), (53 
and if the discontinuity in v is Avn (nv is continuous), then the dis 8, 
continuity in this expression (53) is 
\ 
sAvnw .n+$Aw (A (v An) = $Arnw.n+}Aw.nv—tAw.vn "2 
sArnv.® .n—ftArn.®.vn = 0, (54 R 
“3 
since ® is self-conjugate, Aw = Avn.®, and Aw.n = 0. Thus 
‘ On 
| (Vvw.n—j}v.wn) dS = 0 (55) 5, 
. x the e 
oF 
ae — a edges 
even if S encloses characteristic surfaces on which v is discontinuous. |. | 
eae tion 
This is the first integral identity. 

The second integral identity is obtained from the first by putting v 
successively equal to v,+-v, and v,—Vv,, where v, and v, both satisfy the ~ 
linearized equations of flow, and subtracting the two identities. Since Th 
V,.W,. = V,.®.vV. = Vo.W, this gives 

1- V2 1 2 2° Wy sides 

| (Vv, W..n+Vv, wW,.n—v,.w,n)dS = 0, (56) | Integ 
R Th 
which is a relation between two different flows. In supersonic flows, V, ondi 
and v, may be discontinuous on characteristic surfaces inside S without 
affecting the validity of the result. 
6. Uniqueness s R 

There is. at most. one solution of the linearized equations (16), (17), (18), 
and (19) of section 2 which satisfies the boundary conditions (i) to (ix) of 
section 3. 

Let v, and v, be any two solutions, different if this be possible. and let Be 
v= v,—V,. Then v satisfies the linearized equations (16) to (19) and the since 
boundary conditions n.v = 0 on the mean body surface ©, U.v — 0 on 


mean free boundaries, and the conditions (iii) to (ix) of section 3. 
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Consider the identity obtained by multiplying (55) scalarly by U, with 


v=V,—V~ 
| (U.vw.n—iv.wU.n)dS = 0. (57) 
59 : 
“_ The proof is now in two sections, (i) for subsonic flow (M < 1), and (ii) for 
ys supersonic flow (MV 1). 
si i) Subsonic flow 
loner 
= Take the closed surface S in equation (57) to consist of five parts: 
f th © = both sides of the mean body surface plus small tubular surfaces 
surrounding its edges 
+= the surface of the vortex sheet plus small cylindrical surfaces 
si surrounding its edges 
e dis §, = a cylindrical surface of large radius R, whose axis is parallel to U 
and passes near to the body, 
§, = a plane R.U R,U at a large distance R, upstream from the 
f body 
ie §, = a planeR.U = R, U at a large distance R, downstream from the 
body 
7 On XX. w.n v.n 0, and U.n on the mean body surface, and 
°°) the contribution from the tubular surfaces is non-zero only on the leading 
i lves from boundary condition (iii). By using the result (45), the contribu- 
ous oe 
n from & is * aK? cosy ds 
l ment atta. (58) 
ing \ . (] M? sins)! 
y the the integration being over the leading edge. 
oon The integral over o vanishes since the contributions from the opposite 
ides of the vortex sheet cancel, and the contribution to the momentum 
(56 tegral from the cylinder at the edges is perpendicular to U. 
The integrals over S, and S, are O(1/ RZ) as R, > «© by the boundary 
ts ndition (vii), and, from the condition (viii), the integral over S, is 
thout Om < 
= (2+-42) dS+ O(1/ R?) (59) 
2. 
sR, oo. On proceeding to the limit R, — 0 the identity (57) gives 
(18 


P 1T “2 08 Us ds ’ i 9 9 y 
x) of aul eee Ll (¢? +42) dS 0. (60) 
(1— M?sin?y Fe ‘ 


Ss 
id let | Both the integrals in (60) converge by the boundary condition (ix), and 
the . : was 
d th ‘ince the integrands are essentially positive, 
0) oT 


K 0, and 4¢, d. 0 on Ss. (61) 
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Therefore v is continuous on S, and there is no vortex sheet in the differ. 
ence flow, so that Vv is continuous and bounded everywhere except possibly 
at leading edges where v ~ o(r-*). An application of the divergence 
theorem for dw with V.w = 0, namely, 

| dw.ndS = | v.wdr, (62 


to the space between the body and a sphere of very large radius now 


proves that v = 0 everywhere in the same way as in the usual proof of 


the uniqueness theorem for incompressible flow (cf. 14), and the uniqueness 
theorem is proved for the subsonic case. 

The presence of free boundaries only modifies the proof very slightly, 
and it is necessary to take as part of the closed surface the mean surface 
of the free boundary. The contribution to the integral of equation (57 
vanishes since U.v and U.n both vanish on this surface, and the proof 
then follows the same lines as before. 

(ii) Supersonic flow 

Take the closed surface to consist of X, o, and S,, as in the subsonic case, 
plus a closing characteristic surface S, which is the envelope of the down 
stream characteristic semi-cones with vertices on the leading supersonic 
edge of the mean body surface. The contributions from & and o are the 
same as before, and the contribution from S, for finite R,, not necessarily 
large here, is . 

LU | {( M? 1)? pe - p?} ds. (63) 

The contribution from S, vanishes, since v = 0 just upstream from S, 
and the integrand in (57) has been proved to be continuous at characteristic 
surfaces in section 5. Thus the identity (57) gives 


> 


—$U | ((M?—1)¢?+¢2+¢2} dS 0, 64 
(1—M2sin2ys)? * \ IP: $y p:} ¢ (64) 


_ [ wk?\cosd ds 


the integration being over the leading subsonic edge in the first integral. 


Again both these integrals converge by the boundary condition (ix), and, | 


since the integrands are essentially positive, 
K 0, and 4d, ¢,, d. 0 on Ss. (65) 
The latter result is true for any S, which lies downstream from the body 
so that ¢ is constant over the whole of space downstream from any §;, 
and all the derivatives of ¢ vanish in this space. The uniqueness theorem 
for a properly-set Cauchy initial value problem for the space upstream 
from S, (cf. 15) now gives that v = 0 upstream from S, and hence v = 0) 
everywhere, which proves the uniqueness theorem for the supersonic case. 
The modifications necessary when there is a free surface are the same as 


for the subsonic case. 
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FLOW 


7, A general reversed flow theorem 

Let n, and n, be vector functions of position on the surfaces of two thin 
lies having the same mean surface X, the two bodies being respectively in 
pposite flows for which the undisturbed stream velocities are U and —U. 


ty 


, and Vv, be perturbation velocities satisfying the linearized equations 


(18), and (19) and satisfying the boundary conditions 
n.vV, U.n,,. n.v, = U.n, (66) 


YS. and otherwise satisfying the boundary conditions (iil) to (ix) of section 3. 


tp, and p, be the corresponding linearized perturbation pressures, 
Py Py U.V,, De = pa U.Ve. (67) 


en U. | (p,n,+p.n,) dS 0. (68) 
The reversed flow theorem in this form is too general for its significance 
be appreciated immediately, and some special cases of more immediate 
ractical application are derived in the next section. 

In order to prove this theorem, the identity (56), multiplied scalarly 


U. namely, 
| (U.v, w,.n+U.v,w,.n—v,.w, U.n)dS = 0 (69) 


sapplied to the closed surface which consists of the mean surface X, the 
surfaces of any vortex sheets, the surface at infinity, and tubular surfaces 
fvery small radius surrounding the edges of the surface & and the edges 
fthe vortex sheet where v, or v, may be infinite. 
It is a simple consequence of the boundary conditions (iv) to (xiii) that 
the integrals over the surface of the vortex sheet and the surface at infinity 
mish. The integrals over the small tubular surfaces at the edges of = 
so vanish since where v, is infinite, v, is bounded, by the boundary 
ndition (iii), and vice versa. The integrals over the tubular surfaces 
swrounding the edges of the vortex sheet also vanish; for, if the edges in 
the two flows do not coincide over any part of their lengths, then where 
is infinite, v, is bounded and vice versa; if this be not the case, then any 
ntribution from the integral in (56) is a vector normal to U, and vanishes 
n being multiplied scalarly by U. 


Thus the only contribution to the integral in (69) is from =, and hence 


| (U.v, w,.n+U.v,w,.n—Vv,.w, U.n) dS = 0. (70) 
On XL, U.n = 0 by definition of &, and w,.n = v,.n and w,.n = V,.n 
) that, on multiplying (70) by p, and by using the relations (66) and 


67), the relation (68) is proved. 
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It is to be noticed that p, and p, are the linearized perturbation pres- stead) 
sures, so that any aerodynamic forces and moments on the respective surfac 


bodies which are calculated on application of (68) do not include suction 


é . 2 é or, on 
forces on sharp edges, as has been pointed out at the end of section 4. 
The theorem is still true if the flows include free boundaries on which 
U.v, and U.v, vanish. In thi 
8. Some special cases of the general reversed flow theorem 
I. The simplest case occurs when the two bodies are the same and 
n, and n, are the unit normals to the body surface, so that n, a, = n*. 
Then (68) becomes equat 
U. | p,n* dS U. | pon* dS, (71 
and comparison with (51) shows that this result implies that the drag 
forces (neglecting leading edge suction) are of the same magnitude in the 
two flows. When »& is plane this result was first given for supersonic flow 
; On 
in (3) and (4). Paty 
II. With the same bodies, as in I, let a portion =f of the body surfacet | ‘ 
. ; - ‘ _— oes pv t 
in the flow be given an extra small constant vector incidence a,. Then if 
. — . P F the Cc 
n* was the normal to the original body surface, the new normal is approxi- ’ 

- . roe mr the oO 
mately n*+a, An, neglecting second-order quantities. Take now rV 
n, = @,An on %&,, n, = 0 on the rest of &, n, = n* on x; (72) ” 

l 
then p, is the excess pressure consequent upon the change in the body  q, ar 
surface, and (68) gives vecte 


U. | p,n*ds U. | pa, nds, 


or U. | p,n* dS Una. | ponds. (73) 


N 


On comparison with (49) and (51) it is seen that the left-hand side of 
(73) is proportional to the extra drag force (neglecting edge suction force) Ther 
on the original body due to the extra incidence of 2*, and the right-hand char 
side is proportional to the component perpendicular to U and in the plane 
of incidence of the original pressure in the reversed flow acting on =f. 

III. Again with the same bodies, let a portion =f of the body surface 
be given a small angular velocity Q, about a point whose position vector 


relative to some origin of vectors is R,. In the usual theory of such quasi- 


: ay : ‘ Q 
steady motions, it is assumed that the flow and forces are given correctly | 
. . . : , ; ‘ tion 
in linearized theory by warping the body suitably and considering the | | : 
LO 
+ The reader is reminded that ©* denotes the actual body surface and ¥ denotes the ict 


corresponding mean body surface on which the boundary conditions are applied (see 
section 4). 
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steady flow past the warped body. The boundary condition on the body 


surface at R is 


n*.(U+v) = n*.Q, \(R—R,), (74) 
on the mean body surface ©, by neglecting second-order quantities, 
n.v U .n*¥+n.Q, A (R—R,). (75) 
this case, take 
n.v, U.n, ='n.Q, \(R—R,) on &,, 
n, = on the rest of &, 
a= ar on Z; (76) 


lation (68) then gives 


ad. [ p,n* dS ( p.n.Q, \(R—R,) dS, 
a, ( p,n* dS =Q,. ( p(R—R,)Ands. (77) 


‘ s 


= ‘1 
On comparison with (50) and (51), it is seen that the left-hand side of 
is proportional to the extra drag force on the original body caused 
the angular motion of = and the right-hand side is proportional to 
the component in the direction of the axis of rotation of the moment of 
the original (reversed flow) pressures on Xf about the point R,. 
IV. With any two bodies having the same mean surface ©, let a part 
‘¥ of the body surface in the flow U be given a small extra vector incidence 
a, and a part XF of the body surface in the flow —U be given a small 


ector incidence a,, and take (cf. case IT) 


n, a,A\n on da» 

n, 0 on the rest of &, 

n,—a,A\n on Xz, 

n, = 0 on the rest of &. (78) 


(hen p, and p, are the pressure changes consequent upon the surface 


hanges. and (68) gives 


U. ( Pp, % Ands U. ( poe, \nds, 
U ij. ( p,n ds Us a,. ( pon ds. (79) 


s 


Ys a 
On comparison with (49), the left-hand side of (79) is seen to be propor- 
mal to the component in the plane of the incidence a, and perpendicular 

to U of the force on &, due to the extra pressures caused by the extra 

neidence a, of &,, and similarly for the right-hand side. 


The special case of (79), when ©, = ¥, = ¥ and Y is a plane, and when 
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a, and a, are equal vectors perpendicular to U and in the plane of ¥, 
gives the lift slope reversal theorem of (9). 

V. As in IV, let the area =¥ in the flow U be given a small extra 
incidence a,, let the area =¥ in the flow —U be given a small angular 
velocity Q, about the point R, and take 

nm, =@,An on 2, 
n, 0 onthe rest of &, (80) 
n.vV; U.n, = n.Q,A(R—R,) on %, 
n, 0 onthe rest of E. 
The p, and p, are the pressure changes consequent upon the surface 
changes, and (68) gives 


| p,n.Q, \(R—R,) dS U. | pra, And, 


N 
2 1 


or 2. | p,(R—R,) \ndsS Una,. | ponds. (81) 
x1 

Comparison with (49) and (50) give an interpretation of (81) similar to 
those of the previous cases. A special case of this result, as in case IV, 
is given in (9). 

VI. Let the area S¥ in the flow U be given a small angular velocity Q, 
about the point R,, let the area =F in the flow —U be given a small 
angular velocity 2, about the point R,, and take n, as in (76) and n, as 
in (80). Then p, and p, are the pressure changes consequent upon the 
surface changes, and (68) gives 

2,. | p(R—R,) AndS = Q,. | p.(R—R,)AndsS, (82) 
Ys x1 
for which the interpretation is similar to those given above. Again special 
cases of this result are given in (9). 

VII. The six special cases of the general equation (68) given above are 
intended to be illustrative and not exhaustive; for instance the surface Lf 
in the flow U has been taken to have constant vector incidence in cases II, 
IV, and V, whereas it might have been divided into two parts, say, with 
incidences in opposite directions, thus giving an aileron-like movement; 
many other such cases can easily be envisaged. The theorems are true for 
flow in wind tunnels with either closed or open working sections of suffi- 
cient length. 

[It is worth pointing out that not all the special cases I-VI are inde- 
pendent and some can be deduced from the others by superposition, since 
the pressures p, and p, are linear functions of the fluid velocity. The more 
general cases mentioned in the previous paragraph can also be obtained 
by superposition of a suitable choice of I-VI. 
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of 2) 9 A special reversed flow theorem 
There is another reversed flow theorem of which a special case has been 
extra yiven in (3), and incorrectly in (7), which is not deducible from the general 
gular theorem of section 7. 
When there are no vortex sheets in the flow so that the line integral 
$1) for the aerodynamic force vanishes, there is left the integral over S* 
(50 hich depends quadratically on the perturbation velocities, so that the 
mponent of force perpendicular to U is of second order (the drag force 
salways of the second order, of course). If in addition there are no infinite 
irface elocities in the fluid, then the surface S* may be taken to be the mean 
ody surface %, so that 
= me ( (vw.n—}v.wn) ds. (83) 
81 Now take the mean body surface to consist of planes parallel to z = 0, 
, The lateral component of force perpendicular to U and parallel to 
lar to the mean surface & is , . 
e IV. F po | d,w.nds. (84) 
[he y-component of the vector integral identity (56) taken over a closed 
a surface Which consists of © and the surface at infinity gives 
a. | (f,, W..n+¢,, W,.n)dS = 0, (85) 
n the 


the integral over the surface at infinity vanishing by the boundary con- 
litions (vii) to (ix). Since the boundary conditions on = for the flows U 


/ *) + . 
(9a nd —U give 


necial n.W, = N.V, U.n n.V, n.W, (86) 
equation (S85) becomes 
ee | (¢,,W,-n+¢4,, w,.n) dS = 0 (87) 
ce | 
l >» 
es IT, r, from (84) F, F,,, (88) 
with ry _ ‘ ; : . 
hus the y-component of force is reversed in direction and is unchanged 
ment La ere : , 
n magnitude. This result is trivial for subsonic flow, since all the forces 
ue fo! nee. 
suffi- 
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ai SUMMARY 
- {n account is given of the methods available to improve the accuracy of arith- 
al solutions on a square mesh for two-dimensional partial differential equations. 
ee Squaring’ (1) and * Relaxation’ (5) depend for their accuracy on the fineness of the 
1 which the finite difference equations replace the differential equations. The 
accuracy can be obtained on a coarse mesh, with less labour, by the use of more 
rate difference formulae, rather than by refining the mesh. There are two ways 
ing this: 
by using a higher order difference equation as a closer approximation to the 
lifferential equation, o1 
by calculating a correction term from the higher differences to apply to the 
results obtained by using the first approximation. 
w method of calculating the correction term is developed. 


1, Introduction 
lHis paper describes two methods of improving the accuracy of arith- 
tical solutions obtained by ‘Squaring’ (1,2) or ‘Relaxation’ (5), to 
iptic partial differential equations on a square mesh. The types of 
lations considered in this paper are 
(a) V2 constant, (b) Vf d(x, ¥), (c) V4f A(x, y) 
(d) Vif = Of, f..Fy.%¥), 
here suffixes denote derivatives. The arithmetical solutions are usually 
lound by using the first finite difference approximations to the differential 
equations. The values of some wanted function at the centres of the 
squares of the mesh are expressed in terms of the four values at the corners 
the squares. The accuracy of this first approximation can be improved 
either using a higher order difference equation as a closer approximation 
the differential equation (4), or calculating a correction term from the 
igher differences to apply to the results obtained using the first approxi- 
ation (3). 
Che first method results in formulae expressing the value of the function 
the centre of mesh points covering several squares. Bickley’s term 
molecules’ (4) will be used occasionally to describe the square or sets of 


(Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 3 (1950)] 











L. C. WOODS 





350 


squares over which the finite difference equation spreads. This paper gives 


the large molecule formulae of Thom and Bickley, with some extensions, | 


In the second method two ways are available for calculating the error 
terms. For example, the largest error term, usually the only one of any 
importance can be obtained by 

(a) calculating the expression Ay = ;4(53,+-8%,,), where 93., 


fourth differences of the function f at the point 0, in the OX and 


64, are the 


OY directions respectively; this is Fox’s method. 


1 asf | 
4 a : Cc ° ° 
(6) calculate instead A, = Vif—2 J .}, and inserting the value of 
12 Ox"0y* 
V‘*f from the field equation, e.g. for Laplace's V4f = Oand A, a aaa 
) cox-oy* 


where a is the mesh size. The second of these two methods, which involves 
the use of the field equation in the calculation of the error term, requires 
a spread of only three grid points in each direction, as against five in Fox’s 
method. The method, which is developed in section 4, proves easier to 
apply numerically, avoids extrapolation near the boundaries, an undesir- 


able feature of Fox’s method, and provides a criterion of suitable mesh size. 


2. The finite difference formulae 


Adopting and extending Bickley’s notation (4) we write 


eS ’ Beas , 12 
S; hh +-fo 1 Js Hh di So za Ss DSi 
5 9 
20 24 e 
S, =), S, > f., Pf _—— . 
T3. De Cxcy 
where f; (2 1, 2,...) are the values of a wanted function at the lattice 
points shown in Fig. 1. Bickley then shows that 
ae . : 2H” om — 4 2 
S; \4 -ar?Vv?2 + —_(Vi— 294) 4 re (V§—394V?2) 4- 
t . 
9a8 _ ie A 
+— (V8—4F°V44- 29%)+...) fo, (1) 
8! 
—. 46 4q° s 
S, 4 | 2a2V2 4 1 (V4+-474) 4 a (V8+ 1294V2) + 
4! 1 
4a’ 
(V8 + 24044 169%)+...1for (2) 
S. 
; aosa ; eee 128a* .... = 
S. I | 4a2V2 (V4—2974)4 (V§—37°V?) + 
3 ' al ‘ 
| 4! 6! 
512a8 


+= (V8—4G4V44 298) +... fos (3) 
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and similarly it can be shown that 


at 4a 
8, = [8+ 100°V2 4 2” (1704-41498) 4 (6508+ 10598) + 
4! 6! 
4" (257084 876F4V4— 105498) +...'f,, (4) 


8! 


> 44 256a° _. 
S, = [4+ 802V2 4% (yey sty 4 OO (yes pogayey + 
4! 6! 
1024a® 


5: 


(V8 2A AVIS 16) + fo (5) 


22. 1 '0 4 2) 





16 6 Ms 15 13 
_— 


17__‘(|7 4 18 20 

















the special case V2f = constant k say, then V4f = Vf = 0, etc., and 


aV2f 0, ete., and so (1) becomes 
as . , 408 . 
oF 4 fy a*h : D*fo ' PD f+ O(a"), (6) 
6 , S! 7 


rom which follows the well-known first approximation 


4 . ~ 
$7, a*k-+- > Ie. (7) 
1 


tquation (7) was employed by Thom (1) in 1928, in his method of 
squaring a field’, and by Southwell (5) in 1938 in the related process 
termed ‘relaxation’. 


The accuracy obtained by employing (7) on a given mesh can be 


uproved by four distinct methods :— 


a) Reduce the mesh size a, which, from (6), obviously reduces the size 


ithe error term. This, which is the normal procedure in relaxing or 
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squaring, is repeated until a subdivision does not alter the values of { 
found at the lattice points of the previous mesh. 

(b) Use a larger molecule obtained by use of (1) and (2) to eliminate the | 
a*@* error term (four squares); or use of (3), (4), and (5) to eliminate both 
the at74 and a8 terms (sixteen squares). All the V? error terms disappear | 


of course, as in (6). 


(c) Caleulate the largest error term using Fox’s difference method, and 


correct the results for this error. From the general equation (1), the 
largest error term is 

1 a4 (ef, ert 
= (V4 2F")f, | JO | =*) 


12 12\dx* ay! 


2 9 (Oke -§4,)+ O(a%), 
but, as already mentioned, the calculation of the term in this form does 
not involve the field equation. 

(d) If the field equation is used, in this case V?f = k, then V4f = 0, and 
the error term reduces to —(a*/6)7#f, as in (6). 

The larger molecules used by Thom amount to a repeated application 
of (7) over a mesh of four or sixteen squares, carried out in one operation. 
The a* error terms are not eliminated, and thus it is essentially a ‘speeding- 
up’ process. However, accuracy is in fact improved, because the possible 
cumulative errors, arising in the application of (7) square by square to the 


mesh do not occur. 


3. Large molecule formulae 
(a) V?f = constant z. 


Thom’s formulae (2), with the addition of the largest error terms, are: ' 


4 squares: 12f, = 28,+-S,—4a*k—(a'/3)74f,+ O(a) ; (8) 
16 squares: 476f) = 46S,+-32S,+-9S;—576a*k—48a1@4f,+ O(a). (9) 
Alternatively: 
' 22 
L00fo LOS, + 78S,+S;+ 118a°k — — at Gf, + O(a’), (10) | 


~» 
an equation with convenient coefficients, and incidentally a lower error 
term than (9) which is based on (7). All these may be verified by adding 
(1) to (5), putting V?f = /, in the correct proportion. They are particularly 
useful for rectangular fields with given boundary values. 
Eliminating a*Z@*f, from (1) and (2) yields Bickley’s ‘four-square’ 
formula 


> . 80 ee 
20f, 1S, +S,—6a2k ai S*fo | O(a}?), (11) 


In a similar way a*Z*f, and a&Zf, can be eliminated from (3), (4), and (5) 


to give an expression with an error term O(a!2), but with very large and 
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impractical coefficients. The high accuracy of such an expression would 
be generally unwarranted. 

b) Vf b( v,Y). 

In this case (1) becomes: 


r 0 


a! ; seins 28 om ; 
Sy 4) +a" by 4 15 (¥ "Po 2D*fo) + 6! (V4¢9—3F4h9)+ O(a'). (12) 


from the corresponding equations ‘for S;, S,, and S; the a* term can be 


eliminated to give 
oF ' Wn2h 1 6 ‘ 
12f, 8S,—9S,—4S,— 12a*o_)+- O(a). (13) 
fin addition V#f constant, i.e. V2d constant B, say, the a* and 
terms can be eliminated to give 


7S; — 1560a*6,— 456 B-+- O(a), (14) 


848, f 1288, 


result of little practical value. Bickley’s procedure implies that we should 
se the form corresponding to (11), i.e. 


20f,, 1S, +S,—6a*d,.— 4a! V4f,+ O(a), (15) 
hich eliminates the a4 term, to calculate an approximate f,, by ignoring 


the a#V4 term. This term is then found approximately from 
> wit 9 4 
atv "fo a*(a*V"do) = a*(> $;—445), 
I 


nd applied as a correction term. Bickley does not state this directly, his 
rk being expressed as if the computor desired to calculate the values 
fV* from a given f distribution, rather than seeking an f distribution 
) satisfy a relationship involving V?f. He also incorrectly states that 
the above procedure is superior to Fox’s difference method, since ‘it is 
yo-dimensional in character’. Bickley has used a larger molecule to 
eliminate the a4@* term, and a correction term to allow for the a*V‘ term. 
[his is equivalent to (13) which eliminates both terms. The difference 
ethod is exactly equivalent in its elimination of the largest error term, 
treating both the parts a*7* and a*V# as correction terms, since 
, (88, +-55,) . ( Io a) “vy “9 ‘fo: 
12 \ art" dy} 12° °° «6 
c) Vif = d(x, y). 
The first approximation can be obtained by two methods. The first 
as introduced by Thom (2) in connexion with viscous flow. This consists 
husing an intermediate function C(x, y), so that V2f = f, and V7Z = d(x, y). 
When the ¢ field is settled (or relaxed) with approximate boundary con- 
litions, the f field is settled from V?f = ¢; this alters the boundary 


nditions for the ¢ field, which then needs to be recalculated. The fields 


092.1] 
«il 


Aa 
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are settled alternately in this way until the change from one round to the 
next is negligible. The final accuracy can be improved by any of the; 
methods tabulated in section 2. 
Thom also established the fundamental equation of the second method, | 
which can be obtained by using (1), (2), and (3) to eliminate V?f, and 9f,: | 
20fo = 8S,—28,—S?+-a4V“fo+ Ja°V “fo— sa DAW "f+ O(a*). 
The error term, ja°(V°f,—2F4V?f,), could be computed by methods similar 
to those in section 4 below. Alternatively (4) could be used to eliminate | 
this term, but the resulting molecule would be cumbersome. | 
(d) V2f = Of. fesfy,%Y): 
For special simple cases of this general equation, larger molecules could 
be derived, but as the elementary molecules themselves are often complex 
enough, the practical procedure is to apply the correction term as described 
below. 


4. The correction term method 
When V?f = (2, y) as in (12), then 


4fo SS —a*dy ae . (V*bo— 2D fo) — = (V4do— 3D4*$9)+ O(a’), 

(16) 
which covers the types of equations (a) and (b) of section 1. From 
(1) and (2), 
as 


_ as ‘ 4 
ag0 Jota 2 do +- O(a"), (17) 


6 
a 
Y 99 4 IGAf PAV 
S,—28,+4fy = atD*fy+ T toad V*fot+ 
a result which holds irrespective of the value of V?f. Hence 


494f — § 99 14 as | 4 - Av2 a® - 8f | O(ql0 
UL Jo ——S = 1 fo 12” do— L Po— aan BD fo 1 (a im 


360 
Also from (1) applied to ¢, 


Oh wh 4 a‘ ‘ 
a°v"by p $; a 4g¢— 12 (V4d)—22 ‘dy) + O(a"). 


Using these in (16) gives 


py 4 9 Y Y ! a* $. ! 
4fy 2 fi— Oo T (Sp is 2; +4fo) ao 12 >: $; ~~ 4h) i 
6 Ly f 
eee .— a& 44,4 O(a*) (18) 
' 240 360 


or in the special case ¢ = 0, or a constant 


as 


G8 12 9 
3004 Df.+O(a'?). (19) 


4 
4fo 2 fi— 0b t+ §(S3— 28; 4fy)— 





(18) : 
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othe (4g) and (19) will be written in the form 


f the ‘ 
4fy = > fi—4*bo+ Ant 8ot O(a*), (20) 
thod , 
Df, _— i a? /4 
I where Ay = 3(S2—28,+ 4fo) — 12 | > $;—4o), (21) 
2\T 
‘milar | a’ v4 7 yegs 25 
| and 3, : V4do— ag 22 
inate »= 249" P0— 360% Po on 
or if 6 = 0, or a constant, 
—_— , " 
could Ao §(Sp— 28, + 4fo). (23) 
a’ . 
_— By = — 55 Yo (24) 
ribed 3024 
Ignoring the secondary error term 4p, it is required to solve 
" 3 > 9 
4fy = > fi—a*hot Ad. 
1 
90 . P ° ° bina 7 4 7 9  * 
, If f, is a first approximation satisfying 4f, = > f;—a*¢o, A; is calculated 
T 


(16 from (21), i.e. 


Fron [ep ioee A 6M sh. 
j Ao 6( > Ji - Ji T 4f5) * 19 (2 $;—4h9)- 
5 1 ~- 


4 
), (17) Thesecond approximation is found from 4f5 = > f7—a*¢ y+ Ao, or putting 
1 
, , . ; a , , . . 
fi+g), the equation 4g, = > g;+ A) is solved in the usual manner. 
] 


~ 1 4 
. We now put Aj = 4(> g;—2 ¥ g}+4g,), and find 495 = > gj+Aj, etc., so 
5 1 1 


that finally, f, = f;+gi+g;+..., and A; = Aj;+Aj-+.... The boundary 


values for f; are the same as those for f;, while g; = gj =... = 0 on 
these boundaries. The method converges rapidly and often it is unneces- 
sary to proceed beyond g;. The method differs from that of Fox by the 
way in which A, is obtained. 


Equation (23) can be written 


5 - = 2 
6A, (S fi— dfi)- (> fi—4fo) (25) 
(18 5 J 
8 ic 
(Sfi- Yfij- Ro, 
5 1 
(19 vhere ‘R,’ is the ‘residual’ of the relaxation process. If residuals are 


known then the calculation of A, is simplified. The author’s method is to 
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check the residuals at the end of the usual relaxation process, and to 
compute 6A, at the same time by setting up the difference table: 


4f, 
cers 
(fi 4f,) 
a Bi 4 4 
afi (Dd fi 2 Si) (2fi- 4f,) ( 6A,) 
8 4 
(> fi > fi) 
5 1 
Sf, 


Assuming that the residuals need checking in any case. this method of 
calculating A, is quicker than Fox’s method of working out the fourth 
differences, particularly when ¢ is zero or a constant. If 5, is required, it 


4 
can easily be obtained by noting that a®°V44, = a?(¥ a?V°4;—4a°V"4q), 
1 


Ax 


4 
and a°'*, = a*(> ¢;—2 > $;+4¢ 9). a?V74; is computed when finding A,, 
1 


ot 


If d constant, 


ag : ; (oe 3 

8p ano, Sef, and a& Df, = atP*(a*D")f, = 6(> A;—2 > A;+4<A,). 
3024 ' 5 1 

(26) 


5. Boundary conditions 
In this section methods are given for calculating A, on the boundaries 
on which are specified (a) boundary values, (6) normal gradients. 


(a) Values specified 


Sas of fy, 
since —" 4 5 (2, 9), 
on" Oy* 
if a2 nlf a2 pat 
o*f o* of ap af - 
then —. = —|¢ - awe (27 
Cx-oy* Ox" Pc Ce" o£ 
asf 24, otf 
1s . C C Cc ¢ 
Similarly I _=— I (28) 
Ox*dy2 ss Oy” _—s oy" 
These equations give the values of 7*f along boundaries x — constant and 


y = constant. For irregular boundaries (5), values at ‘fictitious nodes’ can 
be obtained by extrapolation in the usual way, although it is often possible 
to avoid ‘irregular stars’ by a conformal transformation (2). This is, of 


course, not an algebraic transformation applicable only to special cases, 


but an arithmetical transformation applicable (at least in principle) to all | 


shapes of regions. An example very suitable for this type of approach is 

. ‘ . bors Bie ’ ; = 
the incompressible flow (V27 = 0) about an aerofoil. The calculations ar 
carried out in the square mesh formed by the velocity equipotentials 
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id to » — constant) and the streamlines (¢ = constant). The aerofoil boundary 
i, a slit along » = 0. The boundary conditions, e say, are specified as 
functions of position s (¢ e(s)) on the aerofoil, but from an assumed rela- 
tion s = 8(7), we can find approximate boundary conditions for the (7, ¢) 
plane. € = e(7), and hence a relaxation solution can be found in this plane. 

| By an integration along the boundary this solution yields a new relation 

‘g) . , . . . ‘ ’ 

_ = 8,(y), and hence a new relation e, = e,(7). Relaxation yields a new 
solution in the (, ¢) plane, and hence a new relation s, = 8,(7), and so on. 
[his iterative process converges rapidly. 
If d constant, or varies slow ly, 
1 1 pa 
a " a* Cc 
od of A, DY. Jo. 
“baat, and 
: 6 6 C2 
ourt! one 
Lj e, A, oc 54.. where 44, are the fourth differences along the boundary. 

"ed, If . . 
[hus for constant accuracy, the mesh length @ must be adjusted so that 

was ais ' 

V"@p) the fourth differences remain approximately constant along the boundary. 

[his proves a useful guide for mesh size variation near the nose of an 

o A 

1B Ao rofoil, sharp corners, etc 

(b) Normal gradients Spe cified 
From Fig. 1 
4A,). (i F Corr . o*\ . 
Js lo a| \ fo i at@ T ° To To sees etc., 
(26 C2 cy) ai\Ox* CHCY cy 
nd thus 
2 2f 3 ast ast 
: f / of oy lo a2 ( ly To | a 6 c to | 9° To 
daries 16 J0 ey ex? ' ey?) ' 31\ ax®ey cy? 
3f 3 ASsf i 
Ch, —— aad | a? ¢ a a: : ‘ 
2] 4 r+-a*V 4 a®_—_-"_ 4 So += (V fot 4F*fo) + O(a*). 
J0 0 9 ‘ 3! 5 ‘ j I 
cy cu-cy 3 cy 12 
(29) 
. 2F . 4 o4f 9776 pet 
Asi , f | 2 12° fo , 2a* O*fy , 2a* Af 
sv } 13 a "he 9 | T > . 
= :, ; =" Se oa 4! dat 6! o26 
(2) 
| / / “Jo a* e*fo a® Hf 
ni ) a ' 9 | « - 
(28 7 cy 2! dy? 3! oy* 
. Hence 
3 p3f 4 
| Cc] ee a” ¢ a ; , . “ 
nt and I Is 2fs t/ 2a - aX ‘Jo ‘ Jo T (V fo 22 *fo) T O(a*). (30) 
13’ Cal ‘is cy 3 cy 12 
3 of 1 
a é a —— 
, ) ‘ D4 i of = : J0 : >» 
yssible From (29) and (30), f,+-/3+2f.—(fs+/6)—2fo = as +0 _ = D*f., i.e. if 
. f Cx*0y 2 
; 1S. 0 , . 
Jo, OF . . . 
cases : » — db», equation (30) can be written 
- - cu- 
) tO al 
. 9 ») CJo aA 
ach 1s 4], hy Is ole al a-Do 
= cy 
S are 37 4 
yn A a Cdby a “ 


entials 
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which corresponds to (18). Thus in this case 


i oo om a3 6 B 0 
Be = Hetle— (ithe +2fe)+%o—F PT Vy | 


the last two terms of which can be calculated when ¢ = ¢(x, y) is specified. | 
If V?f = 0, a large molecule formula corresponding to (11) can be found 








by substituting V?f = 0 in (29) and (30). This leads to 
é 3 
fstfet+2(fitSst 2fe)—10f,—6a Io + O(a*). i. 
cy 
g2 
a 2 Dy 
a 
64 
Fia. 2. 





(c) Irregular stars 

If irregular stars are encountered they can be handled by a method 
described by Southwell (5). The results, however, have an error term of | 
O(a?), whereas the error term away from the boundaries is O(a‘). Extra- 
polation formulae giving greater accuracy near the boundaries are 
unnecessarily complicated (3). The following method is relatively simple 
and reduces the error term on the boundary to O(a). 

From Fig. 2 we have 


faths vail see 
iy 
fo ' pao ia fe 


"7% 


f,: 5 + O(a), 


, f,,a eh 
fs = fo—a + 2 —s + O(a). 


a : te ea ; 
Hence Ath th +fs (3 >" i)los av *fy+ > (b—a)}y+ O(a’), 


from which follows 


a+ a (a+b)? 


(fo +fa)+ fitts a 


ig. 3 shows the alain relaxation pattern for 6 = 


fo +5 (b—a)bo+ O(a*). (33) 
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6. Advantages of the method 

1. The calculation of fourth differences requires a spread of five lattice 
| points in each direction to obtain one central value ci the primary error 
term, which forces the computor to extrapolate near the boundaries. 
Equation (21) shows that the present method requires a spread of only 
three lattice points in each direction. 
2. The large molecule formulae; e.g. (11), are very useful under certain 
(32) | conditions, but it is usually quicker to keep to the simplest relaxation 
pattern and use the correction term method, particularly when working 


sified. | 


found | 


through graded meshes. 

3. It would be possible to caleulate and include the second error term, 
§,, by the method described in section 4. The calculation of 5) by the use 
of sixth differences would be possible, but clumsy. However, near the 
boundaries the excessive extrapolation required would give doubtful 
results, a disadvantage not possessed by the author’s method. 

4. For a given number of significant figures a solution of (7) can be 
improved by either 

(a) advancing to the next grade of mesh of side a/2 in the normal 

manner, or 

(b) applying the first correction term as described in section 4, to the 


soOarTse sn. 
othod coarse mesh 


— The question arises which of these results will have the greatest accuracy. 


~— The methods of Fox and Bickley give no direct answer to this question, 
ie which is certainly important if it is decided for some reason or other that 
imple | results are required on a finer mesh. If, for instance, it can be determined 
that method (b) gives the more accurate results, then during the advance 
to the finer mesh these accurate results could be maintained constant, 
just as if they were fixed boundary values. This would obviously con- 
siderably reduce the work involved in the advance to the finer mesh. 
Of course, if the results on the finer mesh are then improved by use of the 
correction term, the values maintained constant during the advance must 
be permitted to change in the usual way. 

If in (16) a is replaced by a/n, each mesh side is subdivided into n parts, 
and on the resulting finer mesh the first error term is A,/n*. If, on the other 
hand, A, is eliminated on the coarse mesh, the magnitude of the remaining 
error term is 6). When 

A,/n* = 85, ie. m = (Ag/do)*, 
(33) then n can be termed the ‘equivalent mesh size reduction’. Such a mesh 
size reduction, in conjunction with the ordinary methods, is equivalent in 
accuracy to use of the correction term on the coarse mesh; 5, can be 
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calculated at a few points from equation (26). The values of n at these | Thu 
points will give an indication of the comparative accuracy of the correction | cort 
term method. For instance, if n = 2, the values found using the correc- (28) 
tion term could be maintained constant during a reduction of the mesh 
length to a/2. This is illustrated in example (i) below. 

Equations of the form V4f = constant are amenable to a similar treat- 
ment as can be seen from (16); but equations of the type V?f = 4(f.f,.f,,.7,y) Her 
are too complex to receive a general treatment. Simple cases, such as 
V*f+Af = 0, can be treated by the author’s method as illustrated by Thu 
example (ii) below. 





mes 
7. Examples mn 
mr . . + . . . 3 
rhe following two examples were treated by Fox and form a basis for ;' 
comparison. In this section, as in the previous sections, subscript ‘0’ will ng 
be used to denote a central value, subscript i (¢ |e 4) will indicate vale 
1Y 
| Fl ' 
O| 0 
768i y , , 
= Of 58 
3 / ¢ So 
13 12% «A ‘7 
107| 72 g \4. 
640|9427 FLA X 
a eo 75 
79 
Fic. 4 = 
values on the four nearest mesh points, and i (7 = 5,...,8) will indicate 
f values at the next four mesh points. 
(i) Solution of V2f 2560, with f = 0 on x +l,y +], ( 
By symmetry only one octant need be considered, although the figures the 
are drawn for one quadrant for convenience. The first approximation is | obt 
: val 
4f, > fi +2560a?. (34) stal 
i 
. . « : 7 > tl 7 be. 
(a) a 1. Equation (34) gives fp 640 (see Fig. 4). Fox was unable : 
to improve this value. From (23) 
: 4 640 b> ; -_ , 
Ao - 427, andthus 4g, = >} g;+A, = 427, as g;=0. | Th 
) i 
427 g 4 ent 
3 a , - = j =. ” 1 >) , 9 , A) 
Hence 4g ; 107. Again Aj 6(> gj; “2294 495) 12. rels 
vo 
wh?) 4 x 18 
” la ” < 
Hence gj - = 18, AG == 12. 
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Thus fy = fo+90+90+-90 768, which is 1-8 per cent. in error from the 
correct analytical result. A value of n can be calculated, for using (27) and 
98), A, = 0 (¢ = 1...., 8) on the boundary, and hence, from (26), 

it £8 ded © (4511) 
Oo 2 = xX oo ° 
. 3024 \2 0] 3024 
Ay\t . . | 
Hence 7 | -3 (taking the lower nearest integer). 
O 


Thus the value 768 is more accurate than that obtained by dividing the 
mesh length by 2, and employing (7). In fact, by using ordinary relaxation 
na mesh a = 3, the error is 4:7 per cent., and with a = } the error is 
|-2 per cent., which confirms the intermediate value of n = 3 correspond- 
ing to 1-8 per cent. error. Thus we are justified in retaining 768 as a fixed 
alue during the advance to a . 


| | 























0 O 6] 
0 0 O T 
, 587\27  464\44 
P 1| 1 
587 464 260 23) 4 
560 440) 560\27___ 440/40 0 
576 448 O . Bl 
754;30 587.27 
C4 — ; u 0 
30 56 33:3 26 
Je 56 | 720lé7___ 560\27____— ol 
68 | | ae. Sn A 
Fic. 5. Fic. 6. 
(b) a 1 Fig. 5 gives the results obtained using (7) on this mesh, with 


the central value fixed. The middle entries at each point are the results 
btained by using ordinary methods (equation (7), without the central 
alue fixed), while the top figures are the correct values. In Fig. 6 we have 
started with the results found by ordinary methods so that the values can 
be compared with those found by Fox. 


Dealing with point B as an example: 


Ap lf4 x 440— 2(560+560+-0+-0)+(720+-0+-0+0)! = 40. 


( ; 
This calculation is repeated at each lattice point, the values so found being 
entered to the right of the f values. The value Jp is found by the usual 
relaxation process, giving g’, 23, shown entered above 440. Again 

A = 44x 23—2(0+ 26+ 26+0)+(33+0+0+0)} = 4, 
which, together with the other A’ gives rise to gj, = 1. The final value 
tf; = fit+gi+..., A; = Aj+At+..., are shown above the line at each 
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point in Fig. 6. As before, values of n can be calculated and they are 
found to be n, = 3:7, np = 3:8. Thus these results would be quite 
accurate up to a = 1/(2 3-5) = }. 

(c) a = }. The results given in Fig. 6 are now entered on Fig. 7 (shown 
in squares) and maintained constant. Intermediate values are entered 
using (34), and the only large residues arising are shown in circles. Because 


| 
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L 0 0 
| | 
| | 3 
| __19642)_0 184|7__0 
| 
| | foal 2arl-1 0 464|-6 287/30 
re | t 
| | 
| | | og / 
678(-8)558\0 3435) 0 |_ 675|3 5562 340/20 
| y 1 0 
54| yay fee7|_3ss\¢ ol Ss 754L:5 7i2|2_se7tS 3570 Oo] 


Fic. 7. Fie. 8. 

of these fixed values this field is easily relaxed and the results are given 
in Fig. 8. The values of A,, plus small residuals remaining from the 
relaxation obtained from this field, are shown to the right of the lattice 
points in Fig. 8, and the corresponding values of g; are entered above the 
f;. The solution f; = f;+g; differs from the analytical solution quoted 
by Fox at the most by unity in each entry. At the centre n is now greater 
than 3, an indication that the accuracy on this mesh is at least as good as 
that obtained by ordinary methods on a mesh size a = 1/12. 


(ii) An Eigenvalue problem V2w+dAw = 0. (35) 
If w is an approximate solution of (35), Rayleigh’s Principle gives 


( [ wV2w dady 


A= -— (36) 
| w? dady 
, _— ‘ as 
From (16) a?V2w, y w;,—4wy)— Tha —2F*)w)+ O(a*), 
but V4w = —AV2w = A*w, and thus 
9 9 ‘ at 9 GA = 
&V2wy = > w;—4uy— 5 (A2w y>—2F4wy), (37) 
1 - 
4 € ) 
or, less accurately, @V2wy = > w;— 4p. (38) 
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are Starting with an approximate solution to the mode w, V*w is calculated 
juite from (38) and A found from (36). This value of A is used in (35) and from 

aw) > w,—4w, = 0, a new w field is calculated. From this a new A 
Own 1 
ered is calculated, and so on, the process being repeated until no further 
ause improvement is possible. 
In the example given by Fox the problem j 
isto determine the fundamental vibration | 
of a square membrane of unit side. The | -6 74 (-620) 
analytical solution is | -57 
Ww = COS 7X COS 7rY, A lll = 19-74, | 500 Lote 
Fig. 9 shows one quadrant, with a i, 4 
giving the first approximation to w, ob- [1229 (1240) -868 (-876) 
tained by the process mentioned above. Ke 99 
The values of a?V2w found from (38) are /000(1I72 _707|-828. sj; 
shown to the right of the lattice points. Fic. 9. 
The corresponding value of aA calculated 
from (36) is 1-172, i.e. A 18-75. Then from (37) 
seleees 4 
eV2w = ¥ w;,—4wy—0-115w9+ fat Fwy. 
given I 
n the The two correction terms, ‘—0-115w,’ and ‘}a*@4w,’, are easily found 
attice and are shown in Fig. 9 above the approximate values of a?V2w. Using 
re the these new values of a*V2w, in (36) we find a7A = 1-226, therefore A = 19-62. 
uoted With this new value of A the correction term ‘0-115w,’ becomes ‘0-125w,’, 
reater giving slightly altered values of a?V?2w;, shown in brackets in Fig. 9. 
vodas Using these we find A= 19-87. Repeating this process again gives 
\= 19-94, which is 1-0 per cent. above the correct result. Fox’s result 
(35) is \= 19-66, which is 0-4 per cent. below the correct result, but since 
Rayleigh’s Principle should over-estimate A, his result is not necessarily 
- the more accurate. 
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SYSTEMS 
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SUMMARY 
In this paper we discuss the differential equation 
t+ wr kf (x, #, $i, 
where w, k are parameters, k being small (dots denote time-derivatives), and f(., #, t 
is periodic with period 27 in ¢t. Using an approximate solution due to Kryloff and 
Bogoliuboff it is shown that, if f(x, a, t) is a polynomial in a, a, cos ¢, and sin ¢, the 
‘amplitude’ (a?+ #?/w?)) settles down to an asymptotic range of variation of order 
k as t—> «© which depends on the initial conditions, unless w differs only by order k 
from any one of a set of critical rational values. For these values, we shall say that 
‘subharmonic resonance’ occurs. In this case the upper bound for the asymptotic 
range of variation of the amplitude is of order k/|w—w»|, where wp, is the critical 
frequency nearest to w. The extension of the argument to an equation in which 
f(x, 2, t) is not a polynomial in x, #, cos ¢, and sin ¢, but still periodic in ¢, is indicated. 
The asymptotic behaviour of the ‘phase’ tan~'(xw/#) is also considered briefly. 
The reduction of the case of ‘subharmonic resonance’ to an auxiliary differential 
equation of the first order is outlined, and it is shown that if either a,» min, where 
m and n are large relatively prime integers, or if k/|w—q@,)| is small, the results 
obtained in this way are in agreement with those obtained when subharmonic 


resonance does not occur. 


1. NON-LINEAR differential equations of the type 

¢+wx kf (x, a, t), (1) 
where dots indicate differentiation with respect to the time f, w is a con- 
stant, / a small parameter, and f(x, #,t) a periodic function of t, say with 
period 27, have been investigated extensively. They arise in the discussion 
of the forced vibrations of quasi-linear systems of one degree of freedom. 
In particular, two lines of attack have proved fruitful; the first is an 
application of a method due to Poincaré for the determination of periodic 
solutions (1); the second is a method of successive approximations 
developed by Kryloff and Bogoliuboff (2), which is closely related to 
other extensions of the method of ‘variation of parameters’ which have 
also been employed. 

The object of this paper is to examine to what extent the properties of 
exact solutions of (1) can be deduced from the ‘first approximation’ of 
Kryloff and Bogoliuboff, to indicate how their theory can be extended to 
the case when f(x,#,t) is not a polynomial in 2, #, cost, and sinf (as 


[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 3 (1950)] 
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sumed by them), and to consider in what way their method fails when 
» is near to certain critical values for which ‘subharmonic resonance’ 
curs. It should be noted that an equation such as 

£+w§+kF(é,&) = Esint (2) 
the more usual form of the equation for the forced vibrations of a quasi- 
inear system) can be reduced to the form (1) by putting 


EE, 
x & +-— sin f, (3) 
w-— | 
and there is an analogous substitution if the right-hand side of (2) is any 
neriodie function of t. 


2. We begin by replacing (1) by two simultaneous equations of the 
first order. Let x y sin 6, ee wr cos @. (4) 
Then it is easily shown that 

: l ‘ Fuk. ae 
} (xv-+ w*x)eos 0, G Ww (7+ w*ax)sin 8, 


W wr 


so that (1) is equivalent to 


dr Bs : . ! 
f(r sin 6, wr cos 6, t)eos 6 kF(r, 6, t) 
dt Ww ~ 
(5) 
dé 


w f f(r sin 6, wr cos 8, t)sin 8 wtkGi(r, 6, t) 
dt wr 
It may be noted at this stage that the considerations which will be put 
forward apply to any system of the form (5) provided that F and G are 
periodic functions of @ and ¢, with period 27 in both variables. It is 
bvious that F(r,6,t) and rG(r,@,t) satisfy the same conditions of con- 
tinuity, boundedness, and differentiability with respect to r, 0, and ¢ as 
(7,¢,t) does with respect to x, #, and t. 
A question of some importance in applications is how any solution 
r(t), 6 = @(t) of (5), which takes given initial values 7), 0, when t = 0, 
behaves asymptotically as ¢-> 00. We shall be mainly concerned with this 


juestion. Let us first consider the behaviour of 7, defined by 


x?\} 
(x*4 )* (4a) 
@w* 


which reduces to the amplitude of the simple harmonic vibration satisfying 
|) when & = 0. A simple method of some generality for attacking this 
problem has been given elsewhere (3); here we shall follow Kryloff and 
Bogoliuboff, adopting a slightly different point of view. 
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3. Consider a function p defined by 
p = r—ku(r, 6, t), (6) 
where u is periodic in @ with period 27 and remains bounded as t¢ > 00 for 
fixed r and 6. Then 
dp ™ (F-S— wa) mle +03). (7) 
dt et 06 or 
As the sign of dp/dt is essentially given by the term of order k on the 
right-hand side, when & is sufficiently small, an obvious device is to 
choose uw in such a way that this term assumes a suitable simple form. 


Now as 
zo ; 
lim ; pI (," + +0) dtdé 
Tox 277 ot 00 
00 
Qn T 
= lim - (I u(9, T’)—u(0, 0)| d6+w } [ w(2z,, t)- w(0,t)] a 
T-« 277 ‘ 
= lim # [u(O, T')—u(0, 0)| dé = 0, 
Tox 2QrT7 ‘ 


using the periodicity and boundedness of u, the mean value of the 
coefficient of k in (7) is equal to that of F, i.e. to 

2m 2m 

| f(rsin 6, wr cos 8, t)cos 6 dédt. 


Fir) = sf fr 6,1) dédt = 


el “| 





Tw a 
0 0 (8) 
We may therefore try to choose a function u for which this coefficient 
reduces to /j(7), or differs from F,(7) by less than any pre-assigned « > 0. 
Let us assume, in the first place, that f(x,#,t) is a polynomial in 2, 4, 
cost, and sint, so that F(r,@,t) (and likewise rG(r, 0, t)) is a trigonometric 
polynomial in @ and ¢ whose coefficients are polynomials in r. Thus 
F(r,0,t) = Fo(r)-+ >’ Fup etn o+m, (9) 
where the summation extends over a finite set of pairs of integers (m,n), 
the pair n = m = 0 being excluded. 


esa bm Fim(") eil(nO+ml). (10) 
u(nw-+-m) 


Then, if we can take 


the equation (7) becomes 


: > = kRi(r) ae “i @s } 
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and this can be written in the form 


(6) dp : — 
kFy(p)-+F(r, 0,t, k, w), (11) 
© for dt 
where F,(r, 0,t,k,w) = —[F,(p)—R,(r)|-F— —e = (11a) 
k or fale] 


(7) 
and F, remains continuous and bounded as k > 0. Now this choice of u 


1 the is obviously impossible unless, for all pairs of integers (n,m) which occur 
is t he summation in (9), 

“ sae (9) nw+m ~ 0 (12) 
form. 


and we must at present assume that this holds. This condition is not 
satisfied, and (10) becomes meaningless when w = —m/n, where (n,m) is a 
pair of integers for which F 


I ny 


(r) is not identically zero, and these ‘critical’ 
values of w must be treated separately. 
4. The advantage of choosing uw in this manner is that the asymptotic 
of behaviour of p(t) can be deduced from (11) without difficulty, at least with 
an error of order k, and k is assumed to be small, Let 
M M(R’, R’ ky, w) max |F,(r, 0, t, k, w) 
(RP <r<R,0<0<27,0<t<27,0<k<k,), (13) 
f the where R’ 0, and ky is some fixed (small) upper bound for k. (We cannot 
take R’ = 0 as F, depends on G, and by (5) G@ tends to infinity when r 
tends to zero). We may note at once that M — o if w approaches one of 
Iodt. the ‘critical’ values just referred to. More exactly, if 
5 d(w) min |nw+m (14) 


F cient where (n,m) runs through all the pairs of integers which occur in the 
summation in (9), then we can obviously take 


b 2, &, C(R’, Rk C . 

M sal Piast. (15) 
metric 0(w) fe) 
If w is kept fixed, and it is assumed that 6(w) > 0, then this can be 
ignored for the present. Now it follows at once from (11) that dp/dt and 


(m,n F\(p) have the same sign provided that 
Fp) _ kM. (16) 
To discuss the consequences of this result (which is the key to the 
(10 asymptotic behaviour of p(t) and hence also of r(t)) we shall begin by 
considering a simple, but typical, case. Suppose that /j(p) has a single 
positive zero p a, such that 
Kip) > 0 (0< p <a), Fi(p) <9 (p >a), Fy(a) < 0. 


(17) 
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Then the inequality (16) will hold in intervals (R’, a’) and (a”, R”) such that 
a<a<a'’, F(a’) — F(a") = kM (18a) 

and a”—a’ = O(k). (18b) 

Then clearly 

dp i.» 

P=0 iv <a < #"). 

dt 

Hence if, at any time ¢ = fy, a’ < p(t)) < a’, then also a’ < p(t) < a" for 

all ¢ > ¢,. For p(t) can only leave the interval (a’,a”) either by decreasing 

below a’ or by increasing above a”; but as it is an increasing function of { 

whenever it is less than a’, and a decreasing function of ¢ whenever it 

exceeds a”, this is obviously impossible. 

Again, if initially p(0) = p, is, for example, in (R’,a’), then p(t) must 
increase with ¢ as long as p < a’; hence either p(t) a’ as t — ©, or, for 
some t = ft, p(t) = a’, and in that case a’ < p <a” for all t > to, as has 
just been proved. Similarly if py is in (a”, R") then either p(t) —>a” as 
t—o,ora’ < p <a’ from some finite ¢ = tf), and onwards. Thus asymp- 
totically p(t) varies only in the interval (a’,a") whose length is of order k. 
We may express this by saying that, given any e > 0 there is a f, such that 
a’—e < p(t)<a"+e (t >t) (19) 


(provided that R’ < py < R’) or, more concisely, by 


a’ < lim p(t) < lim p(t) < a’. (20) 
ta oe 
In terms of r and @ this result becomes, by (6) and (10), 
, Bale ‘ 
a —e < r k S am ) pi(nb+mt) F a dae (t ’ ty). (21) 


aw 1(nw+m) 

For small & this result is obviously of considerable practical interest, as 
it gives the asymptotic range of amplitude variation with only a small 
error. It must be noted, however, that the ¢f, in (19) is a function of k, 
and tends to infinity as k + 0; this is to be expected, as for k = 0 the 


solution of (1) is = Alelel-+o), 


where A, a are arbitrary constants, and so p(t) = A is constant and 
arbitrary, and in no way related to the zero p = a of Fi(p). Thus the 
smaller the departure from linearity, the greater the time required for the 
system to ‘settle’ to its asymptotic state. 
We must now consider how this result depends on w. From (15) it is 
obvious at once that (16) must be replaced by 
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that nd henee the breadth of the interval of asymptotic variation, as deter- 
(18 a) mined by the method employed here, is given by 
I8b k 
a” —a’ o ). (22) 
O(w) 
| The result, (21), or (20), is unaltered, but if it is to be of any practical use, 
\(w) must be small. Hence, in particular, it is of little value if w varies 
n” for with / in such a way that 
‘asing m. ‘ 
Mane WwW t nk, (23) 
n of { n 
ver it where again (”,m) is a pair of integers for which F,,,(r) is not identically 
zero, and 7 is a constant of order unity. For this type of ‘critical range’ 
must fw a different approach is oby iously required. 
yr, for 
~ 5. In order to complete the discussion of our special example (17) we 
a os ust also consider the choice of the lower limit R’. Since M depends 
symp nearly on G, and G = O(1/r), it is clear that M = O(1/R’). Hence, if the 
der I nequality (16) is to hold for small values of 7 and so also of R’, we require 
hthat | ‘at kM 
F(p) ms (24) 
(19 R 


ynere M, is a constant, and by (6), 
(20 p R’ kM, (25) 


vhere ./, is the upper bound of |w!, uw being defined by (10). If (0) + 0, 


24) is roughly equivalent to 
(2 R kM, /F,(9), (26) 


0 that our conclusions concerning the asymptotic behaviour of p hold 


rest, as ee : ‘ : : 
7 lown to initial values of r of order k. If, as is more usually the case, 
» sma apa : ‘ . : 
cz. | A(0) = 0, some modification is required. Suppose, for example, that 
n OF &, , . 
0 4 F,(0) 4 0, so that (24) is approximately 
ie 
F’(0)p > kM,/R’. 
t ombining this with (25) it follows that we may take R’ to be of order k! 
nt and : ; arm : , 3 
an de B® this case; similarly, if F,(p) ~ p# for small p, R’ can be taken to be 
lus tne ETE . . ° , 
for th forder k10+, For initial values of r in the interval (0, R’) there may 
or the ‘ ; 
solutions which remain small for all t > 0; these are unstable in the 
15) it i sense that if the system is disturbed it will in general drift away from 
») IU : y 


this state of motion and settle down to a motion whose amplitude is in 
1,a"). These unstable, permanently small, solutions, if periodic, have the 
period 27 of the driving force, i.e. they correspond to the ‘fundamental’ 
is distinct from the ‘subharmonic’ response of the system. 


092.11 


Bb 
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6. It is easy to extend the results obtained to the general case, when 
the polynomial /j(p) has several zeros. Suppose for instance that F,(p) 
has ” positive simple zeros p = 44, dy,..., @,, 80 that /(p) is alternately 
positive and negative in the intervals (a;, a;,,). Let a; be a zero of F,(p) 
for which 
Fi(p) > 9 (aj. < p <a;), Fy(p) < 9 (a; < p <@j44), Fo(a;) <0. (27) 
For sufficiently small /, the inequalities 
Fy(p)| > kM, Fi(p)| < kM 
determine three intervals (a4_,,a;), (a, a), (aj,@;,,) such that 
F(a; 1) F,(a;) — F(a) F(a‘, 1) = kM 
Pip)! >kM (a4 <p <a, af <p <aiy) 
F(p)|} << kM (a; < p< a’) 
@;-4—;_4 O(k), aj—a, = O(k),  aj44—j44 = OC) 


u u 


(28) 


It then follows, as in the special case which has just been discussed, that if 


Qi-1 < Po < Gry (29a) 
then a’, < lim p(t) < lim p(t) < aj. (29b) 

to rae 
Thus the initial values of (r,@) can be classified. If we interpret (7,6) as 
polar coordinates in a plane, then the curves 

r—ku(r,@,0) = p 


are a congruence of closed curves. To every zero a; for which (27) holds 
corresponds a domain .%, defined by 


a;_, <r—ku(r, 0,0) < aj... (30a) 


t 
with the property that the asymptotic values of a solution whose initial 
values lie in A, are in the closed domain 7; defined by 


’ ” 


a, < r—ku(r, 0,t) < a, (30b) 


u 


and this is an annulus with width of order k, so that we may also write 
lim |r —ku/(r, 0,t)—a,; O(k). (31) 
t= 
The solutions whose initial values lie in ZY; possess a type of stability 
analogous to ‘orbital stability’ which may, perhaps, be called ‘amplitude 
stability’. In the case envisaged here, the domains .9/, are separated by 
annular domains of width O(/) which are completely unstable in the sense 
that, though certain solutions may remain in them for all ¢, if such a 
solution is disturbed sufficiently (though only to order k) it drifts away 
from the domain and settles down in one of the ‘stable’ annuli Z;. The 
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neighbourhood of the origin must be treated separately. The ‘unstable’ 
case has already been discussed, and the extension to the ‘stable’ case 
is obvious. 
7. A similar argument can be developed for the coordinate 6. Let 
d = 6—kv(r, 0, t) (32) 
ind suppose that rG, like F, is ‘a trigonometric polynomial in @ and ¢ 
whose coefficients are polynomials in r. Then 
G = G(r)+ DY’ Gam(ryetre+™, (33) 
If we assume once more that m-+-nw + 0 for all pairs (n,m) which occur 


n this summation, then we may put 


ee ‘ 
Vir. 6.t) N\ ; at) ei(nO+mt) (34) 
an U(NW+M) 
. dd 
und obtain w+kGy(r)+ O(k), 

dt 

dd 
T= w+kGy(p)+O(0). (35) 


Now as, if we consider solutions whose initial values lie in °%,, |p—a;| 


t t 


S asymptotically of order k, it follows that 


4 db | a 
lim |* ee kG(a;)| = O(k?). (36) 
‘— 

If, in particular, the initial values lie in Y%,, this can be replaced by 


} 


@ bd 


0 wy kG(a,)| = O(k). (37) 


From this inequality, a necessary condition for the existence of sub- 
harmonic solutions of the original equation (1) can be derived. Such a 
solution is by definition periodic with least period 2n7, where n is an 
integer greater than one (any periodic solution of (1) has necessarily a 
period which is 27 or an integral multiple of 27). The corresponding p(t) 
also has period 2nz, while 4(¢) is increased by 2mz if t is increased by 2nz, 
m being another integer. Hence, applying (37) with ¢ = 2nz, it follows that 


w+kG(a,)+0(k) = ™, (38) 
n 
so that subharmonics can only exist for given m and n and arbitrarily 
small k, if w differs from the rational number m/n by O(k). (In particular, 
this includes the ‘critical’ ranges of w for which our present method 
breaks down.) Alternatively, x may be of order 1/k or larger. 
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If no subharmonics exist, one would expect that among the solutions | (Wh 
whose initial values lie in 7; there are almost periodic ‘drifts’ of the form 
x fa, +k R(t, «+ot)sinfot+a+kS(t,x«+-ot)}, (39) 

where R, S have period 27 in both arguments, o is an irrational number | We 


depending on a,;, and « is an arbitrary constant (4). Now if p and ¢ are tha 
introduced as dependent variables, then it is found that they satisfy 





differential equations of the form 


dp 
-F 2 i’ . 
dt KF y(p)+k*F(p, , t, k) | whe 
(40) | de 
p 
dd , on 
w-+ kG(p) f k?G',(p, d, t, k) steg 
dt 
- ° ° ° ° ° ° “= 
where F,, G, are periodic in both ¢ and ¢ with period 27, and remain a 
continuous and bounded as k + 0. Starting with these equations, it can but 
be shown that every solution whose initial values lie in 7; tends asymptoti and 
cally to a solution of the form (39) or to a subharmonic solution; the proof, the 
which is lengthy and depends on ideas which have not been introduced obt 
e ) 

in this paper, will be omitted.+ It may be added that if the terms of order |, 
k? are omitted in (40), the equations of the ‘first crude approximation’ of | yp, 

Kryloff and Bogoliuboff (2) are obtained, to whom the equations (40) are 
also due. These approximate equations give an estimate of the manner in ’ 
which an arbitrary solution approaches its asymptotic form, and many | |: 
Fon 


examples of their application will be found in reference (2). 
8. If the function f(x,#,¢) is not a polynomial in 2, #, cost, and sint, 
the method developed here seems to break down. For then the finite sum 


in (9) must be replaced by a Fourier series, and even if w is irrational, so | wh 


that (12) is satisfied for all pairs of integers (n,m), the convergence of car 
(10) is not a necessary consequence of that of (9), and its discussion isa | py 
matter of some difficulty. But, by an extension of Weierstrass’s theorem, | tric 
we can with any preassigned « > 0, find an F*(r,6,t) which is a trigono- anc 
metric polynomial in @ and t with coefficients which are polynomials inr | obj 
such that Fir, 0,t)—F*(r,0,t)\ <e« (R'’ <r< R’), (41) | the 
the order of the trigonometric polynomial F* tending to infinity as « > 0. 
Now if we define p* by p* = r—ku*, (42) 
where u* is the trigonometric polynomial satisfying 

= al F* — Fs (r) (43) 

ot fal] 

+ A result of this type seems to have been proved by Kryloff and Bogoliuboft (5), an 

presumably with f(x, @,t) a polynomial in x, #, cost, and sint. A variant of this proof, wh 
derived independently, in which this restrictive assumption is dispensed with, will, it 1s Sir 


hoped, be published elsewhere. 
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where F*(7) is the mean value of F*), it follows that 





itions 
- form dp* , cu* ou* 
: kF*(r)+k( F—F*)—k(F LG , (44) 
(39) | dt or 06 
imber We can now repeat our arguments with F% instead of F,, and it is found 
p are | that results such as (20) or (21) hold except that 
atisty I: 
a” —a’ Of --el, (45) 
\o* 
where 5* is defined in relation to F* as 5 was defined for F. Thus 8* 
(40 lepends on the order of the polynomial F*, and in general it decreases 
steadily to zero as this order tends to infinity, which happens if we let 
-0. The set of ‘critical values’ also depends on e, and becomes dense 
on 0,00) as e > O unless F is of a special form. Hence we cannot let « > 0; 
+ pay . - ° b 
v can but there will be a value of « for which ¢+4/8* is as small as possible, 
iptoul ind this value will give the best possible bounds. It is not easy to discuss 


prool the details of this process in general, but the character of the results 


duced btained here persists, though with a certain loss of precision. (This 
order result may be compared with that obtained by a different method else- 
on Of | where (3).) 
10) are 
—= 9, Returning to the simpler case. of an f(x,#,t) which is a polynomial 
oe na, #, cost, and sint, we must now consider the ‘critical’ ranges of w. 
For this purpose, let us now assume that 
int ¢ ' 
= Ww f T nk, (46) 
fe sum p 
nal,so | where p, q are integers and » is a parameter of order unity. This case 
7. ot . . . . . : : 
nce ol in be dealt with by using the methods originated by van der Pol and 
m is a by Appleton, and developed by Cartwright and Littlewood (6, 7) for the 
eorem, | triode oscillator synchronization problem. We shall consider it briefly, 
‘1gONno nd show that if either p and q are large, or if 7 is large, the results 
uls in btained will approximate to those which have just been found from 
(4) the first approximation of Kryloff and Bogoliuboff. 
ie =e, Ou equations (5) now become 
/ ) di ’ , ’ 
(42 | ki Fy ) NN F (ie (nd — | 
at sa ~ 
(47) 
dé wa i(nO-+mt) 
(43 = (GIP)FE N+ G(T)+ 2D Gam(reer rm} 
' if si 
off (5), ind these can be considered as defining the motion of a point in a plane 
is proo!, | whose polar coordinates are (r,@) (the phase plane of the equations). 
vill, it 


‘ince the right-hand sides of these equations have period of 27 in ¢, it is 
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sufficient to know the trajectories starting at initial values (75, 4,) in the 
time interval (0,27) in order to deduce the trajectories in the complete 
semi-infinite time interval (0,00). For by the uniqueness theorem relating 
to the solutions of differential equations, if the trajectories are given by 


r= R(ry,9,t), 0 = O(ro, 9%, t) (0 <t< 2zn), (48) | 


then 
r= R(r,,0,,t—27), 0 = O(r,,9,,t—27) (27 - 


rs 
= 


and generally 


r R(r,,0,,t—2s7), 0 = O(r,, 0,,t—2s7) (2s7 < t < 2(s+-1)z), 


where 
P r(2s7) Rir,_.,0,_,, 2), 0. 6(2s7) = O(r,_1, 9,-1, 277). 
Instead of using the interval (0, 27) in this way, we can use the interval 
(0, 2p7). The trajectories will therefore be completely determined if we 
can find the points (7,,,, 9,,,) (8 = 1, 2,:..) corresponding to the point (75, )). 


Now O(r,.,, 9... 2p). (49) 


sp? “sp? 


Tie1)p = R(rep, 5p, 20), 0 


(s+Dp 
We are therefore led to study the mapping of the phase plane into itself 


defined by y” Rr, 8, 2p7), 0’ O(r, 8, 2p7). (50) 


Moreover, we can subtract any integral multiple of 27 from 6’, as (7,4) 
and (7,@+-27) define the same point. Now it follows from (47) that 


2p7r 
rp = Tot k | Fro 99+ Mat) dt+O(kp?), 
) 
é f 
2pT | 
6, = 4-44 2qr- "1 2pmn+ | G(r 6, a 11.) at + O(k?p?). 
; p 
0 


Omitting the term 2q7 on the right-hand side of the last equation, and | 
using the expansions of F,, @ it follows that 


r’ r- 2rpk) Fy(r) 4 > F,,, iq" )é v9! +. O(k2p?) | 
; (51 
g(r err? +0(k2p*)| 


dy 


0 = 04 2xpk\y+G,(r)- x &, 


10. To discuss the mapping (51) for small k, we need the result that, 
if p and & are the solutions of the differential equations 


d : by ; 
E = klBy(p)+ > Rip,-alp)e**| | 


dt a si 
dis ee oe a 
dt kins Go(p) r p Grp, ig( PE ins 
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r(2spm)— p(2spr) = O(kp)) 


then " 
| O(2sprr) —ys(2spr)] (mod 27) - O(kp)) 


The proof of this theorem will be omitted here. Then we can determine 
the points (r,,,9,,), and hence the complete trajectories, with an error of 
order k, if we can solve the equations (52). This is a simpler problem than 
the original one, as the right-hand sides of (52) do not depend on t. The 
asymptotic behaviour of the solutions of (52) is fairly simple. They either 
tend to a ‘critical point’, i.e. to a point at which both dp/dt and dy/dt 
vanish, or they approach a closed trajectory asymptotically. The first 
case corresponds to a subharmonic of order p of the original equation (7), 
while the second one corresponds to a ‘drift’ solution of the type (39) (this 
result is due to Cartwright and Littlewood, but unpublished). Obviously 
each case must be considered on its merits, but there are two generally 
valid results which must be mentioned. 

In the first place, if we suppose that the series expansions in (47) are not 
( 


as n +|m|—> oo, by the Riemann—Lebesgue lemma. Hence if both p and 


finite, but infinite double Fourier series, then the F,,,, G,,, tend to zero 


qare large, and pk small, all the terms in (52) except F, and 7+ G, will be 
negligible, and these equations reduce approximately to a form satisfied 
by the first approximation of Kryloff and Bogoliuboff (with 4 = ¢+-a1). 

Secondly, if » is large, but ky still small, and if we consider the single 
first-order equation equivalent to (52), 


dp Fy(p)+ > Fp,-1(p)exp(ilpy) 
dp 4 +Golp)+ & Giy-rqlp)exp(ilpp) 


then (for sufficiently large ») the only critical point will be the origin, 


(54) 


which corresponds to a fixed point of the mapping 
yr’ R(r, 6, 27), 0’ ©(r, 0, 27) 


and hence to a solution of the original equation of period 27. Also, there 
will be closed trajectories within O(1/n) of the circles 


p = a, F(a;) = 0. (55) 


This is in complete agreement with the results already obtained in the 
non-critical case. To prove it, we write (54) in the form 
dp ] . - . P ] — 
(Fy(p)+ > Fp a(ee@?}4 o( 3): (56) 
dis n n* 
We are assuming that is so large that the denominator of the right-hand 
side of (54), which is proportional to the right-hand side of the second 
equation of (52) is positive, so that the trajectories encircle the origin in 
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the clockwise sense. Thus if p = py for y = 0, then p’ = p(27), where p 


is considered as a function of %, is given by 


p Po {Fo(p) i > Fry, 1q\ Po) pp} dys I ol] 
0 
: 27 ay l 
or p Po 7 Fy( po) 4 of.) (57) 


Hence the closed trajectories of (56), for which p’ = pp, are within O(1/7) 
of the zeros of Fj(p). Moreover, it can be shown from (57) that a closed 
trajectory associated with a zero p = a of F,(p) is stable if F)(a) < 0, and 
completely unstable if F(0) > 0 (if Fo(a) = 0 the trajectory may be semi- 
stable, and in any case the interpretation in terms of the exact solution 
of the mapping equations (51) becomes difficult). The corresponding exact 
solutions of the (7, #)-equations are orbitally stable, or completely unstable 
‘drifts’, respectively. Thus we have, for a stable drift associated with a 
zero p = a of Fi(p), that 


r—a,| = O(k)4 0") Ofk l ) o(5). 
n) fo) fo) 


where 5 nw—m 


and this is obviously in agreement with our earlier results. 
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{ NOTE ON THE DIFFRACTION OF A 
CYLINDRICAL WAVE BY A PERFECTLY 
CONDUCTING HALF-PLANE 
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(Electrical Engineering Department, Imperial College, London) 
Received 7 March 1950) 


SUMMARY 


The most easily derived solutions in diffraction theory are for the case of incident 
ine waves. The solution for a line-source can then be obtained by a further 
gration, but this method, although well known in some fields, does not seem to 
been applied in diffraction theory. Because of its importance in more compli- 
ited problems it is here illustrated by the example of a line-source in the presence 
in infinitely thin, perfectly conducting half-plane. The introduction of Hankel 
netions is avoided, and the solution appears in a useful form, equivalent to that 
given by Macdonald and analogous to Sommerfeld’s famous solution for an incident 
ne Wave 


1. Introduction 

THE problem to be considered is the two-dimensional one illustrated in 
Fig. 1; T is a line-source, situated at (79, 4)), propagating the cylindrical 
wave (omitting the time factor) 


Lre-im4H2)(bR) w 
Nae) ) 0 ( / (AR) ( ) 











[Quart. Journ. Mech. and Applied Math., Vol. III, Pt. 3 (1950)] 
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and the half-plane y = 0, x > 0 is occupied by an infinitely thin, perfectly 
conducting sheet. The problem is essentially scalar, being expressed in 
terms of EL, or H, according to the polarization. For convenience, attention 
will be confined to the case where H = (0, 0, H,). 

Carslaw (1, 2) gave the solution in a form equivalent to 


H, = u(8,)+u(—9), (2) 
where 
u(0,) ein/4 s sin $a Hk fp21 72 __ Opp cos(a—6 Nd 
. 4,/(27) J cos 4a—cos}(6—0,) °* YS 5 9 “OT” x—6o)§] dar, 


(3) 
the integration being over a path in the upper half of the complex plane, 
going from 7+i00 under the branch-point to —72-+io0. 

In his first paper Carslaw (1) prefaced the derivation of his solution with 
the words ‘... the results .. . are not obtained in a workable form owing 
to the necessity for introducing Bessel’s functions into the integrals . . .’, 
But later Macdonald (3) showed that (2) could be usefully transformed into 


. Eo 1 
eit/4 a ; et7/4 i , 
H. p ikR ch € dé ate e ikS ch & dé, (4) 
2 4 (277) J 4 (277) a 


where Rsh&) = 2,/(797)cos $(@—4)), Ssh&, = 2,/(797)cos $(0+-6,), and sh 
and ch denote the hyperbolic sine and cosine respectively. 

In this note the problem is treated by a different method which does not 
involve Hankel functions. Since the final result is a form essentially 
equivalent to (4), the reason for undertaking the analysis must be briefly 
stated. 

The use of integral equations has recently led to the rigorous solution 
of several new problems in diffraction theory, one of which has important 
applications in radio propagation (4). Because the physical discontinuities 
are planar, the fundamentally easiest problems are those for which the 
incident wave is plane: in the radio propagation problem, however, it is 
essential to deal with a transmitter which is conveniently taken as a line- 
source. Now, although the dictum that general solutions can be built up 
from plane-wave solutions is perhaps more often preached than practised, 
the procedure has proved extremely useful in propagation theory, in 
particular in connexion with the problem just mentioned. This latter 
problem is a rather complicated generalization of the present one; but it 
seems desirable that the results to be given in this note should be set out 
independently, both because of their interest in a different field and 
because it is best to illustrate the technique by applying it to a relatively 
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simple example. Furthermore, the solution (4) is perhaps not as well 
known as it deserves to be, although it is quoted by Bateman (5). Baker 
and Copson (6), for instance, only mention Carslaw’s form (2), and indeed 
the author was not aware of Macdonald’s result until after it had been 
obtained independently by the present method. It will be seen later that 
(4) is virtually as simple as Sommerfeld’s famous solution for an incident 


plane wave. 


2. The method 


First, consider the problem when the incident field is the plane wave 
Hi = etkreos@-a) (0 <a < 2). (5) 
The total field may be conveniently written as 
HP = HPG+ Hp, (6) 


where the affix p indicates that the incident wave is plane, and the 
affixes G and D imply ‘geometrical optics’ and ‘diffraction’ respectively. 
Then Sommerfeld’s solution is equivalent to 
eo ee 

é COS ga COS a8 ikr cos(@—B) dp 


? 
COS P- COS & 


for y> 0, 


77 


sie = 
Hp ——- ' (7) 
: D COS $a COS ¢ eee . 
2 a8 ikr cos(0+B) dp, for y < 0, 
x cos B+-cos « j 
S on g 


where S(@#) denotes the path of ‘steepest descents’ passing through the real 
angle #0, as in Fig. 2. 

Now consider the line-source at (79,9) where 0 < 6, < 7. Introducing 
the appropriate phase factor exp| —ikr, cos(8)—«)], Sommerfeld’s integral 
representation of the Hankel function may be written 


. 


H)(kR) e—tkrecosOy—agikrcos—a) dy (y < rysin 6). (8) 


77 


« 
S(in 


The form of (8) indicates that the solution for an incident field (1) can be 
obtained from that for an incident field (5) on multiplying by the factor 


(27) bp —im/4p —ikry cos(Oy—a) (9) 


and integrating with respect to « over the path S(47). 
A point to be noticed is that « is complex on S(47), whereas it is normally 
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taken to be real. But this introduces no difficulty, for (6) and (7) are still 
valid if it is appreciated that, for 0 < rea < z, 


, pikrcos(@—a)_|_ pikrcos(0+a) for « to the left of S(z A), 


eikr cos(d—a) for « to the right of S(7—8@) 
HpG and S(@—7), (10) 
0 for « to the right of S(7—8@) 


and to the left of S(@—z), 


relations which yield the usual ‘geometrical optics’ field when « is real. 














ee ee ee ee es ee ee ee we es ee oe = 


ic. 2. 


For a fixed 6, therefore, the terms in (6) are continuous for all « on S(1z). 


and the suggested procedure leads to the solution of our problem. 


3. The solution 
The solution that we have to consider is 


im/4 % 

¢ , : 

H HE ikrg cos(O,—a) dx : 
. (27) , V(27) 

S(iz) S(im) (11) 


HP: ikrg cos(Gy—a) d x. 


where H?® is given by (10) and H?” by (7). To separate this field into 
‘geometrical optics’ and ‘diffraction’ terms, the path of integration for a 
in the second expression in (11) must be taken through the ‘predominant’ 
angle @,. The poles of the integrand of (7), regarded as a function of x, are 
given by cos « cos 8, where B lies on S(@#) or S(27—86) according as to 
y is positive or negative: in transposing the a-path of integration from 
S(4) to S(8@) these poles may be crossed, and their residues must then 
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» il] . rr . . ~ . . . . 
still be allowed for. The contribution of these residues combines with the first 
expression in (11) to yield the ‘geometrical optics’ term of the solution. 
In fact. if we write (11) as 
H, H¢+ HP, (12) 
10) it is easy to check, from a consideration of the different cases, that the 
correct form of H®% is obtained, namely 
| v (dre 4 AP (AR)+ HP (RS)| for 0<60<27—h, 
real. H® ) (dare 4 HYP (KR) for 7—0,<0<274+6, (13) 
0 for 7+0, <0 < 27; 
and then, for all values of y, the diffraction term may be written 
( 1 P i COS 1\ 1 COS 1 B T 6g > 
HP | | 2\> 0) 3(/ ) ik(rg COS a+r COs f) dad, (14) 
7m. (277) J J cos(a+6y)+cos(B+6) 
SOO S(O 
the symmetry in 79, 6) and 7, @ being demanded by the reciprocity theorem. 
Equations (13) and (14) give one form of the solution. It is convenient 
to express (14) as i 
expre H! 1(8,) 1 J 4), (15) 
eit j in bg ¢ ik(rg COS A+F COS B) 1 
where I(6,) ( rd B. (16) 
477, 277) b i cos 3(« | B | A, +f) 
SCO SOO 
By taking a+ as a new variable it may be seen that this can lead to 
Carslaw’s form (2). 
4. The reduction of the solution 
Multiplying top and bottom of the integrand in (16) by 
Siz 4 cos }(a—B+0,+8) 
and discarding that part which is an odd function of 8, we have 
Hh ‘ 
87r,/( 277) 2 E | COS 4 x4 GO, 1_@)—sin $B 
S(O SiO 
1 ik(r9 ec ‘cos B) 7 
ly. ’ cos(3B)e SE COS SF COR dadB. (17) 
cos 4(a+6,+6)+sin $8 ‘ 
L] ; . 
In the second of these integrals change « to —a and then recombine the 
d into a 
two. This gives 
} tor Y 
; pin/4 
inant I(0,.) 
v. are <7, (=77) 
an te cos 4(6, #)cos } x COS Be ik(ro cosa+rcos B) d vd B 
. fron (cos a—1)+- (cos B—1)—4 sin $(8,+6)sin 4a sin $8+ 2 cos? $(6,+6) 


(18) 
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Now make the ‘steepest descents’ substitutions 





x V2 e-*7/4 sin 4a, y V2 e-'74 sin 48, (19) 
so that 
ei7/4 faa 
I(8,) ~ amend tkro+7) cos $(0,+-8) . 
Tr (277) 
. . ¢ k(ro x? +ry?) 
x | - ‘ , . dudy. (20) 
a+ y?+ 2 sin $(0,+0)xy+ 21 cos? $(0,+6) 


Let R, = ry +r, and make the polar substitutions 


x = /(R,/r9)p cos ¢, y = ,(R,/ro)psin ¢. (21) 
Then we may write ‘ 
I(#,) = — — ekRs eos 4(0,+0) pd (p)e**iP* dp, (22) 
74) (277) : 
where 
J(p) 7 liele.\ere® oe 4 1 . a 
J ply (7/rp)cos*h+,/(r9/r)sin’d+ 2 sin $(0,+4)sin ¢ cos d |4 
. +-{2t,/(797) R,}cos* 4(6,+) 
(23) 
J(p) can be evaluated by the standard technique of putting z exp(id). 


The path of integration then becomes the unit circle, and it is only 
necessary to calculate the residues of the enclosed poles. Without going 
into details, it will be found that there are four poles altogether, given by 
z +p, 2 = +Ppo, say; but |p, p, 1, so that only two poles lie within 
the unit circle. The result of the integration is 
. . 1 i 

0") cost 1(0,+0)| , (24) 
R? 7 


where the square root is the branch in the first quadrant of the complex 


J (p) 2m |sec 4(6,+ 84) G +. 22 p? - 


plane. 
eferring to the geometry of Fig. 1, (22) evidently becomes 


x 


w KR, p* 


1(6 ) T ein 4 ikR, - : . pé " dp, 
' ds | vile? +2((R,—S)/ Ry) p?+-2(( 2, + S)/R)]} 
0 
F for cos $(@,+6) 20. (25) 
| S 
The substitution p?>+i * s—2# 
R, R, 
then gives : 
1(6,) = [Petnte-ws fe f 1(0,-+-6 
( ee iv ik ' + r cos 4(6,-4 Q), 
>) —¢ ¢ | JRL 2) ‘ or COS (4, ) 











or 





(19) 


(20) 


(21) 


(22) 


0) 
(23) 


p(id). 
only 
going 
en by 


vithin 


(24) 


mplex 


(26 
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The corresponding expression for J(—6,) will be 


9 e e—tkR» 


I(—4) F /—efmhe-he dA, + for cos 4$(0,—@) 2 0. 


. V 
V{UR,— R)/R} (27) 
Now we can write the incident field (1) in the form 


2 : e-ikR¥ 
|= et7/4e—ike | gone Os (28) 
N 7 P. 4/(A?+-2) 


and the ‘reflected’ field is given by the same expression with S replacing R. 
Using (28), the diffraction term (15), determined by (26) and (27), can be 
neatly combined with the geometrical optics term (13). In order not to 
have to specify changes in the signs of the square roots it is convenient 
to write 


q=2 mae } cos 3(4 6) ty {k(R,—R)}, + for cos $(6)—8) : | 


“AN \R,+R 
{ kro . ~ 
= 2 r : g} eos 4(0,+6) ty {k(R,—S)}, + for cos $(@,+ 84) = 0. 
N “TA 
(29) 
The total field then becomes 
2 . : id? ; < ¢ iA? 
H. ei/4| o-ikR dA\+e-iks dA}. (30) 
A) 7 J (A?+2kR) J /(A?+2k8) 


€ €2 


5. The nature of the solution 
[fin the first and second terms of (30) we make the respective substitutions 
A (2k R)sh($€), A (2k R)ch(3€) 


the solution (30) goes over to Macdonald’s form (4). For the purpose of 
computation, however, (30) is the most satisfactory form. It is completely 
analogous to that well known for an incident plane wave; and to see how 
30) goes over to this latter solution when the line-source is taken to 
infinity we need only multiply by ,/(krg)exp(ikr,) and let ry 0. This 
gives the famous result 


m4 
(etkreos8—Co) FT —,/(2kr)cos $(8,—8@) |+ e%reos9+0) FT —,/(2kr)cos $(0,+)]}, 
\7r : > 
(31) 
where F(a) ( e-i* dy. (32) 


a 
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The diffraction term (27) may be written 


e—tkR, Ey 


dx: 33 
J (A?+2R/R,) (33) 


v(1— RII) 

and providing only that /R, is large compared to unity, an asymptotic 
approximation to (33) is obtained by putting A equal to its lower limit 
value in the non-exponential factor of the integrand (7). This gives 


¢ ikR 


| = ei7 F[./{k( R,— R)}}. 34 
Na \{k(R,+ R)} I, , | (34) 


The approximation to (26) under the same conditions is obtained on 
replacing R by S in (34). These approximations give the complete field 
in terms of the tabulated Fresnel integral, as for the case of an incident 
plane wave, and they are only invalid if both the source and the point of 
observation are well within a wave-length of the diffracting edge. 
If k( R,—R) is large, the asymptotic expansion 
1 : . 
_ I 
ela k (a) ~— : 

21a 

enables (34) to be written in the form 
e-in/d y ikro y ikr 
9 /(9 :1(A ) kr kr) (35) 
2, (277)cos $(@9—@) 4 (Arg) y (AT) 
This and the analogous expression with —@, replacing 6, represent the 
diffraction field as a cylindrical wave diverging from the diffracting edge. 
It is interesting to note that the representation (35) is valid everywhere 
except in the region close to the shadow edge 6—@, = 7 bounded by a 
hyperbola 
yI R—r = constant. 
This hyperbola is the counterpart of the parabola well known in the case 
I 

of an incident plane wave (see e.g. 6, p. 145). 
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